Neurexin 1 variants as risk factors for suicide death.

Suicide is a significant public health concern with complex etiology. Although the genetic component of suicide is well established, the scope of gene networks and biological mechanisms underlying suicide has yet to be defined. Previously, we reported genome-wide evidence that neurexin 1 (NRXN1), a key synapse organizing molecule, is associated with familial suicide risk. Here we present new evidence for two non-synonymous variants (rs78540316; P469S and rs199784139; H885Y) associated with increased familial risk of suicide death. We tested the impact of these variants on binding interactions with known partners and assessed functionality in a hemi-synapse formation assay. Although the formation of hemi-synapses was not altered with the P469S variant relative to wild-type, both variants increased binding to the postsynaptic binding partner, leucine-rich repeat transmembrane neuronal 2 (LRRTM2) in vitro. Our findings indicate that variants in NRXN1 and related synaptic genes warrant further study as risk factors for suicide death.

Genome-Wide Association Study of Suicide Death and Polygenic Prediction of Clinical Antecedents.

OBJECTIVE: Death by suicide is a highly preventable yet growing worldwide health crisis. To date, there has been a lack of adequately powered genomic studies of suicide, with no sizable suicide death cohorts available for analysis. To address this limitation, the authors conducted the first comprehensive genomic analysis of suicide death using previously unpublished genotype data from a large population-ascertained cohort. METHODS: The analysis sample comprised 3,413 population-ascertained case subjects of European ancestry and 14,810 ancestrally matched control subjects. Analytical methods included principal component analysis for ancestral matching and adjusting for population stratification, linear mixed model genome-wide association testing (conditional on genetic-relatedness matrix), gene and gene set-enrichment testing, and polygenic score analyses, as well as single-nucleotide polymorphism (SNP) heritability and genetic correlation estimation using linkage disequilibrium score regression. RESULTS: Genome-wide association analysis identified two genome-wide significant loci (involving six SNPs: rs34399104, rs35518298, rs34053895, rs66828456, rs35502061, and rs35256367). Gene-based analyses implicated 22 genes on chromosomes 13, 15, 16, 17, and 19 (q<0.05). Suicide death heritability was estimated at an h2SNP value of 0.25 (SE=0.04) and a value of 0.16 (SE=0.02) when converted to a liability scale. Notably, suicide polygenic scores were significantly predictive across training and test sets. Polygenic scores for several other psychiatric disorders and psychological traits were also predictive, particularly scores for behavioral disinhibition and major depressive disorder. CONCLUSIONS: Multiple genome-wide significant loci and genes were identified and polygenic score prediction of suicide death case-control status was demonstrated, adjusting for ancestry, in independent training and test sets. Additionally, the suicide death sample was found to have increased genetic risk for behavioral disinhibition, major depressive disorder, depressive symptoms, autism spectrum disorder, psychosis, and alcohol use disorder compared with the control sample.

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.

Meta-analytic approaches for examining complexity and heterogeneity in studies of adolescent development.

INTRODUCTION: In the field of adolescent development, meta-analysis offers valuable tools for synthesizing and assessing cumulative research evidence on the effectiveness of programs, practices, and policies intended to promote healthy adolescent development. When examining the impact of a program implemented across multiple primary studies, variation is often observed in the methodological attributes of those primary studies, such as their implementation methods, program components, participant characteristics, outcome measurement, and the systems in which programs are deployed. Differences in methodological attributes of primary studies represented in a meta-analysis, referred to as complexity, can yield variation in true effects across primary studies, which is described as heterogeneity. METHODS: We discuss heterogeneity as a parameter of interest in meta-analysis, introducing and demonstrating both graphical and statistical methods for evaluating the magnitude and impact of heterogeneity. We discuss approaches for presenting characteristics of heterogeneity in meta-analytic findings, and methods for identifying and statistically controlling for aspects of methodological complexity that may contribute to variation in effects across primary studies. RESULTS: Topics and methods related to assessing and explaining heterogeneity were contextualized in the field of adolescent development using a sample of primary studies from a large meta-analysis examining the effectiveness of brief alcohol interventions for youth. We highlighted approaches currently underutilized in the field and provided R code for key methods to broaden their use. CONCLUSIONS: By discussing various heterogeneity statistics, visualizations, and explanatory methods, this article provides the applied developmental researcher a foundational understanding of complexity and heterogeneity in meta-analysis.

Novel pedigree analysis implicates DNA repair and chromatin remodeling in multiple myeloma risk.

The high-risk pedigree (HRP) design is an established strategy to discover rare, highly-penetrant, Mendelian-like causal variants. Its success, however, in complex traits has been modest, largely due to challenges of genetic heterogeneity and complex inheritance models. We describe a HRP strategy that addresses intra-familial heterogeneity, and identifies inherited segments important for mapping regulatory risk. We apply this new Shared Genomic Segment (SGS) method in 11 extended, Utah, multiple myeloma (MM) HRPs, and subsequent exome sequencing in SGS regions of interest in 1063 MM / MGUS (monoclonal gammopathy of undetermined significance-a precursor to MM) cases and 964 controls from a jointly-called collaborative resource, including cases from the initial 11 HRPs. One genome-wide significant 1.8 Mb shared segment was found at 6q16. Exome sequencing in this region revealed predicted deleterious variants in USP45 (p.Gln691* and p.Gln621Glu), a gene known to influence DNA repair through endonuclease regulation. Additionally, a 1.2 Mb segment at 1p36.11 is inherited in two Utah HRPs, with coding variants identified in ARID1A (p.Ser90Gly and p.Met890Val), a key gene in the SWI/SNF chromatin remodeling complex. Our results provide compelling statistical and genetic evidence for segregating risk variants for MM. In addition, we demonstrate a novel strategy to use large HRPs for risk-variant discovery more generally in complex traits.

Evidence for Association Between Low Frequency Variants in CHRNA6/CHRNB3 and Antisocial Drug Dependence.

Common SNPs in nicotinic acetylcholine receptor genes (CHRN genes) have been associated with drug behaviors and personality traits, but the influence of rare genetic variants is not well characterized. The goal of this project was to identify novel rare variants in CHRN genes in the Center for Antisocial Drug Dependence (CADD) and Genetics of Antisocial Drug Dependence (GADD) samples and to determine if low frequency variants are associated with antisocial drug dependence. Two samples of 114 and 200 individuals were selected using a case/control design including the tails of the phenotypic distribution of antisocial drug dependence. The capture, sequencing, and analysis of all variants in 16 CHRN genes (CHRNA1-7, 9, 10, CHRNB1-4, CHRND, CHRNG, CHRNE) were performed independently for each subject in each sample. Sequencing reads were aligned to the human reference sequence using BWA prior to variant calling with the Genome Analysis ToolKit (GATK). Low frequency variants (minor allele frequency < 0.05) were analyzed using SKAT-O and C-alpha to examine the distribution of rare variants among cases and controls. In our larger sample, the region containing the CHRNA6/CHRNB3 gene cluster was significantly associated with disease status using both SKAT-O and C-alpha (unadjusted p values <0.05). More low frequency variants in the CHRNA6/CHRNB3 gene region were observed in cases compared to controls. These data support a role for genetic variants in CHRN genes and antisocial drug behaviors.

Mesolimbic transcriptional response to hedonic substitution of voluntary exercise and voluntary ethanol consumption.

The mesolimbic dopaminergic pathway has been implicated in many rewarding behaviors, including the consumption of ethanol and voluntary exercise. It has become apparent that different rewarding stimuli activate this pathway, and therefore it is possible for these behaviors to influence each other, i.e. hedonic substitution. Using adult female C57BL/6J mice, we demonstrate that voluntary access to a running wheel substantially reduces the consumption and preference of ethanol. Furthermore, we examined gene expression of several genes involved in regulating the mesolimbic dopaminergic pathway, which we hypothesized to be the main pathway involved in hedonic substitution. In the striatum, we observed a reduction in mRNA expression of Drd1a due to exercise. Hippocampal Bdnf mRNA increased in response to exercise and decreased in response to ethanol. Furthermore, there was an interaction effect of exercise and ethanol on the expression of Slc18a2 in the midbrain. These data suggest an important role for this pathway, and especially for Bdnf and Slc18a2 in regulating hedonic substitution.