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.

Random allelic expression in the adult human body.

Genes are typically assumed to express both parental alleles similarly, yet cell lines show random allelic expression (RAE) for many autosomal genes that could shape genetic effects. Thus, understanding RAE in human tissues could improve our understanding of phenotypic variation. Here, we develop a methodology to perform genome-wide profiling of RAE and biallelic expression in GTEx datasets for 832 people and 54 tissues. We report 2,762 autosomal genes with some RAE properties similar to randomly inactivated X-linked genes. We found that RAE is associated with rapidly evolving regions in the human genome, adaptive signaling processes, and genes linked to age-related diseases such as neurodegeneration and cancer. We define putative mechanistic subtypes of RAE distinguished by gene overlaps on sense and antisense DNA strands, aggregation in clusters near telomeres, and increased regulatory complexity and inputs compared with biallelic genes. We provide foundations to study RAE in human phenotypes, evolution, and disease.

Noncanonical genomic imprinting in the monoamine system determines naturalistic foraging and brain-adrenal axis functions.

Noncanonical genomic imprinting can cause biased expression of one parental allele in a tissue; however, the functional relevance of such biases is unclear. To investigate ethological roles for noncanonical imprinting in dopa decarboxylase (Ddc) and tyrosine hydroxylase (Th), we use machine learning to decompose naturalistic foraging in maternal and paternal allele mutant heterozygous mice. We uncover distinct roles for the maternal versus paternal alleles on foraging, where maternal alleles affect sons while daughters are under paternal allelic control. Each parental allele controls specific action sequences reflecting decisions in naive or familiar contexts. The maternal Ddc allele is preferentially expressed in subsets of hypothalamic GABAergic neurons, while the paternal allele predominates in subsets of adrenal cells. Each Ddc allele affects distinct molecular and endocrine components of the brain-adrenal axis. Thus, monoaminergic noncanonical imprinting has ethological roles in foraging and endocrine functions and operates by affecting discrete subsets of cells.

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.

Elephant Genomes Reveal Accelerated Evolution in Mechanisms Underlying Disease Defenses.

Disease susceptibility and resistance are important factors for the conservation of endangered species, including elephants. We analyzed pathology data from 26 zoos and report that Asian elephants have increased neoplasia and malignancy prevalence compared with African bush elephants. This is consistent with observed higher susceptibility to tuberculosis and elephant endotheliotropic herpesvirus (EEHV) in Asian elephants. To investigate genetic mechanisms underlying disease resistance, including differential responses between species, among other elephant traits, we sequenced multiple elephant genomes. We report a draft assembly for an Asian elephant, and defined 862 and 1,017 conserved potential regulatory elements in Asian and African bush elephants, respectively. In the genomes of both elephant species, conserved elements were significantly enriched with genes differentially expressed between the species. In Asian elephants, these putative regulatory regions were involved in immunity pathways including tumor-necrosis factor, which plays an important role in EEHV response. Genomic sequences of African bush, forest, and Asian elephant genomes revealed extensive sequence conservation at TP53 retrogene loci across three species, which may be related to TP53 functionality in elephant cancer resistance. Positive selection scans revealed outlier genes related to additional elephant traits. Our study suggests that gene regulation plays an important role in the differential inflammatory response of Asian and African elephants, leading to increased infectious disease and cancer susceptibility in Asian elephants. These genomic discoveries can inform future functional and translational studies aimed at identifying effective treatment approaches for ill elephants, which may improve conservation.

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.

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.

Parallel Accelerated Evolution in Distant Hibernators Reveals Candidate Cis Elements and Genetic Circuits Regulating Mammalian Obesity.

Obesity is a clinical problem and an important adaptation in many species. Hibernating mammals, for example, become obese, insulin resistant, and hyperinsulinemic to store fat. Here, we combine comparative phylogenomics with large-scale human genome data to uncover candidate cis elements regulating mammalian obesity. Our study examines genetic elements conserved across non-hibernating mammals to identify genome-wide patterns of accelerated evolution in hibernators from different clades. The results reveal the existence of parallel accelerated regions (pARs) in distant hibernators. Hibernator pARs are disproportionately located near human obesity susceptibility genes compared to random conserved regions, hibernator ARs that are not parallel, and non-hibernator pARs. We found 364 candidate obesity-regulating cis elements and genetic circuits in different cell types. The Fat Mass and Obesity (FTO) locus, the strongest genetic risk factor for human obesity, is an enriched site for hibernator pARs. Our results uncover noncoding cis elements with putative roles in obesity and hibernation.

Complex Economic Behavior Patterns Are Constructed from Finite, Genetically Controlled Modules of Behavior.

Complex ethological behaviors could be constructed from finite modules that are reproducible functional units of behavior. Here, we test this idea for foraging and develop methods to dissect rich behavior patterns in mice. We uncover discrete modules of foraging behavior reproducible across different strains and ages, as well as nonmodular behavioral sequences. Modules differ in terms of form, expression frequency, and expression timing and are expressed in a probabilistically determined order. Modules shape economic patterns of feeding, exposure, activity, and perseveration responses. The modular architecture of foraging changes developmentally, and different developmental, genetic, and parental effects are found to shape the expression of specific modules. Dissecting modules from complex patterns is powerful for phenotype analysis. We discover that both parental alleles of the imprinted Prader-Willi syndrome gene Magel2 are functional in mice but regulate different modules. Our study found that complex economic patterns are built from finite, genetically controlled modules.

Accelerated Evolution in Distinctive Species Reveals Candidate Elements for Clinically Relevant Traits, Including Mutation and Cancer Resistance.

The identity of most functional elements in the mammalian genome and the phenotypes they impact are unclear. Here, we perform a genome-wide comparative analysis of patterns of accelerated evolution in species with highly distinctive traits to discover candidate functional elements for clinically important phenotypes. We identify accelerated regions (ARs) in the elephant, hibernating bat, orca, dolphin, naked mole rat, and thirteen-lined ground squirrel lineages in mammalian conserved regions, uncovering ∼33,000 elements that bind hundreds of different regulatory proteins in humans and mice. ARs in the elephant, the largest land mammal, are uniquely enriched near elephant DNA damage response genes. The genomic hotspot for elephant ARs is the E3 ligase subunit of the Fanconi anemia complex, a master regulator of DNA repair. Additionally, ARs in the six species are associated with specific human clinical phenotypes that have apparent concordance with overt traits in each species.

Diverse Non-genetic, Allele-Specific Expression Effects Shape Genetic Architecture at the Cellular Level in the Mammalian Brain.

Interactions between genetic and epigenetic effects shape brain function, behavior, and the risk for mental illness. Random X inactivation and genomic imprinting are epigenetic allelic effects that are well known to influence genetic architecture and disease risk. Less is known about the nature, prevalence, and conservation of other potential epigenetic allelic effects in vivo in the mouse and primate brain. Here we devise genomics, in situ hybridization, and mouse genetics strategies to uncover diverse allelic effects in the brain that are not caused by imprinting or genetic variation. We found allelic effects that are developmental stage and cell type specific, that are prevalent in the neonatal brain, and that cause mosaics of monoallelic brain cells that differentially express wild-type and mutant alleles for heterozygous mutations. Finally, we show that diverse non-genetic allelic effects that impact mental illness risk genes exist in the macaque and human brain. Our findings have potential implications for mammalian brain genetics. VIDEO ABSTRACT.

Noncanonical Genomic Imprinting Effects in Offspring.

Here, we describe an RNA-sequencing (RNA-seq)-based approach that accurately detects even modest maternal or paternal allele expression biases at the tissue level, which we call noncanonical genomic imprinting effects. We profile imprinting in the arcuate nucleus (ARN) and dorsal raphe nucleus of the female mouse brain as well as skeletal muscle (mesodermal) and liver (endodermal). Our study uncovers hundreds of noncanonical autosomal and X-linked imprinting effects. Noncanonical imprinting is highly tissue-specific and enriched in the ARN, but rare in the liver. These effects are reproducible across different genetic backgrounds and associated with allele-specific chromatin. Using in situ hybridization for nascent RNAs, we discover that autosomal noncanonical imprinted genes with a tissue-level allele bias exhibit allele-specific expression effects in subpopulations of neurons in the brain in vivo. We define noncanonical imprinted genes that regulate monoamine signaling and determine that these effects influence the impact of inherited mutations on offspring behavior.

Selective Oma1 protease-mediated proteolysis of Cox1 subunit of cytochrome oxidase in assembly mutants.

Stalled biogenesis of the mitochondrial cytochrome c oxidase (CcO) complex results in degradation of subunits containing redox cofactors. The conserved Oma1 metalloproteinase mediates facile Cox1 degradation in cells lacking the Coa2 assembly factor, but not in a series of other mutants stalled in CcO maturation. Oma1 is activated in coa2Δ cells, but the selective Cox1 degradation does not arise merely from its activation. Oma1 is also active in cells with dysfunctional mitochondria and cox11Δ cells impaired in CcO maturation, but this activation does not result in Oma1-mediated Cox1 degradation. The facile and selective degradation of Cox1 in coa2Δ cells, relative to other CcO assembly mutants, is likely due to impaired hemylation and subsequent misfolding of the subunit. Specific Cox1 proteolysis in coa2Δ cells arises from a combination of Oma1 activation and a susceptible conformation of Cox1.

Autoregulation of the rsc4 tandem bromodomain by gcn5 acetylation.

An important issue for chromatin remodeling complexes is how their bromodomains recognize particular acetylated lysine residues in histones. The Rsc4 subunit of the yeast remodeler RSC contains an essential tandem bromodomain (TBD) that binds acetylated K14 of histone H3 (H3K14ac). We report a series of crystal structures that reveal a compact TBD that binds H3K14ac in the second bromodomain and, remarkably, binds acetylated K25 of Rsc4 itself in the first bromodomain. Endogenous Rsc4 is acetylated only at K25, and Gcn5 is identified as necessary and sufficient for Rsc4 K25 acetylation in vivo and in vitro. Rsc4 K25 acetylation inhibits binding to H3K14ac, and mutation of Rsc4 K25 results in altered growth rates. These data suggest an autoregulatory mechanism in which Gcn5 performs both the activating (H3K14ac) and inhibitory (Rsc4 K25ac) modifications, perhaps to provide temporal regulation. Additional regulatory mechanisms are indicated as H3S10 phosphorylation inhibits Rsc4 binding to H3K14ac peptides.

The structure of the yFACT Pob3-M domain, its interaction with the DNA replication factor RPA, and a potential role in nucleosome deposition.

We report the crystal structure of the middle domain of the Pob3 subunit (Pob3-M) of S. cerevisiae FACT (yFACT, facilitates chromatin transcription), which unexpectedly adopts an unusual double pleckstrin homology (PH) architecture. A mutation within a conserved surface cluster in this domain causes a defect in DNA replication that is suppressed by mutation of replication protein A (RPA). The nucleosome reorganizer yFACT therefore interacts in a physiologically important way with the central single-strand DNA (ssDNA) binding factor RPA to promote a step in DNA replication. Purified yFACT and RPA display a weak direct physical interaction, although the genetic suppression is not explained by simple changes in affinity between the purified proteins. Further genetic analysis suggests that coordinated function by yFACT and RPA is important during nucleosome deposition. These results support the model that the FACT family has an essential role in constructing nucleosomes during DNA replication, and suggest that RPA contributes to this process.