Integrated methylome and phenome study of the circulating proteome reveals markers pertinent to brain health

Characterising associations between the methylome, proteome and phenome may provide insight into biological pathways governing brain health. Here, we report an integrated DNA methylation and phenotypic study of the circulating proteome in relation to brain health. Methylome-wide association studies of 4058 plasma proteins are performed (N = 774), identifying 2928 CpG-protein associations after adjustment for multiple testing. These are independent of known genetic protein quantitative trait loci (pQTLs) and common lifestyle effects. Phenome-wide association studies of each protein are then performed in relation to 15 neurological traits (N = 1,065), identifying 405 associations between the levels of 191 proteins and cognitive scores, brain imaging measures or APOE e4 status. We uncover 35 previously unreported DNA methylation signatures for 17 protein markers of brain health. The epigenetic and proteomic markers we identify are pertinent to understanding and stratifying brain health.

Complex trait methylation scores in the prediction of major depressive disorder

BACKGROUND

DNA methylation (DNAm) is associated with time-varying environmental factors that contribute to major depressive disorder (MDD) risk. We sought to test whether DNAm signatures of lifestyle and biochemical factors were associated with MDD to reveal dynamic biomarkers of MDD risk that may be amenable to lifestyle interventions.

METHODS

Here, we calculated methylation scores (MS) at multiple p-value thresholds for lifestyle (BMI, smoking, alcohol consumption, and educational attainment) and biochemical (high-density lipoprotein (HDL) and total cholesterol) factors in Generation Scotland (GS) (N=9,502) and in a replication cohort (ALSPAC_adults, N=565), using CpG sites reported in previous well-powered methylome-wide association studies. We also compared their predictive accuracy for MDD to a MDD MS in an independent GS sub-sample (N=4,432).

FINDINGS

Each trait MS was significantly associated with its corresponding phenotype in GS (β_range=0.089-1.457) and in ALSPAC (β_range=0.078-2.533). Each MS was also significantly associated with MDD before and after adjustment for its corresponding phenotype in GS (β_range=0.053-0.145). After accounting for relevant lifestyle factors, MS for educational attainment (β=0.094) and alcohol consumption (MS_p-value<0.01-0.5; β_range=-0.069-0.083) remained significantly associated with MDD in GS. Smoking (AUC=0.569) and educational attainment (AUC=0.585) MSs could discriminate MDD from controls better than the MDD MS (AUC=0.553) in the independent GS sub-sample. Analyses implicating MDD did not replicate across ALSPAC, although the direction of effect was consistent for all traits when adjusting for the MS corresponding phenotypes.

INTERPRETATION

We showed that lifestyle and biochemical MS were associated with MDD before and after adjustment for their corresponding phenotypes (p_nominal<0.05), but not when smoking, alcohol consumption, and BMI were also included as covariates. MDD results did not replicate in the smaller, female-only independent ALSPAC cohort (N_ALSPAC=565; N_GS=9,502), potentially due to demographic differences or low statistical power, but effect sizes were consistent with the direction reported in GS. DNAm scores for modifiable MDD risk factors may contribute to disease vulnerability and, in some cases, explain additional variance to their observed phenotypes.

FUNDING

Wellcome Trust.

DNA methylome-wide association study of genetic risk for depression implicates antigen processing and immune responses

BACKGROUND

Depression is a disabling and highly prevalent condition where genetic and epigenetic, such as DNA methylation (DNAm), differences contribute to disease risk. DNA methylation is influenced by genetic variation but the association between polygenic risk of depression and DNA methylation is unknown.

METHODS

We investigated the association between polygenic risk scores (PRS) for depression and DNAm by conducting a methylome-wide association study (MWAS) in Generation Scotland (N = 8898, mean age = 49.8 years) with replication in the Lothian Birth Cohorts of 1921 and 1936 and adults in the Avon Longitudinal Study of Parents and Children (ALSPAC) (Ncombined = 2049, mean age = 79.1, 69.6 and 47.2 years, respectively). We also conducted a replication MWAS in the ALSPAC children (N = 423, mean age = 17.1 years). Gene ontology analysis was conducted for the cytosine-guanine dinucleotide (CpG) probes significantly associated with depression PRS, followed by Mendelian randomisation (MR) analysis to infer the causal relationship between depression and DNAm.

RESULTS

Widespread associations (NCpG = 71, pBonferroni < 0.05, p < 6.3 × 10-8) were found between PRS constructed using genetic risk variants for depression and DNAm in CpG probes that localised to genes involved in immune responses and neural development. The effect sizes for the significant associations were highly correlated between the discovery and replication samples in adults (r = 0.79) and in adolescents (r = 0.82). Gene Ontology analysis showed that significant CpG probes are enriched in immunological processes in the human leukocyte antigen system. Additional MWAS was conducted for each lead genetic risk variant. Over 47.9% of the independent genetic risk variants included in the PRS showed associations with DNAm in CpG probes located in both the same (cis) and distal (trans) locations to the genetic loci (pBonferroni < 0.045). Subsequent MR analysis showed that there are a greater number of causal effects found from DNAm to depression than vice versa (DNAm to depression: pFDR ranged from 0.024 to 7.45 × 10-30; depression to DNAm: pFDR ranged from 0.028 to 0.003).

CONCLUSIONS

PRS for depression, especially those constructed from genome-wide significant genetic risk variants, showed methylome-wide differences associated with immune responses. Findings from MR analysis provided evidence for causal effect of DNAm to depression.

Epigenome-wide association study of alcohol consumption in N = 8161 individuals and relevance to alcohol use disorder pathophysiology: identification of the cystine/glutamate transporter SLC7A11 as a top target

Alcohol misuse is common in many societies worldwide and is associated with extensive morbidity and mortality, often leading to alcohol use disorders (AUD) and alcohol-related end-organ damage. The underlying mechanisms contributing to the development of AUD are largely unknown; however, growing evidence suggests that alcohol consumption is strongly associated with alterations in DNA methylation. Identification of alcohol-associated methylomic variation might provide novel insights into pathophysiology and novel treatment targets for AUD. Here we performed the largest single-cohort epigenome-wide association study (EWAS) of alcohol consumption to date (N = 8161) and cross-validated findings in AUD populations with relevant endophenotypes, as well as alcohol-related animal models. Results showed 2504 CpGs significantly associated with alcohol consumption (Bonferroni p value < 6.8 × 10-8) with the five leading probes located in SLC7A11 (p = 7.75 × 10-108), JDP2 (p = 1.44 × 10-56), GAS5 (p = 2.71 × 10-47), TRA2B (p = 3.54 × 10-42), and SLC43A1 (p = 1.18 × 10-40). Genes annotated to associated CpG sites are implicated in liver and brain function, the cellular response to alcohol and alcohol-associated diseases, including hypertension and Alzheimer's disease. Two-sample Mendelian randomization confirmed the causal relationship of consumption on AUD risk (inverse variance weighted (IVW) p = 5.37 × 10-09). A methylation-based predictor of alcohol consumption was able to discriminate AUD cases in two independent cohorts (p = 6.32 × 10-38 and p = 5.41 × 10-14). The top EWAS probe cg06690548, located in the cystine/glutamate transporter SLC7A11, was replicated in an independent cohort of AUD and control participants (N = 615) and showed strong hypomethylation in AUD (p < 10-17). Decreased CpG methylation at this probe was consistently associated with clinical measures including increased heavy drinking days (p < 10-4), increased liver function enzymes (GGT (p = 1.03 × 10-21), ALT (p = 1.29 × 10-6), and AST (p = 1.97 × 10-8)) in individuals with AUD. Postmortem brain analyses documented increased SLC7A11 expression in the frontal cortex of individuals with AUD and animal models showed marked increased expression in liver, suggesting a mechanism by which alcohol leads to hypomethylation-induced overexpression of SLC7A11. Taken together, our EWAS discovery sample and subsequent validation of the top probe in AUD suggest a strong role of abnormal glutamate signaling mediated by methylomic variation in SLC7A11. Our data are intriguing given the prominent role of glutamate signaling in brain and liver and might provide an important target for therapeutic intervention.

Methylome-wide association study of antidepressant use in Generation Scotland and the Netherlands Twin Register implicates the innate immune system

Antidepressants are an effective treatment for major depressive disorder (MDD), although individual response is unpredictable and highly variable. Whilst the mode of action of antidepressants is incompletely understood, many medications are associated with changes in DNA methylation in genes that are plausibly linked to their mechanisms. Studies of DNA methylation may therefore reveal the biological processes underpinning the efficacy and side effects of antidepressants. We performed a methylome-wide association study (MWAS) of self-reported antidepressant use accounting for lifestyle factors and MDD in Generation Scotland (GS:SFHS, N = 6428, EPIC array) and the Netherlands Twin Register (NTR, N = 2449, 450 K array) and ran a meta-analysis of antidepressant use across these two cohorts. We found ten CpG sites significantly associated with self-reported antidepressant use in GS:SFHS, with the top CpG located within a gene previously associated with mental health disorders, ATP6V1B2 (β = -0.055, pcorrected = 0.005). Other top loci were annotated to genes including CASP10, TMBIM1, MAPKAPK3, and HEBP2, which have previously been implicated in the innate immune response. Next, using penalised regression, we trained a methylation-based score of self-reported antidepressant use in a subset of 3799 GS:SFHS individuals that predicted antidepressant use in a second subset of GS:SFHS (N = 3360, β = 0.377, p = 3.12 × 10-11, R2 = 2.12%). In an MWAS analysis of prescribed selective serotonin reuptake inhibitors, we showed convergent findings with those based on self-report. In NTR, we did not find any CpGs significantly associated with antidepressant use. The meta-analysis identified the two CpGs of the ten above that were common to the two arrays used as being significantly associated with antidepressant use, although the effect was in the opposite direction for one of them. Antidepressants were associated with epigenetic alterations in loci previously associated with mental health disorders and the innate immune system. These changes predicted self-reported antidepressant use in a subset of GS:SFHS and identified processes that may be relevant to our mechanistic understanding of clinically relevant antidepressant drug actions and side effects.

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.

Epigenetic prediction of major depressive disorder

Variation in DNA methylation (DNAm) is associated with lifestyle factors such as smoking and body mass index (BMI) but there has been little research exploring its ability to identify individuals with major depressive disorder (MDD). Using penalised regression on genome-wide CpG methylation, we tested whether DNAm risk scores (MRS), trained on 1223 MDD cases and 1824 controls, could discriminate between cases (n = 363) and controls (n = 1417) in an independent sample, comparing their predictive accuracy to polygenic risk scores (PRS). The MRS explained 1.75% of the variance in MDD (β = 0.338, p = 1.17 × 10-7) and remained associated after adjustment for lifestyle factors (β = 0.219, p = 0.001, R2 = 0.68%). When modelled alongside PRS (β = 0.384, p = 4.69 × 10-9) the MRS remained associated with MDD (β = 0.327, p = 5.66 × 10-7). The MRS was also associated with incident cases of MDD who were well at recruitment but went on to develop MDD at a later assessment (β = 0.193, p = 0.016, R2 = 0.52%). Heritability analyses found additive genetic effects explained 22% of variance in the MRS, with a further 19% explained by pedigree-associated genetic effects and 16% by the shared couple environment. Smoking status was also strongly associated with MRS (β = 0.440, p ≤ 2 × 10-16). After removing smokers from the training set, the MRS strongly associated with BMI (β = 0.053, p = 0.021). We tested the association of MRS with 61 behavioural phenotypes and found that whilst PRS were associated with psychosocial and mental health phenotypes, MRS were more strongly associated with lifestyle and sociodemographic factors. DNAm-based risk scores of MDD significantly discriminated MDD cases from controls in an independent dataset and may represent an archive of exposures to lifestyle factors that are relevant to the prediction of MDD.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

Birth weight associations with DNA methylation differences in an adult population

The Developmental Origins of Health and Disease (DOHaD) theory predicts that prenatal and early life events shape adult health outcomes. Birth weight is a useful indicator of the foetal experience and has been associated with multiple adult health outcomes. DNA methylation (DNAm) is one plausible mechanism behind the relationship of birth weight to adult health. Through data linkage between Generation Scotland and historic Scottish birth cohorts, and birth records held through the NHS Information and Statistics Division, a sample of 1,757 individuals with available birth weight and DNAm data was derived. Epigenome-wide association studies (EWAS) were performed in two independently generated DNAm subgroups (nSet1 = 1,395, nSet2 = 362), relating adult DNAm from whole blood to birth weight. Meta-analysis yielded one genome-wide significant CpG site (p = 5.97x10-9), cg00966482. There was minimal evidence for attenuation of the effect sizes for the lead loci upon adjustment for numerous potential confounder variables (body mass index, educational attainment, and socioeconomic status). Associations between birth weight and epigenetic measures of biological age were also assessed. Associations between lower birth weight and higher Grim Age acceleration (p(FDR) = 3.6x10-3) and shorter DNAm-derived telomere length (p(FDR) = 1.7x10-3) are described, although results for three other epigenetic clocks were null. Our results provide support for an association between birth weight and DNAm both locally at one CpG site, and globally via biological ageing estimates.

Structural brain correlates of serum and epigenetic markers of inflammation in major depressive disorder

Inflammatory processes are implicated in the aetiology of Major Depressive Disorder (MDD); however, the relationship between peripheral inflammation, brain structure and depression remains unclear, partly due to complexities around the use of acute/phasic inflammatory biomarkers. Here, we report the first large-scale study of both serological and methylomic signatures of CRP (considered to represent acute and chronic measures of inflammation respectively) and their associations with depression status/symptoms, and structural neuroimaging phenotypes (T1 and diffusion MRI) in a large community-based sample (Generation Scotland; NMDD cases = 271, Ncontrols = 609). Serum CRP was associated with overall MDD severity, and specifically with current somatic symptoms- general interest (β = 0.145, PFDR = 6 × 10-4) and energy levels (β = 0.101, PFDR = 0.027), along with reduced entorhinal cortex thickness (β = -0.095, PFDR = 0.037). DNAm CRP was significantly associated with reduced global grey matter/cortical volume and widespread reductions in integrity of 16/24 white matter tracts (with greatest regional effects in the external and internal capsules, βFA= -0.12 to -0.14). In general, the methylation-based measures showed stronger associations with imaging metrics than serum-based CRP measures (βaverage = -0.15 versus βaverage = 0.01 respectively). These findings provide evidence for central effects of peripheral inflammation from both serological and epigenetic markers of inflammation, including in brain regions previously implicated in depression. This suggests that these imaging measures may be involved in the relationship between peripheral inflammation and somatic/depressive symptoms. Notably, greater effects on brain morphology were seen for methylation-based rather than serum-based measures of inflammation, indicating the importance of such measures for future studies.

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.

Epigenome-wide analyses identify DNA methylation signatures of dementia risk

INTRODUCTION

Dementia pathogenesis begins years before clinical symptom onset, necessitating the understanding of premorbid risk mechanisms. Here we investigated potential pathogenic mechanisms by assessing DNA methylation associations with dementia risk factors in Alzheimer's disease (AD)-free participants.

METHODS

Associations between dementia risk measures (family history, AD genetic risk score [GRS], and dementia risk scores [combining lifestyle, demographic, and genetic factors]) and whole-blood DNA methylation were assessed in discovery and replication samples (n = ~400 to ~5000) from Generation Scotland.

RESULTS

AD genetic risk and two dementia risk scores were associated with differential methylation. The GRS associated predominantly with methylation differences in cis but also identified a genomic region implicated in Parkinson disease. Loci associated with dementia risk scores were enriched for those previously associated with body mass index and alcohol consumption.

DISCUSSION

Dementia risk measures show widespread association with blood-based methylation, generating several hypotheses for assessment by future studies.

Characterisation of an inflammation-related epigenetic score and its association with cognitive ability

BACKGROUND

Chronic systemic inflammation has been associated with incident dementia, but its association with age-related cognitive decline is less clear. The acute responses of many inflammatory biomarkers mean they may provide an unreliable picture of the chronicity of inflammation. Recently, a large-scale epigenome-wide association study identified DNA methylation correlates of C-reactive protein (CRP)-a widely used acute-phase inflammatory biomarker. DNA methylation is thought to be relatively stable in the short term, marking it as a potentially useful signature of exposure.

METHODS

We utilise a DNA methylation-based score for CRP and investigate its trajectories with age, and associations with cognitive ability in comparison with serum CRP and a genetic CRP score in a longitudinal study of older adults (n = 889) and a large, cross-sectional cohort (n = 7028).

RESULTS

We identified no homogeneous trajectories of serum CRP with age across the cohorts, whereas the epigenetic CRP score was consistently found to increase with age (standardised β = 0.07 and 0.01) and to do so more rapidly in males compared to females. Additionally, the epigenetic CRP score had higher test-retest reliability compared to serum CRP, indicating its enhanced temporal stability. Higher serum CRP was not found to be associated with poorer cognitive ability (standardised β = - 0.08 and - 0.05); however, a consistent negative association was identified between cognitive ability and the epigenetic CRP score in both cohorts (standardised β = - 0.15 and - 0.08).

CONCLUSIONS

An epigenetic proxy of CRP may provide a more reliable signature of chronic inflammation, allowing for more accurate stratification of individuals, and thus clearer inference of associations with incident health outcomes.

DNA methylation in APOE: The relationship with Alzheimer's and with cardiovascular health

INTRODUCTION

Genetic variation in the apolipoprotein E (APOE) gene is associated with Alzheimer's disease (AD) and risk factors for cardiovascular disease (CVD). DNA methylationat APOE has been associated with altered cognition and AD. It is unclear if epigenetic marks could be used for predicting future disease.

METHODS

We assessed blood-based DNA methylation at 13 CpGs in the APOE gene in 5828 participants from the Generation Scotland (GS) cohort. Using linear mixed models regression, we examined the relationships among APOE methylation, cognition, cholesterol, the family history of AD and the risk for CVD.

RESULTS

DNA methylation at two CpGs was associated with the ratio of total cholesterol and HDL cholesterol, but not with cognition, family history of AD, or the risk of CVD.

DISCUSSION

APOE methylation is associated with the levels of blood cholesterol, but there is no evidence for the utility of APOE methylation as a biomarker for predicting AD or CVD.

An epigenome-wide association study of sex-specific chronological ageing

BACKGROUND

Advanced age is associated with cognitive and physical decline and is a major risk factor for a multitude of disorders. There is also a gap in life expectancy between males and females. DNA methylation differences have been shown to be associated with both age and sex. Here, we investigate age-by-sex differences in blood-based DNA methylation in an unrelated cohort of 2586 individuals between the ages of 18 and 87 years, with replication in a further 4450 individuals between the ages of 18 and 93 years.

METHODS

Linear regression models were applied, with stringent genome-wide significance thresholds (p < 3.6 × 10-8) used in both the discovery and replication data. A second, highly conservative mixed linear model method that better controls the false-positive rate was also applied, using the same genome-wide significance thresholds.

RESULTS

Using the linear regression method, 52 autosomal and 597 X-linked CpG sites, mapping to 251 unique genes, replicated with concordant effect size directions in the age-by-sex interaction analysis. The site with the greatest difference mapped to GAGE10, an X-linked gene. Here, DNA methylation levels remained stable across the male adult age range (DNA methylation by age r = 0.02) but decreased across female adult age range (DNA methylation by age r = - 0.61). One site (cg23722529) with a significant age-by-sex interaction also had a quantitative trait locus (rs17321482) that is a genome-wide significant variant for prostate cancer. The mixed linear model method identified 11 CpG sites associated with the age-by-sex interaction.

CONCLUSION

The majority of differences in age-associated DNA methylation trajectories between sexes are present on the X chromosome. Several of these differences occur within genes that have been implicated in sexually dimorphic traits.

A meta-analysis of genome-wide association studies of epigenetic age acceleration

'Epigenetic age acceleration' is a valuable biomarker of ageing, predictive of morbidity and mortality, but for which the underlying biological mechanisms are not well established. Two commonly used measures, derived from DNA methylation, are Horvath-based (Horvath-EAA) and Hannum-based (Hannum-EAA) epigenetic age acceleration. We conducted genome-wide association studies of Horvath-EAA and Hannum-EAA in 13,493 unrelated individuals of European ancestry, to elucidate genetic determinants of differential epigenetic ageing. We identified ten independent SNPs associated with Horvath-EAA, five of which are novel. We also report 21 Horvath-EAA-associated genes including several involved in metabolism (NHLRC, TPMT) and immune system pathways (TRIM59, EDARADD). GWAS of Hannum-EAA identified one associated variant (rs1005277), and implicated 12 genes including several involved in innate immune system pathways (UBE2D3, MANBA, TRIM46), with metabolic functions (UBE2D3, MANBA), or linked to lifespan regulation (CISD2). Both measures had nominal inverse genetic correlations with father's age at death, a rough proxy for lifespan. Nominally significant genetic correlations between Hannum-EAA and lifestyle factors including smoking behaviours and education support the hypothesis that Hannum-based epigenetic ageing is sensitive to variations in environment, whereas Horvath-EAA is a more stable cellular ageing process. We identified novel SNPs and genes associated with epigenetic age acceleration, and highlighted differences in the genetic architecture of Horvath-based and Hannum-based epigenetic ageing measures. Understanding the biological mechanisms underlying individual differences in the rate of epigenetic ageing could help explain different trajectories of age-related decline.

A meta-analysis of genome-wide association studies of epigenetic age acceleration

'Epigenetic age acceleration' is a valuable biomarker of ageing, predictive of morbidity and mortality, but for which the underlying biological mechanisms are not well established. Two commonly used measures, derived from DNA methylation, are Horvath-based (Horvath-EAA) and Hannum-based (Hannum-EAA) epigenetic age acceleration. We conducted genome-wide association studies of Horvath-EAA and Hannum-EAA in 13,493 unrelated individuals of European ancestry, to elucidate genetic determinants of differential epigenetic ageing. We identified ten independent SNPs associated with Horvath-EAA, five of which are novel. We also report 21 Horvath-EAA-associated genes including several involved in metabolism (NHLRC, TPMT) and immune system pathways (TRIM59, EDARADD). GWAS of Hannum-EAA identified one associated variant (rs1005277), and implicated 12 genes including several involved in innate immune system pathways (UBE2D3, MANBA, TRIM46), with metabolic functions (UBE2D3, MANBA), or linked to lifespan regulation (CISD2). Both measures had nominal inverse genetic correlations with father's age at death, a rough proxy for lifespan. Nominally significant genetic correlations between Hannum-EAA and lifestyle factors including smoking behaviours and education support the hypothesis that Hannum-based epigenetic ageing is sensitive to variations in environment, whereas Horvath-EAA is a more stable cellular ageing process. We identified novel SNPs and genes associated with epigenetic age acceleration, and highlighted differences in the genetic architecture of Horvath-based and Hannum-based epigenetic ageing measures. Understanding the biological mechanisms underlying individual differences in the rate of epigenetic ageing could help explain different trajectories of age-related decline.

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.

Identification of novel differentially methylated sites with potential as clinical predictors of impaired respiratory function and COPD

Background

The causes of poor respiratory function and COPD are incompletely understood, but it is clear that genes and the environment play a role. As DNA methylation is under both genetic and environmental control, we hypothesised that investigation of differential methylation associated with these phenotypes would permit mechanistic insights, and improve prediction of COPD. We investigated genome-wide differential DNA methylation patterns using the recently released 850 K Illumina EPIC array. This is the largest single population, whole-genome epigenetic study to date.

Methods

Epigenome-wide association studies (EWASs) of respiratory function and COPD were performed in peripheral blood samples from the Generation Scotland: Scottish Family Health Study (GS:SFHS) cohort (n = 3781; 274 COPD cases and 2919 controls). In independent COPD incidence data (n = 149), significantly differentially methylated sites (DMSs; p < 3.6 × 10⁻⁸) were evaluated for their added predictive power when added to a model including clinical variables, age, sex, height and smoking history using receiver operating characteristic analysis. The Lothian Birth Cohort 1936 (LBC1936) was used to replicate association (n = 895) and prediction (n = 178) results.

Findings

We identified 28 respiratory function and/or COPD associated DMSs, which mapped to genes involved in alternative splicing, JAK-STAT signalling, and axon guidance. In prediction analyses, we observed significant improvement in discrimination between COPD cases and controls (p < .05) in independent GS:SFHS (p = .016) and LBC1936 (p = .010) datasets by adding DMSs to a clinical model.

Interpretation

Identification of novel DMSs has provided insight into the molecular mechanisms regulating respiratory function and aided prediction of COPD risk. Further studies are needed to assess the causality and clinical utility of identified associations.

Fund

Wellcome Trust Strategic Award 10436/Z/14/Z.

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.

Assessment of dried blood spots for DNA methylation profiling

Background: DNA methylation reflects health-related environmental exposures and genetic risk, providing insights into aetiological mechanisms and potentially predicting disease onset, progression and treatment response. An increasingly recognised need for large-scale, longitudinally-profiled samples collected world-wide has made the development of efficient and straightforward sample collection and storage procedures a pressing issue. An alternative to the low-temperature storage of EDTA tubes of venous blood samples, which are frequently the source of the DNA used in such studies, is to collect and store at room temperature blood samples using purpose built filter paper, such as Whatman FTA® cards. Our goal was to determine whether DNA stored in this manner can be used to generate DNA methylation profiles comparable to those generated using blood samples frozen in EDTA tubes.

Methods: DNA methylation profiles were obtained from matched EDTA tube and Whatman FTA® card whole-blood samples from 62 Generation Scotland: Scottish Family Health Study participants using the Infinium HumanMethylation450 BeadChip. Multiple quality control procedures were implemented, the relationship between the two sample types assessed, and epigenome-wide association studies (EWASs) performed for smoking status, age and the interaction between these variables and sample storage method.

Results: Dried blood spot (DBS) DNA methylation profiles were of good quality and DNA methylation profiles from matched DBS and EDTA tube samples were highly correlated (mean r = 0.991) and could distinguish between participants. EWASs replicated established associations for smoking and age, with no evidence for moderation by storage method.

Conclusions: Our results support the use of Whatman FTA® cards for collecting and storing blood samples for DNA methylation profiling. This approach is likely to be particularly beneficial for large-scale studies and those carried out in areas where freezer access is limited. Furthermore, our results will inform consideration of the use of newborn heel prick DBSs for research use.

Abstract P3-08-02: Withdrawn

This abstract was withdrawn by the authors.

Citation Format: Harano K, Wang Y, Masuda H, Lim B, Parinyanitikul N, Lee H-J, Seitz RS, Morris SW, Bailey DB, Hout DR, Rao A, Lucci A, Tripathy D, Krishnamurthy S, Ueno NT. Withdrawn [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P3-08-02.

Abstract P4-02-05: Apocrine morphology and LAR molecular subtype predict prognosis of TNBC patients with residual disease after neoadjuvant chemotherapy

Background: TNBC molecular subtype classification updated by Lehmann et al. includes 4 subtypes: basal-like 1 and 2 (BL1), (BL2), mesenchymal (M), and luminal androgen receptor (LAR), and as a modifier of these subtypes, an Immunomodulatory (IM) gene expression signature. However, molecular subtypes have not been linked to morphological features of TNBC. Apocrine carcinoma has been proposed as a TNBC category that expresses androgen receptor. LAR-subtype TNBC has a poor response to neoadjuvant systemic therapy (NST). We hypothesized that defining the apocrine-featured TNBC by morphology and molecular subtype predict the prognosis of patients with residual disease after NST. Methods: We created the Pan-Pacific TNBC Consortium dataset, which contains paired samples of matched pre and post-NST TNBC tumors from 4 institutions. All patients received NST and didn't have a pathological complete response (pCR). Three pathologists examined hematoxylin and eosin-stained slides of 86 pre-NST samples and determined (1) the presence of apocrine differentiation, (2) the level of tumor-infiltrating lymphocytes (TILs), (3) the histological grade (HG), and (4) the rate of necrosis. These morphological features were compared among the subtypes. For a sample to be considered apocrine positive, apocrine differentiation had to be identified by 2 or more pathologists. Fisher's exact test was used to test the association of subtypes and morphological features. The log-rank test was used to compare disease-free survival (DFS). Results: Twelve of 24 (50%) apocrine-positive tumor samples were LAR subtype, and12 of 17 (70%) LAR-subtype tumor samples exhibited apocrine differentiation. The other subtypes showed following: BL1, 11/44 (25%); BL2, 0/7 (0%); M, 1/10 (10%); unclassified, 0/8 (0%). The median follow-up time was 22 months. In all populations, 2-year DFS rates were higher in patients with apocrine-positive tumors than in those whose tumors did not exhibit apocrine differentiation (P = .027; 2-year DFS, 85% vs 54%). The LAR subtype was also associated with lower HG, although LAR tumors had a similar prognosis to the other subtypes. In the combined analysis of subtypes and apocrine differentiation, patients with apocrine-positive LAR tumors had a higher 2-year DFS rate than did those with apocrine-negative LAR tumors (P = .044; 2-year DFS, 88% vs. 30%). However, patients with apocrine-positive BL1 tumors had no better DFS than did those with apocrine-negative BL1 tumors (P = .133). TIL levels and the presence of the IM signature were positively associated (P = .01), and apocrine differentiation positivity tended to be negatively associated with TIL level (P = .06). Neither TIL level nor IM signature was associated with survival. Conclusion: Apocrine differentiation was associated with the LAR subtype of TNBC and better prognosis in patients who did not have a pCR. The LAR subtype alone did not predict DFS; however, LAR tumors with apocrine differentiation had a better prognosis than did LAR tumors without apocrine differentiation. Using a combination of morphologic and genomic testing may be helpful in determining the prognosis of patients with apocrine-positive TNBC tumors who have residual disease after NST.

Citation Format: Masuda H, Miura S, Harano K, Wang Y, Hirota Y, Matsunaga Y, Lim B, Lucci A, Parinyanitikul N, Lee HJ, Gong G, Rao A, Seitz RS, Morris SW, Hout DR, Nakamura S, Tripathy D, Harada O, Krishnamurthy S, Ueno NT. Apocrine morphology and LAR molecular subtype predict prognosis of TNBC patients with residual disease after neoadjuvant chemotherapy [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P4-02-05.

An epigenetic score for BMI based on DNA methylation correlates with poor physical health and major disease in the Lothian Birth Cohort

BACKGROUND

The relationship between obesity and adverse health is well established, but little is known about the contribution of DNA methylation to obesity-related health outcomes. This study tests associations between an epigenetic score for body mass index (BMI) and health-related, cognitive, psychosocial and lifestyle outcomes in the Lothian Birth Cohort 1936. This study also tests whether these associations are independent of phenotypic BMI.

METHOD

Analyses were conducted using data from the Lothian Birth Cohort 1936 (n = 892). Weights for the epigenetic BMI score were derived using penalised regression on methylation data from unrelated Generation Scotland participants (n = 2562). Associations were tested for replication in an independent sample: the Lothian Birth Cohort 1921 (n = 433).

RESULTS

A higher epigenetic BMI score was associated with higher BMI (R² = 0.1), greater body weight (R² = 0.06), greater time taken to walk 6 m, poorer lung function and poorer general physical health (all R² = 0.02), greater levels of triglycerides (R² = 0.09), greater %total HbA1c (R² = 0.06), lower levels of high-density lipoprotein cholesterol (HDL; R² = 0.08), higher HDL ratio (HDL/total cholesterol; R² = 0.03), lower health-related quality of life, physical inactivity, and greater social deprivation (all R² = 0.02). The epigenetic BMI score (per SD) was also associated with type 2 diabetes (OR 2.17, 95% CI 1.67, 2.84), cardiovascular disease (OR 1.45, 95% CI 1.24, 1.71) and high blood pressure (OR 1.30, 95% CI 1.13, 1.49; all p < 0.00026 after Bonferroni correction). Associations were replicated for BMI (R² = 0.06), body weight (R² = 0.04), health-related quality of life (R² = 0.02), HbA1c (R² = 0.07) and triglycerides (R² = 0.07; all p < 0.0045 after Bonferroni correction).

CONCLUSIONS

We observed and replicated associations between an epigenetic score for BMI and variables related to poor physical health and metabolic syndrome. Regression models with both epigenetic and phenotypic BMI scores as predictors accounted for a greater proportion of variance in all outcome variables than either predictor alone, demonstrating independent and additive effects of epigenetic and phenotypic BMI scores.

Epigenetic prediction of complex traits and death

BACKGROUND

Genome-wide DNA methylation (DNAm) profiling has allowed for the development of molecular predictors for a multitude of traits and diseases. Such predictors may be more accurate than the self-reported phenotypes and could have clinical applications.

RESULTS

Here, penalized regression models are used to develop DNAm predictors for ten modifiable health and lifestyle factors in a cohort of 5087 individuals. Using an independent test cohort comprising 895 individuals, the proportion of phenotypic variance explained in each trait is examined for DNAm-based and genetic predictors. Receiver operator characteristic curves are generated to investigate the predictive performance of DNAm-based predictors, using dichotomized phenotypes. The relationship between DNAm scores and all-cause mortality (n = 212 events) is assessed via Cox proportional hazards models. DNAm predictors for smoking, alcohol, education, and waist-to-hip ratio are shown to predict mortality in multivariate models. The predictors show moderate discrimination of obesity, alcohol consumption, and HDL cholesterol. There is excellent discrimination of current smoking status, poorer discrimination of college-educated individuals and those with high total cholesterol, LDL with remnant cholesterol, and total:HDL cholesterol ratios.

CONCLUSIONS

DNAm predictors correlate with lifestyle factors that are associated with health and mortality. They may supplement DNAm-based predictors of age to identify the lifestyle profiles of individuals and predict disease risk.

Investigating the relationship between DNA methylation age acceleration and risk factors for Alzheimer's disease

Introduction

The “epigenetic clock” is a DNA methylation–based estimate of biological age and is correlated with chronological age—the greatest risk factor for Alzheimer's disease (AD). Genetic and environmental risk factors exist for AD, several of which are potentially modifiable. In this study, we assess the relationship between the epigenetic clock and AD risk factors.

Methods

Multilevel models were used to assess the relationship between age acceleration (the residual of biological age regressed onto chronological age) and AD risk factors relating to cognitive reserve, lifestyle, disease, and genetics in the Generation Scotland study (n = 5100).

Results

We report significant associations between age acceleration and body mass index, total cholesterol to high-density lipoprotein cholesterol ratios, socioeconomic status, high blood pressure, and smoking behavior (Bonferroni-adjusted P < .05).

Discussion

Associations are present between environmental risk factors for AD and age acceleration. Measures to modify such risk factors might improve the risk profile for AD and the rate of biological ageing. Future longitudinal analyses are therefore warranted.

Altered DNA methylation associated with a translocation linked to major mental illness

Recent work has highlighted a possible role for altered epigenetic modifications, including differential DNA methylation, in susceptibility to psychiatric illness. Here, we investigate blood-based DNA methylation in a large family where a balanced translocation between chromosomes 1 and 11 shows genome-wide significant linkage to psychiatric illness. Genome-wide DNA methylation was profiled in whole-blood-derived DNA from 41 individuals using the Infinium HumanMethylation450 BeadChip (Illumina Inc., San Diego, CA). We found significant differences in DNA methylation when translocation carriers (n = 17) were compared to related non-carriers (n = 24) at 13 loci. All but one of the 13 significant differentially methylated positions (DMPs) mapped to the regions surrounding the translocation breakpoints. Methylation levels of five DMPs were associated with genotype at SNPs in linkage disequilibrium with the translocation. Two of the five genes harbouring significant DMPs, DISC1 and DUSP10, have been previously shown to be differentially methylated in schizophrenia. Gene Ontology analysis revealed enrichment for terms relating to neuronal function and neurodevelopment among the genes harbouring the most significant DMPs. Differentially methylated region (DMR) analysis highlighted a number of genes from the MHC region, which has been implicated in psychiatric illness previously through genetic studies. We show that inheritance of a translocation linked to major mental illness is associated with differential DNA methylation at loci implicated in neuronal development/function and in psychiatric illness. As genomic rearrangements are over-represented in individuals with psychiatric illness, such analyses may be valuable more widely in the study of these conditions.

Abstract P1-07-14: Rates of immune infiltration in patients with triple-negative breast cancers by molecular subtype and in patients with inflammatory and non-inflammatory breast cancers

Background

In patients with triple-negative breast cancer (TNBC), tumor-infiltrating lymphocytes (TILs) have been reported to be associated with improved survival. Lehmann et al. identified 6 molecular subtypes of TNBC [basal-like (BL) 1, BL2, mesenchymal (M), mesenchymal stem like (MSL), immunomodulatory (IM), and luminal androgen receptor (LAR)], and we previously reported that TNBC subtype is a predictor of pathologic complete response (pCR). Recently, the IM gene expression signature has been shown to be indicative of the presence of TILs and has been incorporated into TNBC subtyping as a modifier of the other groups rather than a separate subtype. However, the association between TNBC subtype and the presence of TILs is not known. We hypothesized that the BL2 and LAR subtypes, which have low pCR rates, have low rates of immune infiltration. Inflammatory breast cancer (IBC) is an aggressive cancer that is frequently triple-negative. The association between IBC and the presence of TILs also is not known. In this study, we analyzed the association between TNBC molecular subtype and the IM signature and determined whether the IM signature differed between patients with IBC and non-IBC.

Methods

We retrospectively analyzed 88 patients with TNBC from the World IBC Consortium dataset for whom IBC status was known (IBC, n=39; non-IBC, n=49) and tumor gene expression data were available. TNBC specimens were classified using the TNBCtype algorithm (Insight Genetics, Inc., TN, USA), which uses a 101-gene signature. For each tumor, the TNBCtype algorithm reports the TNBC molecular subtype (BL1, BL2, M, MSL, or LAR) and the IM status, which is described as positive (IM+) or negative (IM-). Recently, Fisher's exact test was used to analyze differences in subtype distribution between the IM+ and IM- tumors.

Results

The subtype distribution differed significantly between the IM+ and IM- tumors

IM signature in TNBC subtypesSubtypeTotal (n=88)IM+ (n=32)IM- (n=56)BL13015 (50)15 (50)BL2202 (100)M808 (100)MSL3113 (42)18 (58)LAR121 (8)11 (92)Not determined53 (60)2 (40) (p=0.0087). The majority of IM+ cases occurred in the BL1 and MSL subtypes. No IM+ cases were observed in the BL2 or M subtypes, and only 1 was observed in the LAR subtype. IM+ cases occurred at roughly the same frequency in patients with IBC (33%) and non-IBC (37%, p=0.73).

Conclusions

TNBC molecular subtypes differ in their degree of immune infiltration, and most IM+ TNBCs are of the BL1 and MSL subtypes. Our finding that the proportion of IM+ cases was not different between IBC and non-IBC indicates that TILs are recruited to the tumor microenvironment similarly in IBC and non-IBC tumors. Further, Pietenpol et al recently showed that the MSL signature represents normal stromal cells rather than tumor cells by performing laser-capture microdissection of TNBC specimen. Validation studies are needed to corroborate and further expand upon our findings.

Citation Format: Harano K, Wang Y, Lim B, Seitz RS, Morris SW, Bailey DB, Hout DR, Skelton RL, Ring BZ, Masuda H, Rao AUK, Woodward WA, Reuben JM, Ueno NT. Rates of immune infiltration in patients with triple-negative breast cancers by molecular subtype and in patients with inflammatory and non-inflammatory breast cancers [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P1-07-14.

Balanced translocation linked to psychiatric disorder, glutamate, and cortical structure/function

Rare genetic variants of large effect can help elucidate the pathophysiology of brain disorders. Here we expand the clinical and genetic analyses of a family with a (1;11)(q42;q14.3) translocation multiply affected by major psychiatric illness and test the effect of the translocation on the structure and function of prefrontal, and temporal brain regions. The translocation showed significant linkage (LOD score 6.1) with a clinical phenotype that included schizophrenia, schizoaffective disorder, bipolar disorder, and recurrent major depressive disorder. Translocation carriers showed reduced cortical thickness in the left temporal lobe, which correlated with general psychopathology and positive psychotic symptom severity. They showed reduced gyrification in prefrontal cortex, which correlated with general psychopathology severity. Translocation carriers also showed significantly increased activation in the caudate nucleus on increasing verbal working memory load, as well as statistically significant reductions in the right dorsolateral prefrontal cortex glutamate concentrations. These findings confirm that the t(1;11) translocation is associated with a significantly increased risk of major psychiatric disorder and suggest a general vulnerability to psychopathology through altered cortical structure and function, and decreased glutamate levels.

DNA methylation in a Scottish family multiply affected by bipolar disorder and major depressive disorder

BACKGROUND

Bipolar disorder (BD) is a severe, familial psychiatric condition. Progress in understanding the aetiology of BD has been hampered by substantial phenotypic and genetic heterogeneity. We sought to mitigate these confounders by studying a multi-generational family multiply affected by BD and major depressive disorder (MDD), who carry an illness-linked haplotype on chromosome 4p. Within a family, aetiological heterogeneity is likely to be reduced, thus conferring greater power to detect illness-related changes. As accumulating evidence suggests that altered DNA methylation confers risk for BD and MDD, we compared genome-wide methylation between (i) affected carriers of the linked haplotype (ALH) and married-in controls (MIs), (ii) well unaffected haplotype carriers (ULH) and MI, (iii) ALH and ULH and (iv) all haplotype carriers (LH) and MI.

RESULTS

Nominally significant differences in DNA methylation were observed in all comparisons, with differences withstanding correction for multiple testing when the ALH or LH group was compared to the MIs. In both comparisons, we observed increased methylation at a locus in FANCI, which was accompanied by increased FANCI expression in the ALH group. FANCI is part of the Fanconi anaemia complementation (FANC) gene family, which are mutated in Fanconi anaemia and participate in DNA repair. Interestingly, several FANC genes have been implicated in psychiatric disorders. Regional analyses of methylation differences identified loci implicated in psychiatric illness by genome-wide association studies, including CACNB2 and the major histocompatibility complex. Gene ontology analysis revealed enrichment for methylation differences in neurologically relevant genes.

CONCLUSIONS

Our results highlight altered DNA methylation as a potential mechanism by which the linked haplotype might confer risk for mood disorders. Differences in the phenotypic outcome of haplotype carriers might, in part, arise from additional changes in DNA methylation that converge on neurologically important pathways. Further work is required to investigate the underlying mechanisms and functional consequences of the observed differences in methylation.

SuRFing the genomics wave: an R package for prioritising SNPs by functionality

Identifying functional non-coding variants is one of the greatest unmet challenges in genetics. To help address this, we introduce an R package, SuRFR, which integrates functional annotation and prior biological knowledge to prioritise candidate functional variants. SuRFR is publicly available, modular, flexible, fast, and simple to use. We demonstrate that SuRFR performs with high sensitivity and specificity and provide a widely applicable and scalable benchmarking dataset for model training and validation.

Website: http://www.cgem.ed.ac.uk/resources/

708 Common and 2010 rare DISC1 locus variants identified in 1542 subjects: analysis for association with psychiatric disorder and cognitive traits

A balanced t(1;11) translocation that transects the Disrupted in schizophrenia 1 (DISC1) gene shows genome-wide significant linkage for schizophrenia and recurrent major depressive disorder (rMDD) in a single large Scottish family, but genome-wide and exome sequencing-based association studies have not supported a role for DISC1 in psychiatric illness. To explore DISC1 in more detail, we sequenced 528 kb of the DISC1 locus in 653 cases and 889 controls. We report 2718 validated single-nucleotide polymorphisms (SNPs) of which 2010 have a minor allele frequency of <1%. Only 38% of these variants are reported in the 1000 Genomes Project European subset. This suggests that many DISC1 SNPs remain undiscovered and are essentially private. Rare coding variants identified exclusively in patients were found in likely functional protein domains. Significant region-wide association was observed between rs16856199 and rMDD (P=0.026, unadjusted P=6.3 × 10(-5), OR=3.48). This was not replicated in additional recurrent major depression samples (replication P=0.11). Combined analysis of both the original and replication set supported the original association (P=0.0058, OR=1.46). Evidence for segregation of this variant with disease in families was limited to those of rMDD individuals referred from primary care. Burden analysis for coding and non-coding variants gave nominal associations with diagnosis and measures of mood and cognition. Together, these observations are likely to generalise to other candidate genes for major mental illness and may thus provide guidelines for the design of future studies.

Association analysis of Neuregulin 1 candidate regions in schizophrenia and bipolar disorder

Schizophrenia (SCZ) and bipolar disorder (BPD) are severe heritable psychiatric disorders involving a complex genetic aetiology. Neuregulin 1 (NRG1) is a leading candidate gene for SCZ, and has recently been implicated in BPD. We previously reported association of two NRG1 haplotypes with SCZ and BPD in a Scottish case-control sample. One haplotype is located at the 5' end of the gene (region A), and the other is located at the 3' end (region B). Here, association to haplotypes within regions A and B was assessed in patients with SCZ and BPD in a second Scottish case-control sample and in the two Scottish samples combined. Association to region B was also assessed in patients with SCZ and BPD in a German case-control sample, and in all three samples combined. No evidence was found for association in the new samples when analysed individually; however, in the joint analysis of the two Scottish samples, a region B haplotype comprising two SNPs (rs6988339 and rs3757930) was associated with SCZ and the combined case group (SCZ: p=0.0037, OR=1.3, 95% CI: 1.1-1.6; BPD+SCZ: p=0.0080, OR=1.2, 95% CI: 1.1-1.5), with these associations withstanding multiple testing correction at the single-test level (SCZ: p(st)=0.022; BPD+SCZ: p(st)=0.044). This study supports the involvement of NRG1 variants in the less well studied 3' region in conferring susceptibility to SCZ and BPD in the Scottish population.

Association analysis of the chromosome 4p15-p16 candidate region for bipolar disorder and schizophrenia

Several independent linkage studies have identified chromosome 4p15-p16 as a putative region of susceptibility for bipolar disorder (BP), schizophrenia (SCZ) and related phenotypes. Previously, we identified two subregions (B and D) of the 4p15-p16 region that are shared by three of four 4p-linked families examined. Here, we describe a large-scale association analysis of regions B and D (3.8 and 4.5 Mb, respectively). We selected 408 haplotype-tagging single nucleotide polymorphisms (SNPs) on a block-by-block basis from the International HapMap project and tested them in 368 BP, 386 SCZ and 458 control individuals. Nominal significance thresholds were determined using principal component analysis as implemented in the program SNPSpD. In region B, overlapping SNPs and haplotypes met the region-wide threshold (P<or=0.0005) at the global and individual haplotype test level and clustered in two regions. In region D, no individual SNPs were nominally significant, but multiple global and individual haplotypes were associated with BP and/or SCZ (region-wide threshold, P<or=0.0003). These overlapping haplotypes fell into two regions. Within each of these four clusters, at least one globally significant haplotype withstood permutation testing (P(gp)<or=0.05). Five predicted genes were found within these associated regions, while Known/RefSeq genes, including KIAA0746 and PPARGC1A, mapped nearby. There were also nine other clusters within regions B and D with nominally significant haplotypes, but only at the individual haplotype level. KIAA0746, PPARGC1A, GPR125, CCKAR and DKFZp761B107 overlapped with these regions. This study has identified significant associations between BP and SCZ within the chromosome 4p linkage region, resulting in candidate regions worthy of further investigation.

The PDE4B gene confers sex-specific protection against schizophrenia

BACKGROUND

Phosphodiesterase 4B (PDE4B) is a candidate gene for schizophrenia and affective disorders through its disruption by a chromosomal translocation in an individual with schizophrenia, its inhibition by the antidepressant rolipram, and its physical interaction with another key candidate, Disrupted in Schizophrenia (DISC1).

OBJECTIVE

To determine the contribution made by PDE4B to the population risk of schizophrenia and bipolar disorder by carrying out a case-control association study.

METHODS

Twenty-six tagging single nucleotide polymorphisms were selected across the PDE4B gene and genotyped in DNA samples from 386 schizophrenia cases, 368 bipolar disorder cases and 455 controls.

MAIN RESULTS

Single single nucleotide polymorphisms and a resulting haplotype conferred a protective effect against schizophrenia in the female population. The haplotype result remained significant after correction for multiple testing (P=0.012).

CONCLUSION

The observation that a PDE4B haplotype alters the genetic risk of schizophrenia in the Scottish population complements the known participation of this gene in biological processes associated with mental illness. Further studies are needed to replicate this finding and identify underlying sequence variants.

Haplotype analysis and a novel allele-sharing method refines a chromosome 4p locus linked to bipolar affective disorder

BACKGROUND

Bipolar affective disorder (BPAD) and schizophrenia (SCZ) are common conditions. Their causes are unknown, but they include a substantial genetic component. Previously, we described significant linkage of BPAD to a chromosome 4p locus within a large pedigree (F22). Others subsequently have found evidence for linkage of BPAD and SCZ to this region.

METHODS

We constructed high-resolution haplotypes for four linked families, calculated logarithm of the odds (LOD) scores, and developed a novel method to assess the extent of allele sharing within genes between the families.

RESULTS

We describe an increase in the F22 LOD score for this region. Definition and comparison of the linked haplotypes allowed us to prioritize two subregions of 3.8 and 4.4 Mb. Analysis of the extent of allele sharing within these subregions identified 200 kb that shows increased allele sharing between families.

CONCLUSIONS

Linkage of BPAD to chromosome 4p has been strengthened. Haplotype analysis in the additional linked families refined the 20-Mb linkage region. Development of a novel allele-sharing method allowed us to bridge the gap between conventional linkage and association studies. Description of a 200-kb region of increased allele sharing prioritizes this region, which contains two functional candidate genes for BPAD, SLC2A9, and WDR1, for subsequent studies.

Association of Neuregulin 1 with schizophrenia and bipolar disorder in a second cohort from the Scottish population

Neuregulin 1 (NRG1) is a strong candidate for involvement in the aetiology of schizophrenia. A haplotype, initially identified as showing association in the Icelandic and Scottish populations, has shown a consistent effect size in multiple European populations. Additionally, NRG1 has been implicated in susceptibility to bipolar disorder. In this first study to select markers systematically on the basis of linkage disequilibrium across the entire NRG1 gene, we used haplotype-tagging single-nucleotide polymorphisms to identify single markers and haplotypes associated with schizophrenia and bipolar disorder in an independently ascertained Scottish population. Haplotypes in two regions met an experiment-wide significance threshold of P=0.0016 (Nyholt's SpD) and were permuted to correct for multiple testing. Region A overlaps with the Icelandic haplotype and shows nominal association with schizophrenia (P=0.00032), bipolar disorder (P=0.0011), and the combined case group (P=0.0017). This region includes the 5' exon of the NRG1 GGF2 isoform and overlaps the expressed sequence tag (EST) cluster Hs.97362. However, no haplotype in Region A remains significant after permutation analysis (P>0.05). Region B contains a haplotype associated with both schizophrenia (P=0.00014), and the combined case group (P=0.000062), although it does not meet Nyholt's threshold in bipolar disorder alone (P=0.0022). This haplotype remained significant after permutation analysis in both the schizophrenia and combined case groups (P=0.024 and P=0.016, respectively). It spans a approximately 136 kb region that includes the coding sequence of the sensory and motor neuron derived factor (SMDF) isoform and 3' exons of all other known NRG1 isoforms. Our study identifies a new of NRG1 region involved in schizophrenia and bipolar disorder in the Scottish population.

Cytogenetic and genetic evidence supports a role for the kainate-type glutamate receptor gene, GRIK4, in schizophrenia and bipolar disorder

In the search for the biological causes of schizophrenia and bipolar disorder, glutamate neurotransmission has emerged as one of a number of candidate processes and pathways where underlying gene deficits may be present. The analysis of chromosomal rearrangements in individuals diagnosed with neuropsychiatric disorders is an established route to candidate gene identification in both Mendelian and complex disorders. Here we describe a set of genes disrupted by, or proximal to, chromosomal breakpoints (2p12, 2q31.3, 2q21.2, 11q23.3 and 11q24.2) in a patient where chronic schizophrenia coexists with mild learning disability (US: mental retardation). Of these disrupted genes, the most promising candidate is a member of the kainate-type ionotropic glutamate receptor family, GRIK4 (KA1). A subsequent systematic case-control association study on GRIK4 assessed its contribution to psychiatric illness in the karyotypically normal population. This identified two discrete regions of disease risk within the GRIK4 locus: three single single nucleotide polymorphism (SNP) markers with a corresponding underlying haplotype associated with susceptibility to schizophrenia (P=0.0005, odds ratio (OR) of 1.453, 95% CI 1.182-1.787) and two single SNP markers and a haplotype associated with a protective effect against bipolar disorder (P=0.0002, OR of 0.624, 95% CI 0.485-0.802). After permutation analysis to correct for multiple testing, schizophrenia and bipolar disorder haplotypes remained significant (P=0.0430, s.e. 0.0064 and P=0.0190, s.e. 0.0043, respectively). We propose that these convergent cytogenetic and genetic findings provide molecular evidence for common aetiologies for different psychiatric conditions and further support the 'glutamate hypothesis' of psychotic illness.

Association analysis of the chromosome 4p-located G protein-coupled receptor 78 (GPR78) gene in bipolar affective disorder and schizophrenia

The orphan G protein-coupled receptor 78 (GPR78) gene lies within a region of chromosome 4p where we have previously shown linkage to bipolar affective disorder (BPAD) in a large Scottish family. GPR78 was screened for single-nucleotide polymorphisms (SNPs) and a linkage disequilibrium map was constructed. Six tagging SNPs were selected and tested for association on a sample of 377 BPAD, 392 schizophrenia (SCZ) and 470 control individuals. Using standard chi(2) statistics and a backwards logistic regression approach to adjust for the effect of sex, SNP rs1282, located approximately 3 kb upstream of the coding region, was identified as a potentially important variant in SCZ (chi(2) P=0.044; LRT P=0.065). When the analysis was restricted to females, the strength of association increased to an uncorrected allele P-value of 0.015 (odds ratios (OR)=1.688, 95% confidence intervals (CI): 1.104-2.581) and uncorrected genotype P-value of 0.015 (OR=5.991, 95% CI: 1.545-23.232). Under the recessive model, the genotype P-value improved further to 0.005 (OR=5.618, 95% CI: 1.460-21.617) and remained significant after correcting for multiple testing (P=0.017). No single-marker association was detected in the SCZ males, in the BPAD individuals or with any other SNP. Haplotype analysis of the case-control samples revealed several global and individual haplotypes, with P-values <0.05, all but one of which contained SNP rs1282. After correcting for multiple testing, two haplotypes remained significant in both the female BPAD individuals (P=0.038 and 0.032) and in the full sample of affected female individuals (P=0.044 and 0.033). Our results provide preliminary evidence for the involvement of GPR78 in susceptibility to BPAD and SCZ in the Scottish population.

The emerging normal and disease-related roles of anaplastic lymphoma kinase

Anaplastic lymphoma kinase (ALK) is a receptor tyrosine kinase, the normal role of which remains to be completely elucidated. Although work carried out in mammals suggests a function in neural development, results from studies in Drosophila indicate an additional role in visceral muscle differentiation. The aberrant expression of full-length ALK receptor proteins has been reported in neuroblastomas and glioblastomas, while the occurrence of ALK fusion proteins in anaplastic large cell lymphoma (ALCL) has resulted in the identification of the new tumor entity, ALK-positive ALCL. ALK represents one of few examples of a receptor tyrosine kinase implicated in oncogenesis in both haematopoietic and non-haematopoietic tumors, given that ALK fusions also occur in the mesenchymal tumor known as inflammatory myofibroblastic tumor (IMT). The study of ALK fusion proteins, besides demonstrating their importance in tumor development, has also raised the possibility of new therapeutic treatments for patients with ALK-positive malignancies.

Anaplastic lymphoma kinase proteins in growth control and cancer

The normal functions of full-length anaplastic lymphoma kinase (ALK) remain to be completely elucidated. Although considered to be important in neural development, recent studies in Drosophila also highlight a role for ALK in gut muscle differentiation. Indeed, the Drosophila model offers a future arena for the study of ALK, its ligands and signalling cascades. The discovery of activated fusion forms of the ALK tyrosine kinase in anaplastic large cell lymphoma (ALCL) has dramatically improved our understanding of the pathogenesis of these lymphomas and enhanced the pathological diagnosis of this subtype of non-Hodgkin's lymphoma (NHL). Likewise, the realisation that a high percentage of inflammatory myofibroblastic tumours express activated-ALK fusion proteins has clarified the causation of these mesenchymal neoplasms and provided for their easier discrimination from other mesenchymal-derived inflammatory myofibroblastic tumour (IMT) mimics. Recent reports of ALK expression in a range of carcinoma-derived cell lines together with its apparent role as a receptor for PTN and MK, both of which have been implicated in tumourigenesis, raise the possibility that ALK-mediated signalling could play a role in the development and/or progression of a number of common solid tumours. The therapeutic targeting of ALK may prove to have efficacy in the treatment of many of these neoplasms.

Disruption of the RanBP17/Hox11L2 region by recombination with the TCRdelta locus in acute lymphoblastic leukemias with t(5;14)(q34;q11)

The t(5;14)(q33-34;q11) translocation constitutes a recurrent rearrangement in acute lymphoblastic leukemia involving the T cell receptor (TCR) delta locus on chromosome 14. Breakpoint sequences of the derivative chromosome 5 were isolated by application of a ligation-mediated PCR technique using TCR delta-specific primers to amplify genomic DNA from the leukemic cells of a patient with t(5;14). Through exon trap analysis, we identified various putative exons of the chromosome 5 target gene of the translocation; compilation of sequence information of trapped exons and available expressed sequence tags (ESTs) from the GenBank database allowed us to assemble 1.2 kb of the cDNA. Full-length cDNAs were isolated from a human testis cDNA library and sequence analysis predicted a putative Ran binding protein, a novel member of the importin-beta superfamily of nuclear transport receptors, called RanBP17. The t(5;14) breakpoint maps to the 3' coding region of the gene. The breakpoint of a second t(5;14) positive patient was mapped about 8 kb downstream of the most 3' RanBP17 exon and 2 kb upstream of the first exon of the orphan homeobox gene, Hox11L2. In both cases TCR delta enhancer sequences are juxtaposed downstream of the truncated or intact RanBP17 gene, respectively on the derivative chromosome.

SNP genotyping on pooled DNAs: comparison of genotyping technologies and a semi automated method for data storage and analysis

We have compared the accuracy, efficiency and robustness of three methods of genotyping single nucleotide polymorphisms on pooled DNAs. We conclude that (i) the frequencies of the two alleles in pools should be corrected with a factor for unequal allelic amplification, which should be estimated from the mean ratio of a set of heterozygotes (k); (ii) the repeatability of an assay is more important than pinpoint accuracy when estimating allele frequencies, and assays should therefore be optimised to increase the repeatability; and (iii) the size of a pool has a relatively small effect on the accuracy of allele frequency estimation. We therefore recommend that large pools are genotyped and replicated a minimum of four times. In addition, we describe statistical approaches to allow rigorous comparison of DNA pool results. Finally, we describe an extension to our ACeDB database that facilitates management and analysis of the data generated by association studies.

Computational comparison of human genomic sequence assemblies for a region of chromosome 4

Much of the available human genomic sequence data exist in a fragmentary draft state following the completion of the initial high-volume sequencing performed by the International Human Genome Sequencing Consortium (IHGSC) and Celera Genomics (CG). We compared six draft genome assemblies over a region of chromosome 4p (D4S394-D4S403), two consecutive releases by the IHGSC at University of California, Santa Cruz (UCSC), two consecutive releases from the National Centre for Biotechnology Information (NCBI), the public release from CG, and a hybrid assembly we have produced using IHGSC and CG sequence data. This region presents particular problems for genomic sequence assembly algorithms as it contains a large tandem repeat and is sparsely covered by draft sequences. The six assemblies differed both in terms of their relative coverage of sequence data from the region and in their estimated rates of misassembly. The CG assembly method attained the lowest level of misassembly, whereas NCBI and UCSC assemblies had the highest levels of coverage. All assemblies examined included <60% of the publicly available sequence from the region. At least 6% of the sequence data within the CG assembly for the D4S394-D4S403 region was not present in publicly available sequence data. We also show that even in a problematic region, existing software tools can be used with high-quality mapping data to produce genomic sequence contigs with a low rate of rearrangements.

Expression of the oncogenic NPM-ALK chimeric protein in human lymphoid T-cells inhibits drug-induced, but not Fas-induced apoptosis

Anaplastic large cell lymphomas (ALCLs) are frequently associated with the t(2;5)(p23;q35) translocation, leading to the expression of NPM-ALK, a fusion protein linking nucleophosmin and anaplastic lymphoma kinase, a receptor tyrosine kinase. In ALCLs, dimerization of NPM-ALK leads to constitutive autophosphorylation and activation of the kinase, necessary for NPM-ALK oncogenicity. To investigate whether NPM-ALK, like other oncogenic tyrosine kinases, can inhibit drug-induced apoptosis, we permanently transfected NPM-ALK into Jurkat T-cells. As in ALCLs, NPM-ALK was expressed as a constitutively kinase-active 80 kDa protein, and could be detected by immunocytochemistry in nucleoli, nuclei and cytoplasm. Doxorubicin-induced apoptosis (assessed by cell morphology and annexin V-FITC binding) was significantly inhibited in two independent NPM-ALK-expressing clones (5.2+/-1.8 and 7.5+/-0.8% apoptosis), compared to control vector-transduced cells (36+/-6.7%). Similar results were observed with etoposide. In contrast, Fas-induced apoptosis was not inhibited. Cytochrome c release into the cytosol was delayed in doxorubicin-, but not anti-Fas-treated transfectant cells, indicating that apoptosis inhibition occurred upstream of mitochondrial events. Using NPM-ALK mutants, we demonstrated that inhibition of drug-induced apoptosis: (1) requires functional kinase activity, (2) does not involve phospholipase C-gamma, essential for NPM-ALK-mediated mitogenicity and (3) appears to be phosphoinositide 3-kinase independent, despite a strong Akt/PKB activation observed in wild type NPM-ALK-expressing cells. These results suggest that the NPM-ALK antiapoptotic and mitogenic pathways are distinct.

Anaplastic lymphoma kinase (ALK) expression in the inflammatory myofibroblastic tumor: a comparative immunohistochemical study

Inflammatory myofibroblastic tumor (IMT) is an uncommon mesenchymal neoplasm with a variable histologic appearance that may mimic other spindle cell processes, particularly nodular fasciitis, desmoid tumor, and in intra-abdominal locations, gastrointestinal stromal tumor. Recently, gene fusions involving ALK at chromosome 2p23 have been described in IMTs. The resultant ALK protein overexpression in the myofibroblastic component of these tumors is detectable by immunohistochemistry. We examined 73 IMTs, 20 cases of nodular fasciitis, 15 desmoid fibromatoses, and 15 gastrointestinal stromal tumors by immunohistochemistry using ALK-11, a rabbit polyclonal antibody that recognizes the C-terminus of the protein. ALK positivity was detected in 44 of 73 (60%) IMTs. All cases of nodular fasciitis, desmoid fibromatosis, and gastrointestinal stromal tumors were ALK negative (p < 0.001). These findings demonstrate that ALK positivity is common in IMTs, and immunohistochemistry using anti-ALK antibodies can be helpful in the differential diagnosis of these neoplasms. In addition, anti-ALK staining seems to correlate with those IMTs that have the typical tri-patterned histologic appearance and clinical presentation, providing additional support to the premise that IMT is a distinctive clinicopathologic entity within the broad category of inflammatory pseudotumors.