Somatic deletions implicated in functional diversity of brain cells of individuals with schizophrenia and unaffected controls

The gain-of-function UBE3AQ588E variant causes Angelman-like neurodevelopmental phenotypes in mice

Mutations in the E3 ubiquitin ligase UBE3A that cause enzymatic gain-of-function result in disease phenotypes which differ from classic Angelman syndrome. However, these phenotypes are highly heterogeneous raising questions about the mechanistic basis of such phenotypic diversity. Here, we characterize a mouse model harboring a Ube3aQ606E gain of function variant (UBE3AQ588E in humans). Extensive behavioral phenotyping showed that animals possessing a maternally inherited mutation (Ube3amQ606E) paradoxically show many behavioral deficits indicative of overall UBE3A loss-of-function. These included pronounced motor deficits, hypoactivity, and reduced stereotypic behaviors. Moreover, brain weights and MRI analysis revealed global microcephaly with a postnatal onset, consistent with phenotypes described in Angelman syndrome model mice. Additional biochemical analyses demonstrated an increased abundance of UBE3A substrates in brain tissue and immunofluorescence analyses showed that microcephaly is not caused by increased apoptotic cell death. Together, our results strongly suggest a novel mechanism by which the Ube3amQ606E mutation leads to enhanced self-targeted degradation of UBE3A, leading to an overall loss of enzyme activity, resulting in Angelman-like phenotypes in vivo.

Translational Approach to Social Isolation During a Global Pandemic: Hippocampal Somatic Mutation and Stress

OBJECTIVE

The coronavirus disease-2019 (COVID-19) pandemic's social isolation has significantly impacted mental health, increasing depression and anxiety. This study explores the effects of social isolation on both humans and mice, focusing on behavioral changes and hippocampal protein expression. It also investigates genetic alterations through single-cell RNA and whole-genome sequencing (WGS).

METHODS

Here we conducted behavioral studies, protein expression studies, single-nucleus sequencing (snRNAseq), and WGS of the hippocampus of mice that underwent early maternal separation and social isolation, and a demographic study of community populations who had been self-quarantined owing to COVID-19 exposure to investigate the link between somatic mutations and stress due to social isolation.

RESULTS

The demographic study demonstrated more negative mental health findings among individuals who live alone or are single. Mice subjected to early maternal separation and social isolation demonstrated increased anxiety-like behaviors and stress-related corticotropin-releasing hormone receptor 1, and neurogenesis-related sex-determining region Y-box 2 and doublecortin expression. In snRNA-seq, differences, such as transthyretin increase, were observed in the maternal separation group, and somatic mutations, including insertion in the intron site of Tmem267, were observed in the social isolation group on WGS.

CONCLUSION

The results of this study suggest that stress, such as social isolation, can cause changes at the genetic level, as well as behavioral and brain protein changes.

Chromosome instability and aneuploidy in the mammalian brain

This review investigates the role of aneuploidy and chromosome instability (CIN) in the aging brain. Aneuploidy refers to an abnormal chromosomal count, deviating from the normal diploid set. It can manifest as either a deficiency or excess of chromosomes. CIN encompasses a broader range of chromosomal alterations, including aneuploidy as well as structural modifications in DNA. We provide an overview of the state-of-the-art methodologies utilized for studying aneuploidy and CIN in non-tumor somatic tissues devoid of clonally expanded populations of aneuploid cells.CIN and aneuploidy, well-established hallmarks of cancer cells, are also associated with the aging process. In non-transformed cells, aneuploidy can contribute to functional impairment and developmental disorders. Despite the importance of understanding the prevalence and specific consequences of aneuploidy and CIN in the aging brain, these aspects remain incompletely understood, emphasizing the need for further scientific investigations.This comprehensive review consolidates the present understanding, addresses discrepancies in the literature, and provides valuable insights for future research efforts.

Somatic mosaicism in the diseased brain

It is hard to believe that all the cells of a human brain share identical genomes. Indeed, single cell genetic studies have demonstrated intercellular genomic variability in the normal and diseased brain. Moreover, there is a growing amount of evidence on the contribution of somatic mosaicism (the presence of genetically different cell populations in the same individual/tissue) to the etiology of brain diseases. However, brain-specific genomic variations are generally overlooked during the research of genetic defects associated with a brain disease. Accordingly, a review of brain-specific somatic mosaicism in disease context seems to be required. Here, we overview gene mutations, copy number variations and chromosome abnormalities (aneuploidy, deletions, duplications and supernumerary rearranged chromosomes) detected in the neural/neuronal cells of the diseased brain. Additionally, chromosome instability in non-cancerous brain diseases is addressed. Finally, theoretical analysis of possible mechanisms for neurodevelopmental and neurodegenerative disorders indicates that a genetic background for formation of somatic (chromosomal) mosaicism in the brain is likely to exist. In total, somatic mosaicism affecting the central nervous system seems to be a mechanism of brain diseases.

Copy number variations across the blood-brain barrier in multiple sclerosis

Objective

Multiple sclerosis (MS) is a neuroinflammatory disease where immune cells cross the blood–brain barrier (BBB) into the central nervous system (CNS). What predisposes these immune cells to cross the BBB is still unknown. Here, we examine the possibility that genomic rearrangements could predisposespecific immune cells in the peripheral blood to cross the BBB and form sub‐populations of cells involved in the inflammatory process in the CNS.

Methods

We compared copy number variations in paired peripheral blood mononuclear cells (PBMCs) and cerebrospinal fluid (CSF) cells from MS patients. Thereafter, using next generation sequencing, we studied the T‐cell receptor beta (TRB) locus rearrangements and profiled the αβ T cell repertoire in peripheral CD4⁺ and CD8⁺ T cells and in the CSF.

Results

We identified deletions in the T‐cell receptor alpha/delta (TRA/D), gamma (TRG), and TRB loci in CSF cells compared to PBMCs. Further characterization revealed diversity of the TRB locus which was used to describe the character and clonal expansion of T cells in the CNS. T‐cell repertoire profiling from either side of the BBB concluded that the most frequent clones in the CSF samples are unique to an individual. Furthermore, we observed a difference in the proportion of expanded T‐cell clones when comparing samples from MS patients in relapse and remission with opposite trends in CSF and peripheral blood.

Interpretation

This study provides a characterization of the T cells in the CSF and might indicate a role of expanded clones in MS pathogenicity.

Somatic copy number variants in neuropsychiatric disorders

Copy number variants (CNVs) have been implicated in neuropsychiatric disorders, with rare-inherited and de novo CNVs (dnCNVs) having large effects on disease liability. Recent studies started exploring a class of dnCNVs that occur post-zygotically, and are therefore present in some but not all cells of the body. Analogous to conditional mutations in animal models, the presence of risk mutations in a fraction of cells has the potential to enlighten how damaging mutations affect cell-type/cell-circuit specific pathologies leading to neuropsychiatric manifestations. Although mosaic CNVs appear to contribute to a modest fraction of risk (0.3-0.5%), expanding our insights about them with more sensitive experimental and statistical methods, has the potential to help clarify mechanisms of neuropsychiatric disease.

The RNA-Binding Protein HuD Regulates Alternative Splicing and Alternative Polyadenylation in the Mouse Neocortex

The neuronal Hu/ELAV-like proteins HuB, HuC and HuD are a class of RNA-binding proteins that are crucial for proper development and maintenance of the nervous system. These proteins bind to AU-rich elements (AREs) in the untranslated regions (3'-UTRs) of target mRNAs regulating mRNA stability, transport and translation. In addition to these cytoplasmic functions, Hu proteins have been implicated in alternative splicing and alternative polyadenylation in the nucleus. The purpose of this study was to identify transcriptome-wide effects of HuD deletion on both of these nuclear events using RNA sequencing data obtained from the neocortex of Elavl4-/- (HuD KO) mice. HuD KO affected alternative splicing of 310 genes, including 17 validated HuD targets such as Cbx3, Cspp1, Snap25 and Gria2. In addition, deletion of HuD affected polyadenylation of 53 genes, with the majority of significantly altered mRNAs shifting towards usage of proximal polyadenylation signals (PAS), resulting in shorter 3'-UTRs. None of these genes overlapped with those showing alternative splicing events. Overall, HuD KO had a greater effect on alternative splicing than polyadenylation, with many of the affected genes implicated in several neuronal functions and neuropsychiatric disorders.

Genetic factors influencing a neurobiological substrate for psychiatric disorders

A retrospective meta-analysis of magnetic resonance imaging voxel-based morphometry studies proposed that reduced gray matter volumes in the dorsal anterior cingulate and the left and right anterior insular cortex-areas that constitute hub nodes of the salience network-represent a common substrate for major psychiatric disorders. Here, we investigated the hypothesis that the common substrate serves as an intermediate phenotype to detect genetic risk variants relevant for psychiatric disease. To this end, after a data reduction step, we conducted genome-wide association studies of a combined common substrate measure in four population-based cohorts (n = 2271), followed by meta-analysis and replication in a fifth cohort (n = 865). After correction for covariates, the heritability of the common substrate was estimated at 0.50 (standard error 0.18). The top single-nucleotide polymorphism (SNP) rs17076061 was associated with the common substrate at genome-wide significance and replicated, explaining 1.2% of the common substrate variance. This SNP mapped to a locus on chromosome 5q35.2 harboring genes involved in neuronal development and regeneration. In follow-up analyses, rs17076061 was not robustly associated with psychiatric disease, and no overlap was found between the broader genetic architecture of the common substrate and genetic risk for major depressive disorder, bipolar disorder, or schizophrenia. In conclusion, our study identified that common genetic variation indeed influences the common substrate, but that these variants do not directly translate to increased disease risk. Future studies should investigate gene-by-environment interactions and employ functional imaging to understand how salience network structure translates to psychiatric disorder risk.

Genome-Wide DNA Methylation Patterns in Children Exposed to Nonpharmacologically Treated Prenatal Depressive Symptoms: Results From 2 Independent Cohorts

BACKGROUND

Maternal depressive symptoms are a common phenomenon during pregnancy and are related to negative outcomes for child development and health. Modifications in child DNA methylation are discussed as an underlying mechanism for the association between prenatal depressive symptoms and alterations in child outcomes. However, formerly reported genome-wide associations have yet to be replicated.

METHODS

In an epigenome-wide association study (EWAS), alterations of DNA methylation related to maternal prenatal depressive symptoms were investigated in buccal cell samples from 174 children (n = 52 exposed to prenatal depressive symptoms; 6-9 years old) of the German longitudinal study FRAMES-FRANCES. Whole blood samples from the independent, age-comparable ARIES subsample of the ARIES/ALSPAC study (n = 641; n = 159 exposed to prenatal depressive symptoms; 7-8 years old) were examined as a confirmation sample. Depressive symptoms were assessed with the Edinburgh Postnatal Depression Scale. DNA methylation was analyzed with the Infinium Human Methylation 450k BeadChip. Modifications in single CpGs, regions, and biological pathways were investigated. Results were adjusted for age and birth outcomes as well as postnatal and current maternal depressive symptoms. Analyses were performed for the whole sample as well as separated for sex.

RESULTS

The EWAS yielded no differentially methylated CpG or region as well as no accordance between samples withstanding correction for multiple testing. In pathway analyses, no overlapping functional domain was found to be enriched for either sample. A comparison of current and former findings suggests some overlapping methylation modifications from infancy to childhood. Results suggest that there might be sex-specific differential methylation, which should be further investigated in additional studies.

CONCLUSIONS

The current, mainly nonsignificant, results challenge the assumption of consistent modifications of DNA methylation in children exposed to prenatal depressive symptoms. Despite the relatively small sample size used in this study, this lack of significant results may reflect diverse issues of environmental epigenetic studies, which need to be addressed in future research.

Assessment of somatic single-nucleotide variation in brain tissue of cases with schizophrenia

The genetic architecture of schizophrenia (SCZ) includes numerous risk loci across a range of frequencies and sizes, including common and rare single-nucleotide variants and insertions/deletions (indels), as well as rare copy number variants (CNVs). Despite the clear heritability of the disease, monozygotic twins are discordant for SCZ at a significant rate. Somatic variants-genetic changes that arise after fertilization rather than through germline inheritance-are widespread in the human brain and known to contribute to risk for both rare and common neuropsychiatric conditions. The contribution of somatic variants in the brain to risk of SCZ remains to be determined. In this study, we surveyed somatic single-nucleotide variants (sSNVs) in the brains of controls and individuals with SCZ (n = 10 and n = 9, respectively). From each individual, whole-exome sequencing (WES) was performed on DNA from neuronal and non-neuronal nuclei isolated by fluorescence activated nuclear sorting (FANS) from frozen postmortem prefrontal cortex (PFC) samples, as well as DNA extracted from temporal muscle as a reference. We identified an increased burden of sSNVs in cases compared to controls (SCZ rate = 2.78, control rate = 0.70; P = 0.0092, linear mixed effects model), that included a higher rate of non-synonymous and loss-of-function variants (SCZ rate = 1.33, control rate = 0.50; P = 0.047, linear mixed effects model). Our findings suggest sSNVs in the brain may constitute an additional component of the complex genetic architecture of SCZ. This perspective argues for the need to further investigate somatic variation in the brain as an explanation of the discordance in monozygotic twins and a potential guide to the identification of novel therapeutic targets.

Somatic mutations in the human brain: implications for psychiatric research

Psychiatric disorders such as schizophrenia and bipolar disorder are caused by complex gene-environment interactions. While recent advances in genomic technologies have enabled the identification of several risk variants for psychiatric conditions, including single-nucleotide variants and copy-number variations, these factors can explain only a portion of the liability to these disorders. Although non-inherited factors had previously been attributed to environmental causes, recent genomic analyses have demonstrated that de novo mutations are among the main non-inherited risk factors for several psychiatric conditions. Somatic mutations in the brain may also explain how stochastic developmental events and environmental insults confer risk for a psychiatric disorder following fertilization. Here, we review evidence regarding somatic mutations in the brains of individuals with and without neuropsychiatric diseases. We further discuss the potential biological mechanisms underlying somatic mutations in the brain as well as the technical issues associated with the detection of somatic mutations in psychiatric research.

Prioritization and functional assessment of noncoding variants associated with complex diseases

Unraveling functional noncoding variants associated with complex diseases is still a great challenge. We present a novel algorithm, Prioritization And Functional Assessment (PAFA), that prioritizes and assesses the functionality of genetic variants by introducing population differentiation measures and recalibrating training variants. Comprehensive evaluations demonstrate that PAFA exhibits much higher sensitivity and specificity in prioritizing noncoding risk variants than existing methods. PAFA achieves improved performance in distinguishing both common and rare recurrent variants from non-recurrent variants by integrating multiple annotations and metrics. An integrated platform was developed, providing comprehensive functional annotations for noncoding variants by integrating functional genomic data, which can be accessed at http://159.226.67.237:8080/pafa .

Identification of somatic mutations in postmortem human brains by whole genome sequencing and their implications for psychiatric disorders

AIM

Somatic mutations in the human brain are hypothesized to contribute to the functional diversity of brain cells as well as the pathophysiology of neuropsychiatric diseases. However, there are still few reports on somatic mutations in non-neoplastic human brain tissues. This study attempted to unveil the landscape of somatic mutations in the human brain.

METHODS

We explored the landscape of somatic mutations in human brain tissues derived from three individuals with no neuropsychiatric diseases by whole-genome deep sequencing at a depth of around 100. The candidate mutations underwent multi-layered filtering, and were validated by ultra-deep target amplicon sequencing at a depth of around 200 000.

RESULTS

Thirty-one somatic mutations were identified in the human brain, demonstrating the utility of whole-genome sequencing of bulk brain tissue. The mutations were enriched in neuron-expressed genes, and two-thirds of the identified somatic single nucleotide variants in the brain tissues were cytosine-to-thymine transitions, half of which were in CpG dinucleotides.

CONCLUSION

Our developed filtering and validation approaches will be useful to identify somatic mutations in the human brain. The vulnerability of neuron-expressed genes to mutational events suggests their potential relevance to neuropsychiatric diseases.

A novel relationship for schizophrenia, bipolar and major depressive disorder Part 5: a hint from chromosome 5 high density association screen

Familial clustering of schizophrenia (SCZ), bipolar disorder (BPD), and major depressive disorder (MDD) was systematically reported (Aukes, M. F. Genet Med 2012, 14, 338-341) and any two or even three of these disorders could co-exist in some families. In addition, evidence from symptomatology and psychopharmacology also imply that there are intrinsic connections between these three major disorders. A total of 56,569 single nucleotide polymorphism (SNPs) on chromosome 5 were genotyped by Affymetrix Genome-Wide Human SNP array 6.0 on 119 SCZ, 253 BPD (type-I), 177 MDD patients and 1000 controls. Associated SNPs and flanking genes was screen out systematically, and cadherin pathway genes (CDH6, CDH9, CDH10, CDH12, and CDH18) belong to outstanding genes. Unexpectedly, nearly all flanking genes of the associated SNPs distinctive for BPD and MDD were replicated in an enlarged cohort of 986 SCZ patients (P ≤ 9.9E-8). Considering multiple bits of evidence, our chromosome 5 analyses implicated that bipolar and major depressive disorder might be subtypes of schizophrenia rather than two independent disease entities. Also, cadherin pathway genes play important roles in the pathogenesis of the three major mental disorders.

Human brain harbors single nucleotide somatic variations in functionally relevant genes possibly mediated by oxidative stress

Somatic variation in DNA can cause cells to deviate from the preordained genomic path in both disease and healthy conditions. Here, using exome sequencing of paired tissue samples, we show that the normal human brain harbors somatic single base variations measuring up to 0.48% of the total variations. Interestingly, about 64% of these somatic variations in the brain are expected to lead to non-synonymous changes, and as much as 87% of these represent G:C>T:A transversion events. Further, the transversion events in the brain were mostly found in the frontal cortex, whereas the corpus callosum from the same individuals harbors the reference genotype. We found a significantly higher amount of 8-OHdG (oxidative stress marker) in the frontal cortex compared to the corpus callosum of the same subjects (p<0.01), correlating with the higher G:C>T:A transversions in the cortex. We found significant enrichment for axon guidance and related pathways for genes harbouring somatic variations. This could represent either a directed selection of genetic variations in these pathways or increased susceptibility of some loci towards oxidative stress. This study highlights that oxidative stress possibly influence single nucleotide somatic variations in normal human brain.

Gene expression profiling of the dorsolateral and medial orbitofrontal cortex in schizophrenia

Schizophrenia is a complex polygenic disorder of unknown etiology. Over 3,000 candidate genes associated with schizophrenia have been reported, most of which being mentioned only once. Alterations in cognitive processing - working memory, metacognition and mentalization - represent a core feature of schizophrenia, which indicates the involvement of the prefrontal cortex in the pathophysiology of this disorder. Hence we compared the gene expression in postmortem tissue from the left and right dorsolateral prefrontal cortex (DLPFC, Brodmann's area 46), and the medial part of the orbitofrontal cortex (MOFC, Brodmann's area 11/12), in six patients with schizophrenia and six control brains. Although in the past decade several studies performed transcriptome profiling in schizophrenia, this is the first study to investigate both hemispheres, providing new knowledge about possible brain asymmetry at the level of gene expression and its relation to schizophrenia. We found that in the left hemisphere, twelve genes from the DLPFC and eight genes from the MOFC were differentially expressed in patients with schizophrenia compared to controls. In the right hemisphere there was only one gene differentially expressed in the MOFC. We reproduce the involvement of previously reported genes TARDBP and HNRNPC in the pathogenesis of schizophrenia, and report seven novel genes: SART1, KAT7, C1D, NPM1, EVI2A, XGY2, and TTTY15. As the differentially expressed genes only partially overlap with previous studies that analyzed other brain regions, our findings indicate the importance of considering prefrontal cortical regions, especially those in the left hemisphere, for obtaining disease-relevant insights.

Human brain harbors single nucleotide somatic variations in functionally relevant genes possibly mediated by oxidative stress

Somatic variation in DNA can cause cells to deviate from the preordained genomic path in both disease and healthy conditions. Here, using exome sequencing of paired tissue samples, we show that the normal human brain harbors somatic single base variations measuring up to 0.48% of the total variations. Interestingly, about 64% of these somatic variations in the brain are expected to lead to non-synonymous changes, and as much as 87% of these represent G:C>T:A transversion events. Further, the transversion events in the brain were mostly found in the frontal cortex, whereas the corpus callosum from the same individuals harbors the reference genotype. We found a significantly higher amount of 8-OHdG (oxidative stress marker) in the frontal cortex compared to the corpus callosum of the same subjects (p<0.01), correlating with the higher G:C>T:A transversions in the cortex. We found significant enrichment for axon guidance and related pathways for genes harbouring somatic variations. This could represent either a directed selection of genetic variations in these pathways or increased susceptibility of some loci towards oxidative stress. This study highlights that oxidative stress possibly influence single nucleotide somatic variations in normal human brain.

Somatic mutations in disorders with disrupted brain connectivity

Mutations occur during cell division in all somatic lineages. Because neurogenesis persists throughout human life, somatic mutations in the brain arise during development and accumulate with the aging process. The human brain consists of 100 billion neurons that form an extraordinarily intricate network of connections to achieve higher level cognitive functions. Due to this network architecture, perturbed neuronal functions are rarely restricted to a focal area; instead, they are often spread via the neuronal network to affect other connected areas. Although somatic diversity is an evident feature of the brain, the extent to which somatic mutations affect the neuronal structure and function and their contribution to neurological disorders associated with disrupted brain connectivity remain largely unexplored. Notably, recent reports indicate that brain somatic mutations can indeed play a critical role that leads to the structural and functional abnormalities of the brain observed in several neurodevelopmental disorders. Here, I review the extent and significance of brain somatic mutations and provide my perspective regarding these mutations as potential molecular lesions underlying relatively common conditions with disrupted brain connectivity. Moreover, I discuss emerging technical platforms that will facilitate the detection of low-frequency somatic mutations and validate the biological functions of the identified mutations in the context of brain connectivity.

BOD1 Is Required for Cognitive Function in Humans and Drosophila

Here we report a stop-mutation in the BOD1 (Biorientation Defective 1) gene, which co-segregates with intellectual disability in a large consanguineous family, where individuals that are homozygous for the mutation have no detectable BOD1 mRNA or protein. The BOD1 protein is required for proper chromosome segregation, regulating phosphorylation of PLK1 substrates by modulating Protein Phosphatase 2A (PP2A) activity during mitosis. We report that fibroblast cell lines derived from homozygous BOD1 mutation carriers show aberrant localisation of the cell cycle kinase PLK1 and its phosphatase PP2A at mitotic kinetochores. However, in contrast to the mitotic arrest observed in BOD1-siRNA treated HeLa cells, patient-derived cells progressed through mitosis with no apparent segregation defects but at an accelerated rate compared to controls. The relatively normal cell cycle progression observed in cultured cells is in line with the absence of gross structural brain abnormalities in the affected individuals. Moreover, we found that in normal adult brain tissues BOD1 expression is maintained at considerable levels, in contrast to PLK1 expression, and provide evidence for synaptic localization of Bod1 in murine neurons. These observations suggest that BOD1 plays a cell cycle-independent role in the nervous system. To address this possibility, we established two Drosophila models, where neuron-specific knockdown of BOD1 caused pronounced learning deficits and significant abnormalities in synapse morphology. Together our results reveal novel postmitotic functions of BOD1 as well as pathogenic mechanisms that strongly support a causative role of BOD1 deficiency in the aetiology of intellectual disability. Moreover, by demonstrating its requirement for cognitive function in humans and Drosophila we provide evidence for a conserved role of BOD1 in the development and maintenance of cognitive features.

Gene expression profiling of the dorsolateral and medial orbitofrontal cortex in schizophrenia

Schizophrenia is a complex polygenic disorder of unknown etiology. Over 3,000 candidate genes associated with schizophrenia have been reported, most of which being mentioned only once. Alterations in cognitive processing - working memory, metacognition and mentalization - represent a core feature of schizophrenia, which indicates the involvement of the prefrontal cortex in the pathophysiology of this disorder. Hence we compared the gene expression in postmortem tissue from the left and right dorsolateral prefrontal cortex (DLPFC, Brodmann's area 46), and the medial part of the orbitofrontal cortex (MOFC, Brodmann's area 11/12), in six patients with schizophrenia and six control brains. Although in the past decade several studies performed transcriptome profiling in schizophrenia, this is the first study to investigate both hemispheres, providing new knowledge about possible brain asymmetry at the level of gene expression and its relation to schizophrenia. We found that in the left hemisphere, twelve genes from the DLPFC and eight genes from the MOFC were differentially expressed in patients with schizophrenia compared to controls. In the right hemisphere there was only one gene differentially expressed in the MOFC. We reproduce the involvement of previously reported genes TARDBP and HNRNPC in the pathogenesis of schizophrenia, and report seven novel genes: SART1, KAT7, C1D, NPM1, EVI2A, XGY2, and TTTY15. As the differentially expressed genes only partially overlap with previous studies that analyzed other brain regions, our findings indicate the importance of considering prefrontal cortical regions, especially those in the left hemisphere, for obtaining disease-relevant insights.

SoloDel: a probabilistic model for detecting low-frequent somatic deletions from unmatched sequencing data

MOTIVATION

Finding somatic mutations from massively parallel sequencing data is becoming a standard process in genome-based biomedical studies. There are a number of robust methods developed for detecting somatic single nucleotide variations However, detection of somatic copy number alteration has been substantially less explored and remains vulnerable to frequently raised sampling issues: low frequency in cell population and absence of the matched control samples.

RESULTS

We developed a novel computational method SoloDel that accurately classifies low-frequent somatic deletions from germline ones with or without matched control samples. We first constructed a probabilistic, somatic mutation progression model that describes the occurrence and propagation of the event in the cellular lineage of the sample. We then built a Gaussian mixture model to represent the mixed population of somatic and germline deletions. Parameters of the mixture model could be estimated using the expectation-maximization algorithm with the observed distribution of read-depth ratios at the points of discordant-read based initial deletion calls. Combined with conventional structural variation caller, SoloDel greatly increased the accuracy in classifying somatic mutations. Even without control, SoloDel maintained a comparable performance in a wide range of mutated subpopulation size (10-70%). SoloDel could also successfully recall experimentally validated somatic deletions from previously reported neuropsychiatric whole-genome sequencing data.

AVAILABILITY AND IMPLEMENTATION

Java-based implementation of the method is available at http://sourceforge.net/projects/solodel/

CONTACT

swkim@yuhs.ac or dhlee@biosoft.kaist.ac.kr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.