Haematopoietic stem cell-derived immune cells have reduced X chromosome inactivation skewing in systemic lupus erythematosus.

OBJECTIVES: Systemic lupus erythematosus (SLE) shows a marked female bias in prevalence. X chromosome inactivation (XCI) is the mechanism which randomly silences one X chromosome to equalise gene expression between 46, XX females and 46, XY males. Though XCI is expected to result in a random pattern of mosaicism across tissues, some females display a significantly skewed ratio in immune cells, termed XCI-skew. We tested whether XCI was abnormal in females with SLE and hence contributes to sexual dimorphism. METHODS: We assayed XCI in whole blood DNA in 181 female SLE cases, 796 female healthy controls and 10 twin pairs discordant for SLE. Using regression modelling and intra-twin comparisons, we assessed the effect of SLE on XCI and combined clinical, cellular and genetic data via a polygenic score to explore underlying mechanisms. RESULTS: Accommodating the powerful confounder of age, XCI-skew was reduced in females with SLE compared with controls (p=1.3×10-5), with the greatest effect seen in those with more severe disease. Applying an XCI threshold of >80%, we observed XCI-skew in 6.6% of SLE cases compared with 22% of controls. This difference was not explained by differential white cell counts, medication or genetic susceptibility to SLE. Instead, XCI-skew correlated with a biomarker for type I interferon-regulated gene expression. CONCLUSIONS: These results refute current views on XCI-skew in autoimmunity and suggest, in lupus, XCI patterns of immune cells reflect the impact of disease state, specifically interferon signalling, on the haematopoietic stem cells from which they derive.

Exploring the mediation of DNA methylation across the epigenome between childhood adversity and First Episode of Psychosis-findings from the EU-GEI study.

ABTRACT: Studies conducted in psychotic disorders have shown that DNA-methylation (DNAm) is sensitive to the impact of Childhood Adversity (CA). However, whether it mediates the association between CA and psychosis is yet to be explored. Epigenome wide association studies (EWAS) using the Illumina Infinium-Methylation EPIC array in peripheral blood tissue from 366 First-episode of psychosis and 517 healthy controls was performed. Adversity scores were created for abuse, neglect and composite adversity with the Childhood Trauma Questionnaire (CTQ). Regressions examining (I) CTQ scores with psychosis; (II) with DNAm EWAS level and (III) between DNAm and caseness, adjusted for a variety of confounders were conducted. Divide-Aggregate Composite-null Test for the composite null-hypothesis of no mediation effect was conducted. Enrichment analyses were conducted with missMethyl package and the KEGG database. Our results show that CA was associated with psychosis (Composite: OR = 1.68; p = <0.001; abuse: OR = 2.16; p < 0.001; neglect: OR = 2.27; p = <0.001). None of the CpG sites significantly mediated the adversity-psychosis association after Bonferroni correction (p < 8.1 × 10-8). However, 28, 34 and 29 differentially methylated probes associated with 21, 27, 20 genes passed a less stringent discovery threshold (p < 5 × 10-5) for composite, abuse and neglect respectively, with a lack of overlap between abuse and neglect. These included genes previously associated to psychosis in EWAS studies, such as PANK1, SPEG TBKBP1, TSNARE1 or H2R. Downstream gene ontology analyses did not reveal any biological pathways that survived false discovery rate correction. Although at a non-significant level, DNAm changes in genes previously associated with schizophrenia in EWAS studies may mediate the CA-psychosis association. These results and associated involved processes such as mitochondrial or histaminergic disfunction, immunity or neural signalling requires replication in well powered samples. The lack of overlap between mediating genes associated with abuse and neglect suggests differential biological trajectories linking CA subtypes and psychosis.

Can epigenetics shine a light on the biological pathways underlying major mental disorders?

A significant proportion of the global burden of disease can be attributed to mental illness. Despite important advances in identifying risk factors for mental health conditions, the biological processing underlying causal pathways to disease onset remain poorly understood. This represents a limitation to implement effective prevention and the development of novel pharmacological treatments. Epigenetic mechanisms have emerged as mediators of environmental and genetic risk factors which might play a role in disease onset, including childhood adversity (CA) and cannabis use (CU). Particularly, human research exploring DNA methylation has provided new and promising insights into the role of biological pathways implicated in the aetio-pathogenesis of psychiatric conditions, including: monoaminergic (Serotonin and Dopamine), GABAergic, glutamatergic, neurogenesis, inflammatory and immune response and oxidative stress. While these epigenetic changes have been often studied as disease-specific, similarly to the investigation of environmental risk factors, they are often transdiagnostic. Therefore, we aim to review the existing literature on DNA methylation from human studies of psychiatric diseases (i) to identify epigenetic modifications mapping onto biological pathways either transdiagnostically or specifically related to psychiatric diseases such as Eating Disorders, Post-traumatic Stress Disorder, Bipolar and Psychotic Disorder, Depression, Autism Spectrum Disorder and Anxiety Disorder, and (ii) to investigate a convergence between some of these epigenetic modifications and the exposure to known risk factors for psychiatric disorders such as CA and CU, as well as to other epigenetic confounders in psychiatry research.

Genome-wide DNA methylation profiling identifies convergent molecular signatures associated with idiopathic and syndromic autism in post-mortem human brain tissue.

Autism spectrum disorder (ASD) encompasses a collection of complex neuropsychiatric disorders characterized by deficits in social functioning, communication and repetitive behaviour. Building on recent studies supporting a role for developmentally moderated regulatory genomic variation in the molecular aetiology of ASD, we quantified genome-wide patterns of DNA methylation in 223 post-mortem tissues samples isolated from three brain regions [prefrontal cortex, temporal cortex and cerebellum (CB)] dissected from 43 ASD patients and 38 non-psychiatric control donors. We identified widespread differences in DNA methylation associated with idiopathic ASD (iASD), with consistent signals in both cortical regions that were distinct to those observed in the CB. Individuals carrying a duplication on chromosome 15q (dup15q), representing a genetically defined subtype of ASD, were characterized by striking differences in DNA methylationacross a discrete domain spanning an imprinted gene cluster within the duplicated region. In addition to the dramatic cis-effects on DNA methylation observed in dup15q carriers, we identified convergent methylomic signatures associated with both iASD and dup15q, reflecting the findings from previous studies of gene expression and H3K27ac. Cortical co-methylation network analysis identified a number of co-methylated modules significantly associated with ASD that are enriched for genomic regions annotated to genes involved in the immune system, synaptic signalling and neuronal regulation. Our study represents the first systematic analysis of DNA methylation associated with ASD across multiple brain regions, providing novel evidence for convergent molecular signatures associated with both idiopathic and syndromic autism.

Characterizing genetic and environmental influences on variable DNA methylation using monozygotic and dizygotic twins.

Variation in DNA methylation is being increasingly associated with health and disease outcomes. Although DNA methylation is hypothesized to be a mechanism by which both genetic and non-genetic factors can influence the regulation of gene expression, little is known about the extent to which DNA methylation at specific sites is influenced by heritable as well as environmental factors. We quantified DNA methylation in whole blood at age 18 in a birth cohort of 1,464 individuals comprising 426 monozygotic (MZ) and 306 same-sex dizygotic (DZ) twin pairs. Site-specific levels of DNA methylation were more strongly correlated across the genome between MZ than DZ twins. Structural equation models revealed that although the average contribution of additive genetic influences on DNA methylation across the genome was relatively low, it was notably elevated at the highly variable sites characterized by intermediate levels of DNAm that are most relevant for epigenetic epidemiology. Sites at which variable DNA methylation was most influenced by genetic factors were significantly enriched for DNA methylation quantitative trait loci (mQTL) effects, and overlapped with sites where inter-individual variation correlates across tissues. Finally, we show that DNA methylation at sites robustly associated with environmental exposures such as tobacco smoking and obesity is also influenced by additive genetic effects, highlighting the need to control for genetic background in analyses of exposure-associated DNA methylation differences. Estimates of the contribution of genetic and environmental influences to DNA methylation at all sites profiled in this study are available as a resource for the research community (http://www.epigenomicslab.com/online-data-resources).

DNA methylation of FKBP5 and response to exposure-based psychological therapy.

Differential DNA methylation of the hypothalamic-pituitary-adrenal axis related gene FKBP5 has recently been shown to be associated with varying response to environmental influences and may play a role in how well people respond to psychological treatments. Participants (n = 111) received exposure-based cognitive behavioural therapy (CBT) for agoraphobia with or without panic disorder, or specific phobias. Percentage DNA methylation levels were measured for the promoter region and intron 7 of FKBP5. The association between percentage reduction in clinical severity and change in DNA methylation was tested using linear mixed models. The effect of genotype (rs1360780) was tested by the inclusion of an interaction term. The association between change in DNA methylation and FKBP5 expression was examined. Change in percentage DNA methylation at one CpG site of intron 7 was associated with percentage reduction in severity (ß = -4.26, p = 3.90 × 10-4 ), where a decrease in DNA methylation was associated with greater response to therapy. An interaction was detected between rs1360780 and changes in DNA methylation in the promoter region of FKBP5 on treatment outcome (p = .045) but did not survive correction for multiple testing. Changes in DNA methylation were not associated with FKBP5 expression. Decreasing DNA methylation at one CpG site of intron 7 of FKBP5 was strongly associated with decreasing anxiety severity following exposure-based CBT. In addition, there was suggestive evidence that allele-specific methylation at the promoter region may also be associated with treatment response. The results of this study add to the growing literature demonstrating the role of biological processes such as DNA methylation in response to environmental influences.

Methylomic trajectories across human fetal brain development.

Epigenetic processes play a key role in orchestrating transcriptional regulation during development. The importance of DNA methylation in fetal brain development is highlighted by the dynamic expression of de novo DNA methyltransferases during the perinatal period and neurodevelopmental deficits associated with mutations in the methyl-CpG binding protein 2 (MECP2) gene. However, our knowledge about the temporal changes to the epigenome during fetal brain development has, to date, been limited. We quantified genome-wide patterns of DNA methylation at ∼ 400,000 sites in 179 human fetal brain samples (100 male, 79 female) spanning 23 to 184 d post-conception. We identified highly significant changes in DNA methylation across fetal brain development at >7% of sites, with an enrichment of loci becoming hypomethylated with fetal age. Sites associated with developmental changes in DNA methylation during fetal brain development were significantly underrepresented in promoter regulatory regions but significantly overrepresented in regions flanking CpG islands (shores and shelves) and gene bodies. Highly significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small number of regions showing sex-specific DNA methylation trajectories across brain development. Weighted gene comethylation network analysis (WGCNA) revealed discrete modules of comethylated loci associated with fetal age that are significantly enriched for genes involved in neurodevelopmental processes. This is, to our knowledge, the most extensive study of DNA methylation across human fetal brain development to date, confirming the prenatal period as a time of considerable epigenomic plasticity.

Non-replication of the association between 5HTTLPR and response to psychological therapy for child anxiety disorders.

BACKGROUND: We previously reported an association between 5HTTLPR genotype and outcome following cognitive-behavioural therapy (CBT) in child anxiety (Cohort 1). Children homozygous for the low-expression short-allele showed more positive outcomes. Other similar studies have produced mixed results, with most reporting no association between genotype and CBT outcome. AIMS: To replicate the association between 5HTTLPR and CBT outcome in child anxiety from the Genes for Treatment study (GxT Cohort 2, n = 829). METHOD: Logistic and linear mixed effects models were used to examine the relationship between 5HTTLPR and CBT outcomes. Mega-analyses using both cohorts were performed. RESULTS: There was no significant effect of 5HTTLPR on CBT outcomes in Cohort 2. Mega-analyses identified a significant association between 5HTTLPR and remission from all anxiety disorders at follow-up (odds ratio 0.45, P = 0.014), but not primary anxiety disorder outcomes. CONCLUSIONS: The association between 5HTTLPR genotype and CBT outcome did not replicate. Short-allele homozygotes showed more positive treatment outcomes, but with small, non-significant effects. Future studies would benefit from utilising whole genome approaches and large, homogenous samples.

HPA AXIS RELATED GENES AND RESPONSE TO PSYCHOLOGICAL THERAPIES: GENETICS AND EPIGENETICS.

BACKGROUND: Hypothalamic-pituitary-adrenal (HPA) axis functioning has been implicated in the development of stress-related psychiatric diagnoses and response to adverse life experiences. This study aimed to investigate the association between genetic and epigenetics in HPA axis and response to cognitive behavior therapy (CBT). METHODS: Children with anxiety disorders were recruited into the Genes for Treatment project (GxT, N = 1,152). Polymorphisms of FKBP5 and GR were analyzed for association with response to CBT. Percentage DNA methylation at the FKBP5 and GR promoter regions was measured before and after CBT in a subset (n = 98). Linear mixed effect models were used to investigate the relationship between genotype, DNA methylation, and change in primary anxiety disorder severity (treatment response). RESULTS: Treatment response was not associated with FKBP5 and GR polymorphisms, or pretreatment percentage DNA methylation. However, change in FKBP5 DNA methylation was nominally significantly associated with treatment response. Participants who demonstrated the greatest reduction in severity decreased in percentage DNA methylation during treatment, whereas those with little/no reduction in severity increased in percentage DNA methylation. This effect was driven by those with one or more FKBP5 risk alleles, with no association seen in those with no FKBP5 risk alleles. No significant association was found between GR methylation and response. CONCLUSIONS: Allele-specific change in FKBP5 methylation was associated with treatment response. This is the largest study to date investigating the role of HPA axis related genes in response to a psychological therapy. Furthermore, this is the first study to demonstrate that DNA methylation changes may be associated with response to psychological therapies in a genotype-dependent manner.

Measuring adolescents' exposure to victimization: The Environmental Risk (E-Risk) Longitudinal Twin Study.

This paper presents multilevel findings on adolescents' victimization exposure from a large longitudinal cohort of twins. Data were obtained from the Environmental Risk (E-Risk) Longitudinal Twin Study, an epidemiological study of 2,232 children (1,116 twin pairs) followed to 18 years of age (with 93% retention). To assess adolescent victimization, we combined best practices in survey research on victimization with optimal approaches to measuring life stress and traumatic experiences, and introduce a reliable system for coding severity of victimization. One in three children experienced at least one type of severe victimization during adolescence (crime victimization, peer/sibling victimization, Internet/mobile phone victimization, sexual victimization, family violence, maltreatment, or neglect), and most types of victimization were more prevalent among children from low socioeconomic backgrounds. Exposure to multiple victimization types was common, as was revictimization; over half of those physically maltreated in childhood were also exposed to severe physical violence in adolescence. Biometric twin analyses revealed that environmental factors had the greatest influence on most types of victimization, while severe physical maltreatment from caregivers during adolescence was predominantly influenced by heritable factors. The findings from this study showcase how distinct levels of victimization measurement can be harmonized in large-scale studies of health and development.

Epigenome-Wide DNA Methylation Analysis of Monozygotic Twins Discordant for Diurnal Preference.

Diurnal preference is an individual's preference for daily activities and sleep timing and is strongly correlated with the underlying circadian clock and the sleep-wake cycle validating its use as an indirect circadian measure in humans. Recent research has implicated DNA methylation as a mechanism involved in the regulation of the circadian clock system in humans and other mammals. In order to evaluate the extent of epigenetic differences associated with diurnal preference, we examined genome-wide patterns of DNA methylation in DNA from monozygotic (MZ) twin-pairs discordant for diurnal preference. MZ twins were selected from a longitudinal twin study designed to investigate the interplay of genetic and environmental factors in the development of emotional and behavioral difficulties. Fifteen pairs of MZ twins were identified in which one member scored considerably higher on the Horne-Ostberg Morningness-Eveningness Questionnaire (MEQ) than the other. Genome-wide DNA methylation patterns were assessed in twins' buccal cell DNA using the Illumina Infinium HumanMethylation450 BeadChips. Quality control and data pre-processing was undertaken using the wateRmelon package. Differentially methylated probes (DMPs) were identified using an analysis strategy taking into account both the significance and the magnitude of DNA methylation differences. Our data indicate that DNA methylation differences are detectable in MZ twins discordant for diurnal preference. Moreover, downstream gene ontology (GO) enrichment analysis on the top-ranked diurnal preference associated DMPs revealed significant enrichment of pathways that have been previously associated with circadian rhythm regulation, including cell adhesion processes and calcium ion binding.

Genetic and epigenetic associations of MAOA and NR3C1 with depression and childhood adversities.

Monoamine oxidase A (MAOA) harbours a polymorphic upstream variable-number tandem repeat (u-VNTR). The MAOA-L allele of the u-VNTR leads to decreased gene expression levels in vitro and has been found to increase the risk of conduct disorder in males with childhood adversities. Early-life adversities have been associated with hypermethylation of the glucocorticoid receptor (NR3C1). In this study, we first performed a genetic association analysis of the MAOA u-VNTR using individuals with depression (n = 392) and controls (n = 1276). Next, DNA methylation analyses of MAOA and NR3C1 were performed using saliva samples of depressed and control subgroups. Adult MAOA-L females with childhood adversities were found to have a higher risk of developing depression (p = 0.006) and overall MAOA methylation levels were decreased in depressed females compared to controls (mean depressed, 42% vs. mean controls, 44%; p = 0.04). One specific childhood adversity [early parental death (EPD)] was associated with hypermethylation of NR3C1 close to an NGFI-A binding site (mean EPD, 19% vs. mean non-EPD, 14%; p = 0.005). Regression analysis indicated that this association may be mediated by the MAOA-L allele (adjusted R² = 0.24, ANOVA: F = 23.48, p < 0.001). Conclusively: (1) depression in females may result from a gene × childhood-adversity interaction and/or a dysregulated epigenetic programming of MAOA; (2) childhood-adversity subtypes may differentially impact DNA methylation at NR3C1; (3) baseline MAOA-genotypic variations may affect the extent of NR3C1 methylation.

Anxiety sensitivity in adolescence and young adulthood: the role of stressful life events, 5HTTLPR and their interaction.

BACKGROUND: Cognitive biases have long been hypothesized to influence the development and maintenance of symptoms of internalizing problems. Anxiety sensitivity represents one such bias and refers to sensitivity to the physical and emotional symptoms of anxiety and the belief that these are harmful. Twin studies indicate a role for both environmental and genetic influences on anxiety sensitivity. However, little work has been done specifying environments or genes involved in this phenotype. In light of this, we looked at the association between stressful life events, the serotonin transporter gene polymorphism (5HTTLPR), and anxiety sensitivity in a longitudinal sample of adolescents. METHODS: Stressful life events and anxiety sensitivity were measured in over 1,500 individuals at three time points (mean ages 15, 17, and 20 years). 5HTTLPR was genotyped in 1,109 participants. RESULTS: There was consistent evidence for an association between stressful life events and both anxiety sensitivity and change in anxiety sensitivity over time. Although the effect of independent stressful life events was relatively short lived, dependent stressful life events were associated with anxiety sensitivity over time. There was no evidence for a main effect of 5HTTLPR on anxiety sensitivity. 5HTTLPR genotype did not moderate the effect of stressful life events on anxiety sensitivity. CONCLUSIONS: The current study extends previous work by showing that stressful life events, independent of the individual, explained change in cognitions associated with anxiety and depression. This effect does not, however, appear to be moderated by genotype.

Age acquired skewed X chromosome inactivation is associated with adverse health outcomes in humans.

Background: Ageing is a heterogenous process characterised by cellular and molecular hallmarks, including changes to haematopoietic stem cells and is a primary risk factor for chronic diseases. X chromosome inactivation (XCI) randomly transcriptionally silences either the maternal or paternal X in each cell of 46, XX females to balance the gene expression with 46, XY males. Age acquired XCI-skew describes the preferential selection of cells across a tissue resulting in an imbalance of XCI, which is particularly prevalent in blood tissues of ageing females, and yet its clinical consequences are unknown. Methods: We assayed XCI in 1575 females from the TwinsUK population cohort using DNA extracted from whole blood. We employed prospective, cross-sectional, and intra-twin study designs to characterise the relationship of XCI-skew with molecular and cellular measures of ageing, cardiovascular disease risk, and cancer diagnosis. Results: We demonstrate that XCI-skew is independent of traditional markers of biological ageing and is associated with a haematopoietic bias towards the myeloid lineage. Using an atherosclerotic cardiovascular disease risk score, which captures traditional risk factors, XCI-skew is associated with an increased cardiovascular disease risk both cross-sectionally and within XCI-skew discordant twin pairs. In a prospective 10 year follow-up study, XCI-skew is predictive of future cancer incidence. Conclusions: Our study demonstrates that age acquired XCI-skew captures changes to the haematopoietic stem cell population and has clinical potential as a unique biomarker of chronic disease risk. Funding: KSS acknowledges funding from the Medical Research Council [MR/M004422/1 and MR/R023131/1]. JTB acknowledges funding from the ESRC [ES/N000404/1]. MM acknowledges funding from the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. TwinsUK is funded by the Wellcome Trust, Medical Research Council, European Union, Chronic Disease Research Foundation (CDRF), Zoe Global Ltd and the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London.

Sleep quality and diurnal preference in a sample of young adults: associations with 5HTTLPR, PER3, and CLOCK 3111.

Research investigating associations between specific genes and individual differences with regards to the quality and timing of sleep has primarily focussed on serotonin-related and clock genes. However, there are only a few studies of this type and most of those to date have not considered the possibility of gene-environment interaction. Here, we describe associations between sleep quality and diurnal preference and three functional polymorphisms: 5HTTLPR, PERIOD3, and CLOCK 3111. Furthermore, we assessed whether associations between genotypes and sleep phenotypes were moderated by negative life events-a test of gene-environment interaction. DNA from buccal swabs was collected from 947 individuals [mean age = 20.3 years (SD = 1.77), age range = 18-27 years; 61.8% female] and genotyped for the three polymorphisms. Participants completed the Pittsburgh Sleep Quality Index and the Morningness-Eveningness Questionnaire. There was a significant main effect of 5HTTLPR on sleep quality, indicating that "long-long" homozygotes experienced significantly poorer sleep quality (mean = 6.35, SD = 3.36) than carriers of at least one "short" allele (mean = 5.67, SD = 2.96; ß = -0.34, P = 0.005). There were no main effects of 5HTTLPR on diurnal preference; no main effects of PERIOD3 or CLOCK on sleep quality or diurnal preference; and no significant interactions with negative life events. The main effect of the "long" 5HTTLPR allele contradicts previous research, suggesting that perhaps the effects of this gene are heterogeneous in different populations. Failure to replicate previous research in relation to PERIOD3 and CLOCK concurs with previous research suggesting that the effects of these genes are small and may be related to population composition.

DNA methylation signatures of adolescent victimization: analysis of a longitudinal monozygotic twin sample.

Accumulating evidence suggests that individuals exposed to victimization at key developmental stages may have different epigenetic fingerprints compared to those exposed to no/minimal stressful events, however results are inconclusive. This study aimed to strengthen causal inference regarding the impact of adolescent victimization on the epigenome by controlling for genetic variation, age, gender, and shared environmental exposures. We conducted longitudinal epigenome-wide association analyses (EWAS) on DNA methylation (DNAm) profiles of 118 monozygotic (MZ) twin pairs from the Environmental Risk study with and without severe adolescent victimization generated using buccal DNA collected at ages 5, 10 and 18, and the Illumina EPIC array. Additionally, we performed cross-sectional EWAS on age-18 blood and buccal DNA from the same individuals to elucidate tissue-specific signatures of severe adolescent victimization. Our analyses identified 20 suggestive differentially methylated positions (DMPs) (P < 5e-05), with altered DNAm trajectories between ages 10-18 associated with severe adolescent victimization (∆Beta range = -5.5%-5.3%). Age-18 cross-sectional analyses revealed 72 blood (∆Beta range = -2.2%-3.4%) and 42 buccal (∆Beta range = -3.6%-4.6%) suggestive severe adolescent victimization-associated DMPs, with some evidence of convergent signals between these two tissue types. Downstream regional analysis identified significant differentially methylated regions (DMRs) in LGR6 and ANK3 (Sidák P = 5e-09 and 4.07e-06), and one upstream of CCL27 (Sidák P = 2.80e-06) in age-18 blood and buccal EWAS, respectively. Our study represents the first longitudinal MZ twin analysis of DNAm and severe adolescent victimization, providing initial evidence for altered DNA methylomic signatures in individuals exposed to adolescent victimization.

Integrative genomics identifies a convergent molecular subtype that links epigenomic with transcriptomic differences in autism.

Autism spectrum disorder (ASD) is a phenotypically and genetically heterogeneous neurodevelopmental disorder. Despite this heterogeneity, previous studies have shown patterns of molecular convergence in post-mortem brain tissue from autistic subjects. Here, we integrate genome-wide measures of mRNA expression, miRNA expression, DNA methylation, and histone acetylation from ASD and control brains to identify a convergent molecular subtype of ASD with shared dysregulation across both the epigenome and transcriptome. Focusing on this convergent subtype, we substantially expand the repertoire of differentially expressed genes in ASD and identify a component of upregulated immune processes that are associated with hypomethylation. We utilize eQTL and chromosome conformation datasets to link differentially acetylated regions with their cognate genes and identify an enrichment of ASD genetic risk variants in hyperacetylated noncoding regulatory regions linked to neuronal genes. These findings help elucidate how diverse genetic risk factors converge onto specific molecular processes in ASD.

Leveraging epigenetics to examine differences in developmental trajectories of social attention: A proof-of-principle study of DNA methylation in infants with older siblings with autism.

Preliminary evidence suggests that changes in DNA methylation, a widely studied epigenetic mechanism, contribute to the etiology of Autism Spectrum Disorder (ASD). However, data is primarily derived from post-mortem brain samples or peripheral tissue from adults. Deep-phenotyped longitudinal infant cohorts are essential to understand how epigenetic modifications relate to early developmental trajectories and emergence of ASD symptoms. We present a proof-of-principle study designed to evaluate the potential of prospective epigenetic studies of infant siblings of children with ASD. Illumina genome-wide 450 K DNA methylation data from buccal swabs was generated for 63 male infants at multiple time-points from 8 months to 2 years of age (total N = 107 samples). 11 of those infants received a diagnosis of ASD at 3 years. We conducted a series of analyses to characterize DNA methylation signatures associated with categorical outcome and neurocognitive measures from parent-report questionnaire, eye-tracking and electro-encephalography. Effects observed across the entire genome (epigenome-wide association analyses) suggest that collecting DNA methylation samples within infant-sibling designs allows for the detection of meaningful signals with smaller sample sizes than previously estimated. Mapping networks of co-methylated probes associated with neural correlates of social attention implicated enrichment of pathways involved in brain development. Longitudinal modelling found covariation between phenotypic traits and DNA methylation levels in the proximity of genes previously associated with cognitive development, although larger samples and more complete datasets are needed to obtain generalizable results. In conclusion, assessment of DNA methylation profiles at multiple time-points in infant-sibling designs is a promising avenue to comprehend developmental origins and mechanisms of ASD.

Patterns of Reliability: Assessing the Reproducibility and Integrity of DNA Methylation Measurement.

DNA methylation plays an important role in both normal human development and risk of disease. The most utilized method of assessing DNA methylation uses BeadChips, generating an epigenome-wide "snapshot" of >450,000 observations (probe measurements) per assay. However, the reliability of each of these measurements is not equal, and little consideration is paid to consequences for research. We correlated repeat measurements of the same DNA samples using the Illumina HumanMethylation450K and the Infinium MethylationEPIC BeadChips in 350 blood DNA samples. Probes that were reliably measured were more heritable and showed consistent associations with environmental exposures, gene expression, and greater cross-tissue concordance. Unreliable probes were less replicable and generated an unknown volume of false negatives. This serves as a lesson for working with DNA methylation data, but the lessons are equally applicable to working with other data: as we advance toward generating increasingly greater volumes of data, failure to document reliability risks harming reproducibility.

RNA sequencing of identical twins discordant for autism reveals blood-based signatures implicating immune and transcriptional dysregulation.

Background: A gap exists in our mechanistic understanding of how genetic and environmental risk factors converge at the molecular level to result in the emergence of autism symptoms. We compared blood-based gene expression signatures in identical twins concordant and discordant for autism spectrum condition (ASC) to differentiate genetic and environmentally driven transcription differences, and establish convergent evidence for biological mechanisms involved in ASC. Methods: Genome-wide gene expression data were generated using RNA-seq on whole blood samples taken from 16 pairs of monozygotic (MZ) twins and seven twin pair members (39 individuals in total), who had been assessed for ASC and autism traits at age 12. Differential expression (DE) analyses were performed between (a) affected and unaffected subjects (N = 36) and (b) within discordant ASC MZ twin pairs (total N = 11) to identify environmental-driven DE. Gene set enrichment and pathway testing was performed on DE gene lists. Finally, an integrative analysis using DNA methylation data aimed to identify genes with consistent evidence for altered regulation in cis. Results: In the discordant twin analysis, three genes showed evidence for DE at FDR < 10%: IGHG4, EVI2A and SNORD15B. In the case-control analysis, four DE genes were identified at FDR < 10% including IGHG4, PRR13P5, DEPDC1B, and ZNF501. We find enrichment for DE of genes curated in the SFARI human gene database. Pathways showing evidence of enrichment included those related to immune cell signalling and immune response, transcriptional control and cell cycle/proliferation. Integrative methylomic and transcriptomic analysis identified a number of genes showing suggestive evidence for cis dysregulation. Limitations: Identical twins stably discordant for ASC are rare, and as such the sample size was limited and constrained to the use of peripheral blood tissue for transcriptomic and methylomic profiling. Given these primary limitations, we focused on transcript-level analysis. Conclusions: Using a cohort of ASC discordant and concordant MZ twins, we add to the growing body of transcriptomic-based evidence for an immune-based component in the molecular aetiology of ASC. Whilst the sample size was limited, the study demonstrates the utility of the discordant MZ twin design combined with multi-omics integration for maximising the potential to identify disease-associated molecular signals.