Genome-wide association study meta-analysis of suicide death and suicidal behavior.

Suicide is a worldwide health crisis. We aimed to identify genetic risk variants associated with suicide death and suicidal behavior. Meta-analysis for suicide death was performed using 3765 cases from Utah and matching 6572 controls of European ancestry. Meta-analysis for suicidal behavior using data across five cohorts (n = 8315 cases and 256,478 psychiatric or populational controls of European ancestry) was also performed. One locus in neuroligin 1 (NLGN1) passing the genome-wide significance threshold for suicide death was identified (top SNP rs73182688, with p = 5.48 × 10-8 before and p = 4.55 × 10-8 after mtCOJO analysis conditioning on MDD to remove genetic effects on suicide mediated by MDD). Conditioning on suicidal attempts did not significantly change the association strength (p = 6.02 × 10-8), suggesting suicide death specificity. NLGN1 encodes a member of a family of neuronal cell surface proteins. Members of this family act as splice site-specific ligands for beta-neurexins and may be involved in synaptogenesis. The NRXN-NLGN pathway was previously implicated in suicide, autism, and schizophrenia. We additionally identified ROBO2 and ZNF28 associations with suicidal behavior in the meta-analysis across five cohorts in gene-based association analysis using MAGMA. Lastly, we replicated two loci including variants near SOX5 and LOC101928519 associated with suicidal attempts identified in the ISGC and MVP meta-analysis using the independent FinnGen samples. Suicide death and suicidal behavior showed positive genetic correlations with depression, schizophrenia, pain, and suicidal attempt, and negative genetic correlation with educational attainment. These correlations remained significant after conditioning on depression, suggesting pleiotropic effects among these traits. Bidirectional generalized summary-data-based Mendelian randomization analysis suggests that genetic risk for the suicidal attempt and suicide death are both bi-directionally causal for MDD.

Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes.

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

Associations between polygenic risk scores for cardiometabolic phenotypes and adolescent depression and body dissatisfaction.

BACKGROUND: Adolescents with elevated body mass index (BMI) are at an increased risk for depression and body dissatisfaction. Type 2 diabetes (T2D) is an established risk factor for depression. However, shared genetic risk between cardiometabolic conditions and mental health outcomes remains understudied in youth. METHODS: The current study examined associations between polygenic risk scores (PRS) for BMI and T2D, and symptoms of depression and body dissatisfaction, in a sample of 827 community adolescents (Mage = 13.63, SDage = 1.01; 76% girls). BMI, depressive symptoms, and body dissatisfaction were assessed using validated self-report questionnaires. RESULTS: BMI-PRS was associated with phenotypic BMI (ß = 0.24, p < 0.001) and body dissatisfaction (ß = 0.17, p < 0.001), but not with depressive symptoms. The association between BMI-PRS and body dissatisfaction was significantly mediated by BMI (indirect effect = 0.10, CI [0.07-0.13]). T2D-PRS was not associated with depression or body dissatisfaction. CONCLUSIONS: The results suggest phenotypic BMI may largely explain the association between genetic risk for elevated BMI and body dissatisfaction in adolescents. Further research on age-specific genetic effects is needed, as summary statistics from adult discovery samples may have limited utility in youth. IMPACT: The association between genetic risk for elevated BMI and body dissatisfaction in adolescents may be largely explained by phenotypic BMI, indicating a potential pathway through which genetic predisposition influences body image perception. Furthermore, age-specific genetic research is needed to understand the unique influences on health outcomes during adolescence. By identifying BMI as a potential mediator in the association between genetic risk for elevated BMI and body dissatisfaction, the current findings offer insights that could inform interventions targeting body image concerns and mental health in this population.

General v. specific vulnerabilities: polygenic risk scores and higher-order psychopathology dimensions in the Adolescent Brain Cognitive Development (ABCD) Study.

BACKGROUND: Polygenic risk scores (PRSs) capture genetic vulnerability to psychiatric conditions. However, PRSs are often associated with multiple mental health problems in children, complicating their use in research and clinical practice. The current study is the first to systematically test which PRSs associate broadly with all forms of childhood psychopathology, and which PRSs are more specific to one or a handful of forms of psychopathology. METHODS: The sample consisted of 4717 unrelated children (mean age = 9.92, s.d. = 0.62; 47.1% female; all European ancestry). Psychopathology was conceptualized hierarchically as empirically derived general factor (p-factor) and five specific factors: externalizing, internalizing, neurodevelopmental, somatoform, and detachment. Partial correlations explored associations between psychopathology factors and 22 psychopathology-related PRSs. Regressions tested which level of the psychopathology hierarchy was most strongly associated with each PRS. RESULTS: Thirteen PRSs were significantly associated with the general factor, most prominently Chronic Multisite Pain-PRS (r = 0.098), ADHD-PRS (r = 0.079), and Depression-PRS (r = 0.078). After adjusting for the general factor, Depression-PRS, Neuroticism-PRS, PTSD-PRS, Insomnia-PRS, Chronic Back Pain-PRS, and Autism-PRS were not associated with lower order factors. Conversely, several externalizing PRSs, including Adventurousness-PRS and Disinhibition-PRS, remained associated with the externalizing factor (|r| = 0.040-0.058). The ADHD-PRS remained uniquely associated with the neurodevelopmental factor (r = 062). CONCLUSIONS: PRSs developed to predict vulnerability to emotional difficulties and chronic pain generally captured genetic risk for all forms of childhood psychopathology. PRSs developed to predict vulnerability to externalizing difficulties, e.g. disinhibition, tended to be more specific in predicting behavioral problems. The results may inform translation of existing PRSs to pediatric research and future clinical practice.

The benefit of diagnostic whole genome sequencing in schizophrenia and other psychotic disorders.

Schizophrenia has a multifactorial etiology, involving a polygenic architecture. The potential benefit of whole genome sequencing (WGS) in schizophrenia and other psychotic disorders is not well studied. We investigated the yield of clinical WGS analysis in 251 families with a proband diagnosed with schizophrenia (N = 190), schizoaffective disorder (N = 49), or other conditions involving psychosis (N = 48). Participants were recruited in Israel and USA, mainly of Jewish, Arab, and other European ancestries. Trio (parents and proband) WGS was performed for 228 families (90.8%); in the other families, WGS included parents and at least two affected siblings. In the secondary analyses, we evaluated the contribution of rare variant enrichment in particular gene sets, and calculated polygenic risk score (PRS) for schizophrenia. For the primary outcome, diagnostic rate was 6.4%; we found clinically significant, single nucleotide variants (SNVs) or small insertions or deletions (indels) in 14 probands (5.6%), and copy number variants (CNVs) in 2 (0.8%). Significant enrichment of rare loss-of-function variants was observed in a gene set of top schizophrenia candidate genes in affected individuals, compared with population controls (N = 6,840). The PRS for schizophrenia was significantly increased in the affected individuals group, compared to their unaffected relatives. Last, we were also able to provide pharmacogenomics information based on CYP2D6 genotype data for most participants, and determine their antipsychotic metabolizer status. In conclusion, our findings suggest that WGS may have a role in the setting of both research and genetic counseling for individuals with schizophrenia and other psychotic disorders and their families.

Harnessing changes in open chromatin determined by ATAC-seq to generate insulin-responsive reporter constructs.

BACKGROUND: Gene regulation is critical for proper cellular function. Next-generation sequencing technology has revealed the presence of regulatory networks that regulate gene expression and essential cellular functions. Studies investigating the epigenome have begun to uncover the complex mechanisms regulating transcription. Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is quickly becoming the assay of choice for many epigenomic investigations. However, whether intervention-mediated changes in accessible chromatin determined by ATAC-seq can be harnessed to generate intervention-inducible reporter constructs has not been systematically assayed. RESULTS: We used the insulin signaling pathway as a model to investigate chromatin regions and gene expression changes using ATAC- and RNA-seq in insulin-treated Drosophila S2 cells. We found correlations between ATAC- and RNA-seq data, especially when stratifying differentially-accessible chromatin regions by annotated feature type. In particular, our data demonstrated a weak but significant correlation between chromatin regions annotated to enhancers (1-2 kb from the transcription start site) and downstream gene expression. We cloned candidate enhancer regions upstream of luciferase and demonstrate insulin-inducibility of several of these reporters. CONCLUSIONS: Insulin-induced chromatin accessibility determined by ATAC-seq reveals enhancer regions that drive insulin-inducible reporter gene expression.

Polygenic association of glomerular filtration rate decline in world trade center responders.

BACKGROUND: The factors associated with estimated glomerular filtrate rate (eGFR) decline in low risk adults remain relatively unknown. We hypothesized that a polygenic risk score (PRS) will be associated with eGFR decline. METHODS: We analyzed genetic data from 1,601 adult participants with European ancestry in the World Trade Center Health Program (baseline age 49.68 ± 8.79 years, 93% male, 23% hypertensive, 7% diabetic and 1% with cardiovascular disease) with ≥ three serial measures of serum creatinine. PRSs were calculated from an aggregation of single nucleotide polymorphisms (SNPs) from a recent, large-scale genome-wide association study (GWAS) of rapid eGFR decline. Generalized linear models were used to evaluate the association of PRS with renal outcomes: baseline eGFR and CKD stage, rate of change in eGFR, stable versus declining eGFR over a 3-5-year observation period. eGFR decline was defined in separate analyses as "clinical" (> -1.0 ml/min/1.73 m2/year) or "empirical" (lower most quartile of eGFR slopes). RESULTS: The mean baseline eGFR was ~ 86 ml/min/1.73 m2. Subjects with decline in eGFR were more likely to be diabetic. PRS was significantly associated with lower baseline eGFR (B = -0.96, p = 0.002), higher CKD stage (OR = 1.17, p = 0.010), decline in eGFR (OR = 1.14, p = 0.036) relative to stable eGFR, and the lower quartile of eGFR slopes (OR = 1.21, p = 0.008), after adjusting for established risk factors for CKD. CONCLUSION: Common genetic variants are associated with eGFR decline in middle-aged adults with relatively low comorbidity burdens.

Extended familial risk of suicide death is associated with younger age at death and elevated polygenic risk of suicide.

Suicide accounts for >800,000 deaths annually worldwide; prevention is an urgent public health issue. Identification of risk factors remains challenging due to complexity and heterogeneity. The study of suicide deaths with increased extended familial risk provides an avenue to reduce etiological heterogeneity and explore traits associated with increased genetic liability. Using extensive genealogical records, we identified high-risk families where distant relatedness of suicides implicates genetic risk. We compared phenotypic and polygenic risk score (PRS) data between suicides in high-risk extended families (high familial risk (HFR), n = 1,634), suicides linked to genealogical data not in any high-risk families (low familial risk (LFR), n = 147), and suicides not linked to genealogical data with unknown familial risk (UFR, n = 1,865). HFR suicides were associated with lower age at death (mean = 39.34 years), more suicide attempts, and more PTSD and trauma diagnoses. For PRS tests, we included only suicides with >90% European ancestry and adjusted for residual ancestry effects. HFR suicides showed markedly higher PRS of suicide death (calculated using cross-validation), supporting specific elevation of genetic risk of suicide in this subgroup, and also showed increased PRS of PTSD, suicide attempt, and risk taking. LFR suicides were substantially older at death (mean = 49.10 years), had fewer psychiatric diagnoses of depression and pain, and significantly lower PRS of depression. Results suggest extended familiality and trauma/PTSD may provide specificity in identifying individuals at genetic risk for suicide death, especially among younger ages, and that LFR of suicide warrants further study regarding the contribution of demographic and medical risks.

Methylome-wide association findings for major depressive disorder overlap in blood and brain and replicate in independent brain samples.

We present the first large-scale methylome-wide association studies (MWAS) for major depressive disorder (MDD) to identify sites of potential importance for MDD etiology. Using a sequencing-based approach that provides near-complete coverage of all 28 million common CpGs in the human genome, we assay methylation in MDD cases and controls from both blood (N = 1132) and postmortem brain tissues (N = 61 samples from Brodmann Area 10, BA10). The MWAS for blood identified several loci with P ranging from 1.91 × 10-8 to 4.39 × 10-8 and a resampling approach showed that the cumulative association was significant (P = 4.03 × 10-10) with the signal coming from the top 25,000 MWAS markers. Furthermore, a permutation-based analysis showed significant overlap (P = 5.4 × 10-3) between the MWAS findings in blood and brain (BA10). This overlap was significantly enriched for a number of features including being in eQTLs in blood and the frontal cortex, CpG islands and shores, and exons. The overlapping sites were also enriched for active chromatin states in brain including genic enhancers and active transcription start sites. Furthermore, three loci located in GABBR2, RUFY3, and in an intergenic region on chromosome 2 replicated with the same direction of effect in the second brain tissue (BA25, N = 60) from the same individuals and in two independent brain collections (BA10, N = 81 and 64). GABBR2 inhibits neuronal activity through G protein-coupled second-messenger systems and RUFY3 is implicated in the establishment of neuronal polarity and axon elongation. In conclusion, we identified and replicated methylated loci associated with MDD that are involved in biological functions of likely importance to MDD etiology.

A methylation study of long-term depression risk.

Recurrent and chronic major depressive disorder (MDD) accounts for a substantial part of the disease burden because this course is most prevalent and typically requires long-term treatment. We associated blood DNA methylation profiles from 581 MDD patients at baseline with MDD status 6 years later. A resampling approach showed a highly significant association between methylation profiles in blood at baseline and future disease status (P = 2.0 × 10-16). Top MWAS results were enriched specific pathways, overlapped with genes found in GWAS of MDD disease status, autoimmune disease and inflammation, and co-localized with eQTLS and (genic enhancers of) of transcription sites in brain and blood. Many of these findings remained significant after correction for multiple testing. The major themes emerging were cellular responses to stress and signaling mechanisms linked to immune cell migration and inflammation. This suggests that an immune signature of treatment-resistant depression is already present at baseline. We also created a methylation risk score (MRS) to predict MDD status 6 years later. The AUC of our MRS was 0.724 and higher than risk scores created using a set of five putative MDD biomarkers, genome-wide SNP data, and 27 clinical, demographic and lifestyle variables. Although further studies are needed to examine the generalizability to different patient populations, these results suggest that methylation profiles in blood may present a promising avenue to support clinical decision making by providing empirical information about the likelihood MDD is chronic or will recur in the future.

Genome-wide significant regions in 43 Utah high-risk families implicate multiple genes involved in risk for completed suicide.

Suicide is the 10th leading cause of death in the United States. Although environment has undeniable impact, evidence suggests that genetic factors play a significant role in completed suicide. We linked a resource of ~ 4500 DNA samples from completed suicides obtained from the Utah Medical Examiner to genealogical records and medical records data available on over eight million individuals. This linking has resulted in the identification of high-risk extended families (7-9 generations) with significant familial risk of completed suicide. Familial aggregation across distant relatives minimizes effects of shared environment, provides more genetically homogeneous risk groups, and magnifies genetic risks through familial repetition. We analyzed Illumina PsychArray genotypes from suicide cases in 43 high-risk families, identifying 30 distinct shared genomic segments with genome-wide evidence (p = 2.02E-07-1.30E-18) of segregation with completed suicide. The 207 genes implicated by the shared regions provide a focused set of genes for further study; 18 have been previously associated with suicide risk. Although PsychArray variants do not represent exhaustive variation within the 207 genes, we investigated these for specific segregation within the high-risk families, and for association of variants with predicted functional impact in ~ 1300 additional Utah suicides unrelated to the discovery families. None of the limited PsychArray variants explained the high-risk family segregation; sequencing of these regions will be needed to discover segregating risk variants, which may be rarer or regulatory. However, additional association tests yielded four significant PsychArray variants (SP110, rs181058279; AGBL2, rs76215382; SUCLA2, rs121908538; APH1B, rs745918508), raising the likelihood that these genes confer risk of completed suicide.

Exploring the genetic overlap of suicide-related behaviors and substance use disorders.

Suicide-related behaviors are heterogeneous and transdiagnostic, and may demonstrate varying levels of genetic overlap with different substance use disorders (SUDs). We used linkage disequilibrium score regression, genomic structural equation models, and Mendelian randomization to examine the genetic relationships between several SUDs and suicide-related behaviors. Our analyses incorporated summary statistics from the largest genome-wide association studies (GWAS) of problematic alcohol use, the Fagerström test for nicotine dependence, cannabis use disorder, and opioid use disorder (Ns ranging from 46,213-435,563) and GWAS of ever self-harmed, suicide attempt, and suicide death (Ns ranging from 18,223-117,733). We also accounted for genetic liability to depression (N = 500,199) and risk tolerance (N = 315,894). Suicide-related behaviors were significantly genetically correlated with each other and each SUD, but there was little evidence of causal relationships between the traits. Simultaneously correlating a common SUD factor with each specific suicide indicator while controlling for depression and risk tolerance revealed significant, positive genetic correlations between the SUD factor and suicide-related behaviors (rg  = 0.26-0.45, SE = 0.08-0.09). In the model, depression's association with suicide death (ß = 0.42, SE = 0.06) was weaker compared to ever-self harmed and suicide attempt (ß = 0.58, SE = 0.05 and ß = 0.50, SE = 0.06, respectively). We identify a general level of genetic overlap between SUDs and suicide-related behaviors, which is independent of depression and risk tolerance. Additionally, our findings suggest that genetic and behavioral contributions to suicide death may somewhat differ from nonlethal suicide-related behaviors.

Rare protein-coding variants implicate genes involved in risk of suicide death.

Identification of genetic factors leading to increased risk of suicide death is critical to combat rising suicide rates, however, only a fraction of the genetic variation influencing risk has been accounted for. To address this limitation, we conducted the first comprehensive analysis of rare genetic variation in suicide death leveraging the largest suicide death biobank, the Utah Suicide Genetic Risk Study (USGRS). We conducted a single-variant association analysis of rare (minor allele frequency <1%) putatively functional single-nucleotide polymorphisms (SNPs) present on the Illumina PsychArray genotyping array in 2,672 USGRS suicide deaths of non-Finnish European (NFE) ancestry and 51,583 NFE controls from the Genome Aggregation Database. Secondary analyses used an independent control sample of 21,324 NFE controls from the Psychiatric Genomics Consortium. Five novel, high-impact, rare SNPs were identified with significant associations with suicide death (SNAPC1, rs75418419; TNKS1BP1, rs143883793; ADGRF5, rs149197213; PER1, rs145053802; and ESS2, rs62223875). 119 suicide decedents carried these high-impact SNPs. Both PER1 and SNAPC1 have other supporting gene-level evidence of suicide risk, and psychiatric associations exist for PER1 (bipolar disorder, schizophrenia), and for TNKS1BP1 and ESS2 (schizophrenia). Three of the genes (PER1, TNKS1BP1, and ADGRF5), together with additional genes implicated by genome-wide association studies on suicidal behavior, showed significant enrichment in immune system, homeostatic and signal transduction processes. No specific diagnostic phenotypes were associated with the subset of suicide deaths with the identified rare variants. These findings suggest an important role for rare variants in suicide risk and implicate genes and gene pathways for targeted replication.

Building a schizophrenia genetic network: transcription factor 4 regulates genes involved in neuronal development and schizophrenia risk.

The transcription factor 4 (TCF4) locus is a robust association finding with schizophrenia (SCZ), but little is known about the genes regulated by the encoded transcription factor. Therefore, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) of TCF4 in neural-derived (SH-SY5Y) cells to identify genome-wide TCF4 binding sites, followed by data integration with SCZ association findings. We identified 11 322 TCF4 binding sites overlapping in two ChIP-seq experiments. These sites are significantly enriched for the TCF4 Ebox binding motif (>85% having ≥1 Ebox) and implicate a gene set enriched for genes downregulated in TCF4 small-interfering RNA (siRNA) knockdown experiments, indicating the validity of our findings. The TCF4 gene set was also enriched among (1) gene ontology categories such as axon/neuronal development, (2) genes preferentially expressed in brain, in particular pyramidal neurons of the somatosensory cortex and (3) genes downregulated in postmortem brain tissue from SCZ patients (odds ratio, OR = 2.8, permutation P < 4x10-5). Considering genomic alignments, TCF4 binding sites significantly overlapped those for neural DNA-binding proteins such as FOXP2 and the SCZ-associated EP300. TCF4 binding sites were modestly enriched among SCZ risk loci from the Psychiatric Genomic Consortium (OR = 1.56, P = 0.03). In total, 130 TCF4 binding sites occurred in 39 of the 108 regions published in 2014. Thirteen genes within the 108 loci had both a TCF4 binding site ±10kb and were differentially expressed in siRNA knockdown experiments of TCF4, suggesting direct TCF4 regulation. These findings confirm TCF4 as an important regulator of neural genes and point toward functional interactions with potential relevance for SCZ.

Polygenic prediction of PTSD trajectories in 9/11 responders.

BACKGROUND: Genetics hold promise of predicting long-term post-traumatic stress disorder (PTSD) outcomes following trauma. The aim of the current study was to test whether six hypothesized polygenic risk scores (PRSs) developed to capture genetic vulnerability to psychiatric conditions prospectively predict PTSD onset, severity, and 18-year course after trauma exposure. METHODS: Participants were 1490 responders to the World Trade Center (WTC) disaster (mean age at 9/11 = 38.81 years, s.d. = 8.20; 93.5% male; 23.8% lifetime WTC-related PTSD diagnosis). Prospective longitudinal data on WTC-related PTSD symptoms were obtained from electronic medical records and modelled as PTSD trajectories using growth mixture model analysis. Independent regression models tested whether six hypothesized psychiatric PRSs (PTSD-PRS, Re-experiencing-PRS, Generalized Anxiety-PRS, Schizophrenia-PRS, Depression-PRS, and Neuroticism-PRS) are predictive of WTC-PTSD outcomes: lifetime diagnoses, average symptom severity, and 18-year symptom trajectory. All analyses were adjusted for population stratification, 9/11 exposure severity, and multiple testing. RESULTS: Depression-PRS predicted PTSD diagnostic status (OR 1.37, CI 1.17-1.61, adjusted p = 0.001). All PRSs, except PTSD-PRS, significantly predicted average PTSD symptoms (ß = 0.06-0.10, adjusted p < 0.05). Re-experiencing-PRS, Generalized Anxiety-PRS and Schizophrenia-PRS predicted the high severity PTSD trajectory class (ORs 1.21-1.28, adjusted p < 0.05). Finally, PRSs prediction was independent of 9/11 exposure severity and jointly accounted for 3.7 times more variance in PTSD symptoms than the exposure severity. CONCLUSIONS: Psychiatric PRSs prospectively predicted WTC-related PTSD lifetime diagnosis, average symptom severity, and 18-year trajectory in responders to 9/11 disaster. Jointly, PRSs were more predictive of subsequent PTSD than the exposure severity. In the future, PRSs may help identify at-risk responders who might benefit from targeted prevention approaches.

An empirical Bayes approach for multiple tissue eQTL analysis.

Expression quantitative trait locus (eQTL) analyses identify genetic markers associated with the expression of a gene. Most up-to-date eQTL studies consider the connection between genetic variation and expression in a single tissue. Multi-tissue analyses have the potential to improve findings in a single tissue, and elucidate the genotypic basis of differences between tissues. In this article, we develop a hierarchical Bayesian model (MT-eQTL) for multi-tissue eQTL analysis. MT-eQTL explicitly captures patterns of variation in the presence or absence of eQTL, as well as the heterogeneity of effect sizes across tissues. We devise an efficient Expectation-Maximization (EM) algorithm for model fitting. Inferences concerning eQTL detection and the configuration of eQTL across tissues are derived from the adaptive thresholding of local false discovery rates, and maximum a posteriori estimation, respectively. We also provide theoretical justification of the adaptive procedure. We investigate the MT-eQTL model through an extensive analysis of a 9-tissue data set from the GTEx initiative.

RaMWAS: fast methylome-wide association study pipeline for enrichment platforms.

Motivation: Enrichment-based technologies can provide measurements of DNA methylation at tens of millions of CpGs for thousands of samples. Existing tools for methylome-wide association studies cannot analyze datasets of this size and lack important features like principal component analysis, combined analysis with SNP data and outcome predictions that are based on all informative methylation sites. Results: We present a Bioconductor R package called RaMWAS with a full set of tools for large-scale methylome-wide association studies. It is free, cross-platform, open source, memory efficient and fast. Availability and implementation: Release version and vignettes with small case study at bioconductor.org/packages/ramwas Development version at github.com/andreyshabalin/ramwas. Supplementary information: Supplementary data are available at Bioinformatics online.

Enrichment methods provide a feasible approach to comprehensive and adequately powered investigations of the brain methylome.

Methylome-wide association studies are typically performed using microarray technologies that only assay a very small fraction of the CG methylome and entirely miss two forms of methylation that are common in brain and likely of particular relevance for neuroscience and psychiatric disorders. The alternative is to use whole genome bisulfite (WGB) sequencing but this approach is not yet practically feasible with sample sizes required for adequate statistical power. We argue for revisiting methylation enrichment methods that, provided optimal protocols are used, enable comprehensive, adequately powered and cost-effective genome-wide investigations of the brain methylome. To support our claim we use data showing that enrichment methods approximate the sensitivity obtained with WGB methods and with slightly better specificity. However, this performance is achieved at <5% of the reagent costs. Furthermore, because many more samples can be sequenced simultaneously, projects can be completed about 15 times faster. Currently the only viable option available for comprehensive brain methylome studies, enrichment methods may be critical for moving the field forward.

Estimation of cis-eQTL effect sizes using a log of linear model.

The study of expression Quantitative Trait Loci (eQTL) is an important problem in genomics and biomedicine. While detection (testing) of eQTL associations has been widely studied, less work has been devoted to the estimation of eQTL effect size. To reduce false positives, detection methods frequently rely on linear modeling of rank-based normalized or log-transformed gene expression data. Unfortunately, these approaches do not correspond to the simplest model of eQTL action, and thus yield estimates of eQTL association that can be uninterpretable and inaccurate. In this article, we propose a new, log-of-linear model for eQTL action, termed ACME, that captures allelic contributions to cis-acting eQTLs in an additive fashion, yielding effect size estimates that correspond to a biologically coherent model of cis-eQTLs. We describe a non-linear least-squares algorithm to fit the model by maximum likelihood, and obtain corresponding p-values. We perform careful investigation of the model using a combination of simulated data and data from the Genotype Tissue Expression (GTEx) project. Our results reveal little evidence for dominance effects, a parsimonious result that accords with a simple biological model for allele-specific expression and supports use of the ACME model. We show that Type-I error is well-controlled under our approach in a realistic setting, so that rank-based normalizations are unnecessary. Furthermore, we show that such normalizations can be detrimental to power and estimation accuracy under the proposed model. We then show, through effect size analyses of whole-genome cis-eQTLs in the GTEx data, that using standard normalizations instead of ACME noticeably affects the ranking and sign of estimates.

A Whole Methylome Study of Ethanol Exposure in Brain and Blood: An Exploration of the Utility of Peripheral Blood as Proxy Tissue for Brain in Alcohol Methylation Studies.

BACKGROUND: Recent reviews have highlighted the potential use of blood-based methylation biomarkers as diagnostic and prognostic tools of current and future alcohol use and addiction. Due to the substantial overlap that often exists between methylation patterns across different tissues, including blood and brain, blood-based methylation may track methylation changes in brain; however, little work has explored the overlap in alcohol-related methylation in these tissues. METHODS: To study the effects of alcohol on the brain methylome and identify possible biomarkers of these changes in blood, we performed a methylome-wide association study in brain and blood from 40 male DBA/2J mice that received either an acute ethanol (EtOH) or saline intraperitoneal injection. To investigate all 22 million CpGs in the mouse genome, we enriched for the methylated genomic fraction using methyl-CpG binding domain (MBD) protein capture followed by next-generation sequencing (MBD-seq). We performed association tests in blood and brain separately followed by enrichment testing to determine whether there was overlapping alcohol-related methylation in the 2 tissues. RESULTS: The top result for brain was a CpG located in an intron of Ttc39b (p = 5.65 × 10-08 ), and for blood, the top result was located in Espnl (p = 5.11 × 10-08 ). Analyses implicated pathways involved in inflammation and neuronal differentiation, such as CXCR4, IL-7, and Wnt signaling. Enrichment tests indicated significant overlap among the top results in brain and blood. Pathway analyses of the overlapping genes converge on MAPKinase signaling (p = 5.6 × 10-05 ) which plays a central role in acute and chronic responses to alcohol and glutamate receptor pathways, which can regulate neuroplastic changes underlying addictive behavior. CONCLUSIONS: Overall, we have shown some methylation changes in brain and blood after acute EtOH administration and that the changes in blood partly mirror the changes in brain suggesting the potential for DNA methylation in blood to be biomarkers of alcohol use.