A rare missense variant in RCL1 segregates with depression in extended families

A rare variant in GPR156 associated with depression in a Mennonite pedigree causes habenula hyperactivity and stress sensitivity in mice

Major depressive disorder (MDD) is a leading cause of disability worldwide. Risk for MDD is heritable, and the genetic structure of founder populations enables investigation of rare susceptibility alleles with large effect. In an extended Old Order Mennonite family cohort, we identified a rare missense variant in GPR156 (c.1599G>T, p.Glu533Asp) associated with a two-fold increase in the relative risk of MDD. GPR156 is an orphan G protein-coupled receptor localized in the medial habenula, a region implicated in mood regulation. Insertion of a human sequence containing c.1599G>T into the murine Gpr156 locus induced medial habenula hyperactivity and abnormal stress-related behaviors. This work reveals a human variant that is associated with depression, implicates GPR156 as a target for mood regulation, and introduces informative murine models for investigating the pathophysiology and treatment of affective disorders.

Whole-Genome Sequencing Reveals Individual and Cohort Level Insights into Chromosome 9p Syndromes

Previous genomic efforts on chromosome 9p deletion and duplication syndromes have utilized low resolution strategies (i.e., karyotypes, chromosome microarrays). We present the first large-scale whole-genome sequencing (WGS) study of 100 individuals from families with 9p-related syndromes including 85 unrelated probands through the 9P-ARCH (Advanced Research in Chromosomal Health: Genomic, Phenotypic, and Functional Aspects of 9p-Related syndromes) research network. We analyzed the genomic architecture of these syndromes, highlighting fundamental features and their commonalities and differences across individuals. This work includes a machine-learning model that predicts 9p deletion syndrome from gene copy number estimates using WGS data. Two Late Replicating Regions (LRR1 [a previously un-named human fragile site], LRR2) were identified that contain most structural variant breakpoints in 9p deletion syndrome pointing to replication-based issues in structural variant formation. Furthermore, we show the utility of using WGS information to obtain a comprehensive understanding of 9p-related variation in an individual with complex structural variation where chromothripsis is the likely mechanism. Genes on 9p were prioritized based on statistical assessment of human genomic variation. Furthermore, through application of spatial transcriptomics to embryonic mouse tissue we examined 9p-gene expression in craniofacial and brain development. Through these strategies, we identified 24 important genes for the majority (83%) of individuals with 9p deletion syndrome including AK3, BRD10, CD274, CDC37L1, DMRT1, DMRT2, DMRT3, DOCK8, GLIS3, JAK2, KANK1, KDM4C, PLPP6, PTPRD, PUM3, RANBP6, RCL1, RFX3, RIC1, SLC1A1, SMARCA2, UHRF2, VLDLR, and ZNG1A. Two genes (AK3, ZNG1A) are involved in mitochondrial function and testing of the mitochondrial genome revealed excess copy number in individuals with 9p deletion syndrome. This study presents the most comprehensive genomic analysis of 9p-related syndromes to date, with plans for further expansion through our 9P-ARCH research network.

Human mood disorder risk gene Synaptotagmin-14 contributes to mania-like behaviors in mice

Bipolar disorder (BD) and major depressive disorder (MDD) are the most prevalent mood disorders and cause considerable burden worldwide. Compelling evidence suggests a pronounced overlap between these two disorders in clinical symptoms, treatment strategies, and genetic etiology. Here we leverage a BD GWAS (1822 cases and 4650 controls) and a MDD GWAS (5303 cases and 5337 controls), followed by independent replications, to investigate their shared genetic basis among Han Chinese. We have herein identified a lead SNP rs126277 at the 1q32.2 locus, which also exhibited nominal associations with mood disorders and several relevant sub-clinical phenotypes (e.g., mania) in European populations. Bulk tissue and single-cell eQTL analyses suggest that the risk G-allele of rs126277 predicted lower SYT14 mRNA expression in human brains. We generated mice lacking Syt14 (Syt14-/-) and mice with insufficient expression of Syt14 in the hippocampus (Syt14-KD), and found that depletion of Syt14 resulted in mania-like behaviors including hyperactivity and anti-depressive behaviors, resembling aspects of mood disorders. We also confirmed that deficiency of this gene in the hippocampus was sufficient to induce hyperactivity in mice. RNA-sequencing analyses of the hippocampus of Syt14-/- mice revealed significant upregulation of Per1 as well as downregulation of Slc7a11 and Ptprb. Ultrastructural analyses showed significant alteration of the number of vesicles within 50 nm to the active zone and the width of synaptic cleft in the ventral hippocampus of Syt14-/- mice compared with the control mice. Overall, we have identified a novel mood disorder risk gene SYT14, and confirmed its impact on mania-like behaviors. While the current study identifies an essential mood disorder risk gene, further investigations elucidating the detailed mechanisms by which SYT14 contributes to the pathogenesis of the illnesses are needed.

Interlaminar and varicose-projection astrocytes: toward a new understanding of the primate brain

In the last years, science started to move toward a more glio-neurocentric view, in which astrocytes are hypothesized to be directly involved in cognitive functions. Indeed, astrocytes show a variety of shapes with species-specific characteristics, suggesting a specialization of roles during evolution. Interlaminar (ILA) and varicose-projection (VP-As) astrocytes show an anatomical organization that is different compared to the classical horizontal net typically formed by protoplasmic and fibrous astrocytes. ILAs show a modular architecture with the soma in the first cortical layer and processes toward the deep layers with species-specific length. VP-As reside in the deep layers of the cortex, are characterized by varicosities on the longest processes, and are individual-specific. These characteristics suggest roles that are more complex than what was theorized until now. Here, we recapitulate what we know so far from literature from the first time ILAs were described to the most recent discoveries, spanning from morphology description, hypothesis on the development to their features in diseases. For a complete glance on this topic, we included a final paragraph on which techniques and models were used to study ILAs and VP-As, and what new avenues may be opened thanks to more novel methods.

Human Pluripotent Stem Cell-Derived Astrocyte Functionality Compares Favorably with Primary Rat Astrocytes

Astrocytes are essential for the formation and maintenance of neural networks. However, a major technical challenge for investigating astrocyte function and disease-related pathophysiology has been the limited ability to obtain functional human astrocytes. Despite recent advances in human pluripotent stem cell (hPSC) techniques, primary rodent astrocytes remain the gold standard in coculture with human neurons. We demonstrate that a combination of leukemia inhibitory factor (LIF) and bone morphogenetic protein-4 (BMP4) directs hPSC-derived neural precursor cells to a highly pure population of astroglia in 28 d. Using single-cell RNA sequencing, we confirm the astroglial identity of these cells and highlight profound transcriptional adaptations in cocultured hPSC-derived astrocytes and neurons, consistent with their further maturation. In coculture with human neurons, multielectrode array recordings revealed robust network activity of human neurons in a coculture with hPSC-derived or rat astrocytes [3.63 ± 0.44 min-1 (hPSC-derived), 2.86 ± 0.64 min-1 (rat); p = 0.19]. In comparison, we found increased spike frequency within network bursts of human neurons cocultured with hPSC-derived astrocytes [56.31 ± 8.56 Hz (hPSC-derived), 24.77 ± 4.04 Hz (rat); p < 0.01], and whole-cell patch-clamp recordings revealed an increase of postsynaptic currents [2.76 ± 0.39 Hz (hPSC-derived), 1.07 ± 0.14 Hz (rat); p < 0.001], consistent with a corresponding increase in synapse density [14.90 ± 1.27/100 μm2 (hPSC-derived), 8.39 ± 0.63/100 μm2 (rat); p < 0.001]. Taken together, we show that hPSC-derived astrocytes compare favorably with rat astrocytes in supporting human neural network activity and maturation, providing a fully human platform for investigating astrocyte function and neuronal-glial interactions.

Integration of whole-exome sequencing and structural neuroimaging analysis in major depressive disorder: a joint study

Major depressive disorder (MDD) is a common mental illness worldwide and is triggered by an intricate interplay between environmental and genetic factors. Although there are several studies on common variants in MDD, studies on rare variants are relatively limited. In addition, few studies have examined the genetic contributions to neurostructural alterations in MDD using whole-exome sequencing (WES). We performed WES in 367 patients with MDD and 161 healthy controls (HCs) to detect germline and copy number variations in the Korean population. Gene-based rare variants were analyzed to investigate the association between the genes and individuals, followed by neuroimaging-genetic analysis to explore the neural mechanisms underlying the genetic impact in 234 patients with MDD and 135 HCs using diffusion tensor imaging data. We identified 40 MDD-related genes and observed 95 recurrent regions of copy number variations. We also discovered a novel gene, FRMPD3, carrying rare variants that influence MDD. In addition, the single nucleotide polymorphism rs771995197 in the MUC6 gene was significantly associated with the integrity of widespread white matter tracts. Moreover, we identified 918 rare exonic missense variants in genes associated with MDD susceptibility. We postulate that rare variants of FRMPD3 may contribute significantly to MDD, with a mild penetration effect.

Correlation between immune-related genes and depression-like features in an animal model and in humans

A growing body of evidence suggests that immune-related genes play pivotal roles in the pathophysiology of depression. In the present study, we investigated a plausible connection between gene expression, DNA methylation, and brain structural changes in the pathophysiology of depression using a combined approach of murine and human studies. We ranked the immobility behaviors of 30 outbred Crl:CD1 (ICR) mice in the forced swim test (FST) and harvested their prefrontal cortices for RNA sequencing. Of the 24,532 analyzed genes, 141 showed significant correlations with FST immobility time, as determined through linear regression analysis with p ≤ 0.01. The identified genes were mostly involved in immune responses, especially interferon signaling pathways. Moreover, induction of virus-like neuroinflammation in the brains of two separate mouse cohorts (n = 30 each) using intracerebroventricular polyinosinic:polycytidylic acid injection resulted in increased immobility during FST and similar expression of top immobility-correlated genes. In human blood samples, candidate gene (top 5%) expression profiling using DNA methylation analysis found the interferon-related USP18 (cg25484698, p = 7.04 × 10-11, Δβ = 1.57 × 10-2; cg02518889, p = 2.92 × 10-3, Δβ =  - 8.20 × 10-3) and IFI44 (cg07107453, p = 3.76 × 10-3, Δβ =  - 4.94 × 10-3) genes to be differentially methylated between patients with major depressive disorder (n = 350) and healthy controls (n = 161). Furthermore, cortical thickness analyses using T1-weighted images revealed that the DNA methylation scores for USP18 were negatively correlated with the thicknesses of several cortical regions, including the prefrontal cortex. Our results reveal the important role of the interferon pathway in depression and suggest USP18 as a potential candidate target. The results of the correlation analysis between transcriptomic data and animal behavior carried out in this study provide insights that could enhance our understanding of depression in humans.

Genome-wide association studies on Northern Italy isolated populations provide further support concerning genetic susceptibility for major depressive disorder

OBJECTIVES

Major depressive disorder (MDD) is a psychiatric disorder with pathogenesis influenced by both genetic and environmental factors. To date, the molecular-level understanding of its aetiology remains unclear. Thus, we aimed to identify genetic variants and susceptibility genes for MDD with a genome-wide association study (GWAS) approach.

METHODS

We performed a meta-analysis of GWASs and a gene-based analysis on two Northern Italy isolated populations (cases/controls n = 166/472 and 33/320), followed by replication and polygenic risk score (PRS) analyses in Italian independent samples (cases n = 464, controls n = 339).

RESULTS

We identified two novel MDD-associated genes, KCNQ5 (lead SNP rs867262, p = 3.82 × 10^-9) and CTNNA2 (rs6729523, p = 1.25 × 10^-8). The gene-based analysis revealed another six genes (p < 2.703 × 10^-6): GRM7, CTNT4, SNRK, SRGAP3, TRAPPC9, and FHIT. No replication of the genome-wide significant SNPs was found in the independent cohort, even if 14 SNPs around CTNNA2 showed association with MDD and related phenotypes at the nominal level of p (<0.05). Furthermore, the PRS model developed in the discovery cohort discriminated cases and controls in the replication cohort.

CONCLUSIONS

Our work suggests new possible genes associated with MDD, and the PRS analysis confirms the polygenic nature of this disorder. Future studies are required to better understand the role of these findings in MDD.

Modification by genetic polymorphism of lead-induced IQ alteration: a systematic review

As well as a lead-related environmental factor, genetic factors could also corroborate important changes in intelligence quotient (IQ) through single-nucleotide polymorphisms. Thus, a systematic review was carried out to evaluate the possible influence of polymorphism on blood Pb levels and IQ points in pediatric patients (0-19 years old). Following the PRISMA guideline, the studies were systematically collected on PubMed, Scopus, and Embase databases. Six genes (transferrin (TF); glutamate ionotropic receptor NMDA type subunit 2A (GRIN2A); glutamate ionotropic receptor NMDA type subunit 2B (GRIN2B); dopamine receptor D2/ankyrin repeat and kinase domain containing 1 ankyrin repeat and kinase domain containing 1 (DRD2/ANKK1); aminolevulinate dehydratase (ALAD); vitamin D receptor (VDR)) were found in six selected articles. In these genes, 11 single-nucleotide polymorphisms were searched and six different types of variations (missense variant, intron variant, synonymous variant, stop, stop gained) were observed. Due to the few studies in the literature, there is no conclusive data to point out that there is a direct relationship between polymorphisms, Pb levels, and reduction of IQ points.

Novel CAPN1 missense variants in complex hereditary spastic paraplegia with early-onset psychosis

CAPN1-associated hereditary spastic paraplegia (SPG76) is a rare and clinically heterogenous syndrome due to loss of calpain-1 function. Here we illustrate a translational approach to the case of an 18-year-old patient who first presented with psychiatric symptoms followed by spastic gait, intention tremor, and neurogenic bladder dysfunction, consistent with a complex form of HSP. Exome sequencing showed compound-heterozygous missense variants in CAPN1 (NM_001198868.2: c.1712A>G (p.Asn571Ser)/c.1991C>T (p.Ser664Leu)) and a previously reported heterozygous stop-gain variant in RCL1. In silico analyses of the CAPN1 variants predicted a deleterious effect and in vitro functional studies confirmed reduced calpain-1 activity and dysregulated downstream signaling. These findings support a diagnosis of SPG76 and highlight that the psychiatric symptoms can precede the motor symptoms in HSP. Our results also suggest that multiple genes can potentially contribute to complex neuropsychiatric diseases.

The genetic basis of major depression

Major depressive disorder (MDD) is a common, debilitating, phenotypically heterogeneous disorder with heritability ranges from 30% to 50%. Compared to other psychiatric disorders, its high prevalence, moderate heritability, and strong polygenicity have posed major challenges for gene-mapping in MDD. Studies of common genetic variation in MDD, driven by large international collaborations such as the Psychiatric Genomics Consortium, have confirmed the highly polygenic nature of the disorder and implicated over 100 genetic risk loci to date. Rare copy number variants associated with MDD risk were also recently identified. The goal of this review is to present a broad picture of our current understanding of the epidemiology, genetic epidemiology, molecular genetics, and gene-environment interplay in MDD. Insights into the impact of genetic factors on the aetiology of this complex disorder hold great promise for improving clinical care.

Analysis of 200 000 exome-sequenced UK Biobank subjects fails to identify genes influencing probability of developing a mood disorder resulting in psychiatric referral

BACKGROUND

Depression is moderately heritable but there is no common genetic variant which has a major effect on susceptibility. A previous analysis of 50 000 exome-sequenced subjects failed to implicate any genes or sets of genes in which rare variants were associated with risk of affective disorder requiring specialist treatment. A much larger exome-sequenced dataset is now available.

METHODS

Data from 200 632 exome-sequenced UK Biobank participants was analysed. Subjects were treated as cases if they had reported having seen a psychiatrist for 'nerves, anxiety, tension or depression'. Gene-wise weighted burden analysis was performed to see if there were any genes or sets of genes for which there was an excess of rare, functional variants in cases.

RESULTS

There were 22 886 cases and 176 486 controls. There were 22 642 informative genes but no gene or gene set produced a statistically significant result after correction for multiple testing. None of the genes or gene sets with the lowest P values appeared to be an obvious biological candidate.

CONCLUSIONS

The results conform exactly with the expectation under the null hypothesis. It seems unlikely that the use of common, poorly defined phenotypes will produce useful advances in understanding genetic contributions to affective disorder and it might be preferable to focus instead on obtaining large exome-sequenced samples of conditions such as bipolar 1 disorder and severe, recurrent depression. This research has been conducted using the UK Biobank Resource.

Identifying Prognostic Significance of RCL1 and Four-Gene Signature as Novel Potential Biomarkers in HCC Patients

BACKGROUND

Hepatocellular carcinoma (HCC) is a highly malignant disease, and it is characterized by rapid progression and low five-year survival rate. At present, there are no effective methods for monitoring the treatment and prognosis of HCC.

METHODS

The transcriptome and gene expression profiles of HCC were obtained from the Cancer Genome Atlas (TCGA) program, International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) databases. The random forest method was applied to construct a four-gene prognostic model based on RNA terminal phosphate cyclase like 1 (RCL1) expression. The Kaplan-Meier method was performed to evaluate the prognostic value of RCL1, long noncoding RNAs (AC079061, AL354872, and LINC01093), and four-gene signature (SPP1, MYBL2, TRNP1, and FTCD). We examined the relationship between RCL1 expression and immune cells infiltration, tumor mutation burden (TMB), and microsatellite instability (MSI).

RESULTS

The results of multiple databases indicated that the aberrant expression of RCL1 was associated with clinical outcome, immune cells infiltration, TMB, and MSI in HCC patients. Meanwhile, we found that long noncoding RNAs (AC079061, AL354872, and LINC01093) and RCL1 were significantly coexpressed in HCC patients. We also confirmed that the four-gene signature was an independent prognostic factor for HCC patients. Ferroptosis potential index, immune checkpoint molecules, and clinical feature were found to have obvious correlations with risk score. The area under the receiver operating characteristic curve values for the model were 0.7-0.8 in the training set and the validation set, suggesting high robustness of the four-gene signature. We then built a nomogram for facilitating the use in clinical practice.

CONCLUSION

Our study demonstrated that RCL1 and a novel four-gene signature can be used as prognostic biomarkers for predicting clinical outcome in HCC patients; and this model may assist in individualized treatment monitoring of HCC patients in clinical practice.

Rcl1 depletion impairs 18S pre-rRNA processing at the A1-site and up-regulates a cohort of ribosome biogenesis genes in zebrafish

Yeast Rcl1 is a potential endonuclease that mediates pre-RNA cleavage at the A₂-site to separate 18S rRNA from 5.8S and 25S rRNAs. However, the biological function of Rcl1 in opisthokonta is poorly defined. Moreover, there is no information regarding the exact positions of 18S pre-rRNA processing in zebrafish. Here, we report that zebrafish pre-rRNA harbours three major cleavage sites in the 5′ETS, namely –477nt (A′-site), –97nt (A₀-site) and the 5′ETS and 18S rRNA link (A₁-site), as well as two major cleavage regions within the ITS1, namely 208–218nt (site 2) and 20–33nt (site E). We also demonstrate that depletion of zebrafish Rcl1 mainly impairs cleavage at the A₁-site. Phenotypically, rcl1^–/– mutants exhibit a small liver and exocrine pancreas and die before 15 days post-fertilization. RNA-seq analysis revealed that the most significant event in rcl1^–/– mutants is the up-regulated expression of a cohort of genes related to ribosome biogenesis and tRNA production. Our data demonstrate that Rcl1 is essential for 18S rRNA maturation at the A₁-site and for digestive organogenesis in zebrafish. Rcl1 deficiency, similar to deficiencies in other ribosome biogenesis factors, might trigger a common mechanism to upregulate the expression of genes responsible for ribosome biogenesis.

RCL1 copy number variants are associated with a range of neuropsychiatric phenotypes

Mendelian and early-onset severe psychiatric phenotypes often involve genetic variants having a large effect, offering opportunities for genetic discoveries and early therapeutic interventions. Here, the index case is an 18-year-old boy, who at 14 years of age had a decline in cognitive functioning over the course of a year and subsequently presented with catatonia, auditory and visual hallucinations, paranoia, aggression, mood dysregulation, and disorganized thoughts. Exome sequencing revealed a stop-gain mutation in RCL1 (NM_005772.4:c.370 C > T, p.Gln124Ter), encoding an RNA 3'-terminal phosphate cyclase-like protein that is highly conserved across eukaryotic species. Subsequent investigations across two academic medical centers identified eleven additional cases of RCL1 copy number variations (CNVs) with varying neurodevelopmental or psychiatric phenotypes. These findings suggest that dosage variation of RCL1 contributes to a range of neurological and clinical phenotypes.

A linkage and exome study of multiplex families with bipolar disorder implicates rare coding variants of ANK3 and additional rare alleles at 10q11-q21

BACKGROUND

Bipolar disorder is a highly heritable psychiatric condition for which specific genetic factors remain largely unknown. In the present study, we used combined whole-exome sequencing and linkage analysis to identify risk loci and dissect the contribution of common and rare variants in families with a high density of illness.

METHODS

Overall, 117 participants from 15 Australian extended families with bipolar disorder (72 with affective disorder, including 50 with bipolar disorder type I or II, 13 with schizoaffective disorder-manic type and 9 with recurrent unipolar disorder) underwent whole-exome sequencing. We performed genome-wide linkage analysis using MERLIN and conditional linkage analysis using LAMP. We assessed the contribution of potentially functional rare variants using a genebased segregation test.

RESULTS

We identified a significant linkage peak on chromosome 10q11-q21 (maximal single nucleotide polymorphism = rs10761725; exponential logarithm of the odds [LODexp] = 3.03; empirical p = 0.046). The linkage interval spanned 36 protein-coding genes, including a gene associated with bipolar disorder, ankyrin 3 (ANK3). Conditional linkage analysis showed that common ANK3 risk variants previously identified in genome-wide association studies - or variants in linkage disequilibrium with those variants - did not explain the linkage signal (rs10994397 LOD = 0.63; rs9804190 LOD = 0.04). A family-based segregation test with 34 rare variants from 14 genes under the linkage interval suggested rare variant contributions of 3 brain-expressed genes: NRBF2 (p = 0.005), PCDH15 (p = 0.002) and ANK3 (p = 0.014).

LIMITATIONS

We did not examine non-coding variants, but they may explain the remaining linkage signal.

CONCLUSION

Combining family-based linkage analysis with next-generation sequencing data is effective for identifying putative disease genes and specific risk variants in complex disorders. We identified rare missense variants in ANK3, PCDH15 and NRBF2 that could confer disease risk, providing valuable targets for functional characterization.

Targeted sequencing of the BDNF gene in young Chinese Han people with major depressive disorder

BACKGROUND

Adolescence and young adulthood are considered the peak age for the emergence of many psychiatric disorders, in particular major depressive disorder (MDD). Previous research has shown substantial heritability for MDD. In addition, the brain-derived neurotrophic factor (BDNF) gene is known to be associated with MDD. However, there has been no study conducting targeted sequencing of the BDNF gene in young MDD patients so far.

METHOD

To examine whether the BDNF gene is associated with the occurrence of MDD in young patients, we used targeted sequencing to detect the BDNF gene variants in 259 young Chinese Han people (105 MDD patients and 154 healthy subjects).

RESULTS

The BDNF variant rs4030470 was associated with MDD in young Chinese Han people (uncorrected p = 0.046), but this was no longer significant after applying FDR correction (p = 0.552, after FDR correction). We did not find any significant differences in genotype or haplotype frequencies between the case and control groups, and furthermore discovered no rare mutation variants any of the 259 subjects.

CONCLUSION

Our results do not support an association of the BDNF gene variants with MDD in young people in the Chinese Han population.

Looking into the genetic bases of OCD dimensions: a pilot genome-wide association study

The multidimensional nature of obsessive-compulsive disorder (OCD) has been consistently reported. Clinical and biological characteristics have been associated with OCD dimensions in different ways. Studies suggest the existence of specific genetic bases for the different OCD dimensions. In this study, we analyze the genomic markers, genes, gene ontology and biological pathways associated with the presence of aggressive/checking, symmetry/order, contamination/cleaning, hoarding, and sexual/religious symptoms, as assessed via the Dimensional Yale-Brown Obsessive Compulsive Scale (DY-BOCS) in 399 probands. Logistic regression analyses were performed at the single-nucleotide polymorphism (SNP) level. Gene-based and enrichment analyses were carried out for common (SNPs) and rare variants. No SNP was associated with any dimension at a genome-wide level (p < 5 × 10^-8). Gene-based analyses showed one gene to be associated with hoarding (SETD3, p = 1.89 × 10^-08); a gene highly expressed in the brain and which plays a role in apoptotic processes and transcriptomic changes, and another gene associated with aggressive symptoms (CPE; p = 4.42 × 10^-6), which is involved in neurotrophic functions and the synthesis of peptide hormones and neurotransmitters. Different pathways or biological processes were represented by genes associated with aggressive (zinc ion response and lipid metabolism), order (lipid metabolism), sexual/religious (G protein-mediated processes) and hoarding (metabolic processes and anion transport) symptoms after FDR correction; while no pathway was associated with contamination. Specific genomic bases were found for each dimension assessed, especially in the enrichment analyses. Further research with larger samples and different techniques, such as next-generation sequencing, are needed to better understand the differential genetics of OCD dimensions.

Objectives, design and main findings until 2020 from the Rotterdam Study

The Rotterdam Study is an ongoing prospective cohort study that started in 1990 in the city of Rotterdam, The Netherlands. The study aims to unravel etiology, preclinical course, natural history and potential targets for intervention for chronic diseases in mid-life and late-life. The study focuses on cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, otolaryngological, locomotor, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. Since 2016, the cohort is being expanded by persons aged 40 years and over. The findings of the Rotterdam Study have been presented in over 1700 research articles and reports. This article provides an update on the rationale and design of the study. It also presents a summary of the major findings from the preceding 3 years and outlines developments for the coming period.

Exploring genetic variants in obsessive compulsive disorder severity: A GWAS approach

BACKGROUND

The severity of Obsessive-Compulsive Disorder (OCD) varies significantly among probands. No study has specifically investigated the genetic base of OCD severity. A previous study from our group found an OCD polygenic risk score to predict pre- and post-treatment severity. This study explores the genomic bases of OCD severity.

METHODS

We administered the Yale-Brown Obsessive Compulsive Scale (Y-BOCS) to 401 patients at their first visit to our clinic to measure their OCD severity. Genotyping data was collected by using the Infinium PsychArray-24 BeadChip kit (Illumina). We analyzed genetic association with OCD severity in a linear regression analysis at single-nucleotide polymorphism (SNP)- and gene-levels, this last also considering rare variants. Enrichment analyses were performed from gene-based analyses' results.

RESULTS

No SNP reached significant association (p < 10^-8) with the YBOCS. Six markers showed suggestive association (p < 10^-5). The top SNP was an intergenic variant in chromosome 2: rs7578149 (p < 1.89 × 10^-6), located in a region suggestively associated with MDD. Linkage disequilibrium was found for two clusters of SNPs located between SLC16A14 and SP110 in chromosome 2, all of them forming one peak of association. Enrichment analyses revealed OCD genes to be associated with porin activity (FDR = 0.01) and transmembrane structure (FDR = 0.04).

LIMITATIONS

The size of the sample and the transversal nature of the severity measure are limitations of this study.

CONCLUSION

This study contributes to better characterize OCD at an individual level, helping to know more about the prognosis of the disorder and develop more individualized treatments.

A study combining whole-exome sequencing and structural neuroimaging analysis for major depressive disorder

BACKGROUND

Genetic variations associated with major depressive disorder (MDD) may affect the structural aspects of neural networks mediated by the molecular pathways involved in neuronal survival and synaptic plasticity. However, few studies have applied a novel approach such as whole-exome sequencing (WES) analysis to investigate the genetic contribution to the neurostructural changes in MDD.

METHODS

In the first part of the study, we investigated rare variants of selected genes from previous WES studies using a WES analysis in 184 patients with MDD and 82 healthy controls. In the second part of the study, we explored the association between the common genetic variants from the WES analysis and cortical thickness in 91 patients with MDD and 75 healthy controls. The gray-matter thickness of each cortical region was measured using FreeSurfer.

RESULTS

We identified recurrent non-silent variants in 24 MDD-related genes including FASN, MYH13, UNC13D, LILRA1, CACNA1B, TRIO, HOMER3, and BCAR3, and observed eleven recurrently altered copy number alternations where a gain on 15q11.2 and losses on 7q34 and 15q11.1-q11.2 in MDD genomes. We also found that rs11592462 in CDH23, a calcium-dependent cell-adhesion molecule encoding gene, was significantly associated with thinning in the right anterior cingulate cortex.

LIMITATION

The small sample size may lead our findings to be underpowered regarding rare variants.

CONCLUSION

The present study identified that non-synonymous rare variants were significantly associated with risk of MDD and found that genetic contributions to the development of MDD may be mediated by alterations in cortical thickness of emotion-processing neural circuits.

Genetic Advance in Depressive Disorder

Major depressive disorder (MDD) and bipolar disorder (BPD) are both chronic, severe mood disorder with high misdiagnosis rate, leading to substantial health and economic burdens to patients around the world. There is a high misdiagnosis rate of bipolar depression (BD) just based on symptomology in depressed patients whose previous manic or mixed episodes have not been well recognized. Therefore, it is important for psychiatrists to identify these two major psychiatric disorders. Recently, with the accumulation of clinical sample sizes and the advances of methodology and technology, certain progress in the genetics of major depression and bipolar disorder has been made. This article reviews the candidate genes for MDD and BD, genetic variation loci, chromosome structural variation, new technologies, and new methods.

Rare variants in SLC6A4 cause susceptibility to major depressive disorder with suicidal ideation in Han Chinese adolescents and young adults

BACKGROUND

Suicidal ideation (SI) is the most serious symptom of major depressive disorder (MDD) and considered an extreme state. The serotonin transporter gene (SLC6A4) plays a significant role in MDD and suicide pathophysiology. Previous studies have revealed an association between common variants of SLC6A4 with the risk of MDD and suicide. However, very few studies have so far focused on the degree to which rare variants of SLC6A4 are responsible for the depression observed in adolescent and young adult suicide patients. The aim of this study was to examine the impact of common and rare variants of SLC6A4 on the risk of Han Chinese adolescents and young adults suffering MDD with SI.

METHODS

Targeted sequencing of the SLC6A4 gene was conducted using FastTarget technology in Han Chinese adolescents and young adults, of which 74 were MDD patients with SI and 150 were healthy controls. Gene-based association analyses of rare variants were performed using enrichment analysis and a cumulative allele test. An allele association study was performed against common variants.

RESULTS

After sequencing and bioinformatics analysis, a total of 15 single nucleotide variants (SNVs) were detected in the targeted regions from all participants, including 9 common and 6 rare variants. Among these, 5 rare variants were identified within the study group. Enrichment analysis of rare variants demonstrated a statistical difference (p = 0.042) between the study and control groups. Using cumulative allele analysis, alternative alleles in the SLC6A4 gene exhibited an association with MDD patients with SI (cumulative allele: OR = 10.18, 95% CI = 1.18-87.32, p = 0.017). No significant association was found between the 9 common SLC6A4 variants and MDD patients with SI.

CONCLUSIONS

Our results suggest that rare variants of SLC6A4 may contribute to a genetic risk of adolescents and young adults suffering MDD with SI.

Genetic basis of depressive disorders

Depression is a common mental disorder being one of the main causes of disability and mortality worldwide. Despite an intensive research during the past decades, the etiology of depressive disorders (DDs) remains incompletely understood; however, genetic factors are significantly involved in the liability to depression. The present review is focused on the studies based on a candidate gene approach, genome-wide association studies (GWAS) and whole exome sequencing (WES), which previously demonstrated associations between gene polymorphisms and DDs. According to the first approach, DD development is affected by serotonergic (TPH1, TPH2, HTR1A, HTR2A, and SLC6A4), dopaminergic (DRD4, SLC6A3) and noradrenergic (SLC6A2) system genes, and genes of enzymatic degradation (MAOA, COMT). In addition, there is evidence of the involvement of HPA-axis genes (OXTR, AVPR1A, and AVPR1B), sex hormone receptors genes (ESR1, ESR2, and AR), neurotrophin (BDNF) and methylenetetrahydrofolate reductase (MTHFR) genes, neuronal apoptosis (CASP3, BCL-XL, BAX, NPY, APP, and GRIN1) and inflammatory system (TNF, CRP, IL6, IL1B, PSMB4, PSMD9, and STAT3) genes in DD development. The results of the second approach (GWAS and WES) revealed that the PCLO, SIRT1, GNL3, GLT8D1, ITIH3, MTNR1A, BMP5, FHIT, KSR2, PCDH9, and AUTS2 genes predominantly responsible for neurogenesis and cell adhesion are involved in liability to depression. Therefore, the findings discussed suggest that genetic liability to DD is a complex process, which assumes simultaneous functioning of multiple genes including those reported previously, and requires future research in this field.

 

Loss of nuclear UBE3A causes electrophysiological and behavioral deficits in mice and is associated with Angelman syndrome

Mutations affecting the gene encoding the ubiquitin ligase UBE3A cause Angelman syndrome. Although most studies focus on the synaptic function of UBE3A, we show that UBE3A is highly enriched in the nucleus of mouse and human neurons. We found that the two major isoforms of UBE3A exhibit highly distinct nuclear versus cytoplasmic subcellular localization. Both isoforms undergo nuclear import through direct binding to PSMD4 (also known as S5A or RPN10), but the amino terminus of the cytoplasmic isoform prevents nuclear retention. Mice lacking the nuclear UBE3A isoform recapitulate the behavioral and electrophysiological phenotypes of Ube3a^m-/p+ mice, whereas mice harboring a targeted deletion of the cytosolic isoform are unaffected. Finally, we identified Angelman syndrome-associated UBE3A missense mutations that interfere with either nuclear targeting or nuclear retention of UBE3A. Taken together, our findings elucidate the mechanisms underlying the subcellular localization of UBE3A, and indicate that the nuclear UBE3A isoform is the most critical for the pathophysiology of Angelman syndrome.

Genome-wide association analysis uncovers variants for reproductive variation across dog breeds and links to domestication

BACKGROUND AND OBJECTIVES

The diversity of eutherian reproductive strategies has led to variation in many traits, such as number of offspring, age of reproductive maturity and gestation length. While reproductive trait variation has been extensively investigated and is well established in mammals, the genetic loci contributing to this variation remain largely unknown. The domestic dog, Canis lupus familiaris is a powerful model for studies of the genetics of inherited disease due to its unique history of domestication. To gain insight into the genetic basis of reproductive traits across domestic dog breeds, we collected phenotypic data for four traits, cesarean section rate, litter size, stillbirth rate and gestation length, from primary literature and breeders' handbooks.

METHODOLOGY

By matching our phenotypic data to genomic data from the Cornell Veterinary Biobank, we performed genome-wide association analyses for these four reproductive traits, using body mass and kinship among breeds as covariates.

RESULTS

We identified 12 genome-wide significant associations between these traits and genetic loci, including variants near CACNA2D3 with gestation length, MSRB3 and MSANTD1 with litter size, SMOC2 with cesarean section rate and UFM1 with stillbirth rate. A few of these loci, such as CACNA2D3 and MSRB3, have been previously implicated in human reproductive pathologies, whereas others have been associated with domestication-related traits, including brachycephaly (SMOC2) and coat curl (KRT71).

CONCLUSIONS AND IMPLICATIONS

We hypothesize that the artificial selection that gave rise to dog breeds also influenced the observed variation in their reproductive traits. Overall, our work establishes the domestic dog as a system for studying the genetics of reproductive biology and disease.

LAY SUMMARY

The genetic contributors to variation in mammalian reproductive traits remain largely unknown. We took advantage of the domestic dog, a powerful model system, to test for associations between genome-wide variants and four reproductive traits (cesarean section rate, litter size, stillbirth rate and gestation length) that vary extensively across breeds. We identified associations at a dozen loci, including ones previously associated with domestication-related traits, suggesting that selection on dog breeds also influenced their reproductive traits.

DNA methylation signatures of monozygotic twins clinically discordant for multiple sclerosis

Multiple sclerosis (MS) is an inflammatory, demyelinating disease of the central nervous system with a modest concordance rate in monozygotic twins, which strongly argues for involvement of epigenetic factors. We observe highly similar peripheral blood mononuclear cell-based methylomes in 45 MS-discordant monozygotic twins. Nevertheless, we identify seven MS-associated differentially methylated positions (DMPs) of which we validate two, including a region in the TMEM232 promoter and ZBTB16 enhancer. In CD4 + T cells we find an MS-associated differentially methylated region in FIRRE. Additionally, 45 regions show large methylation differences in individual pairs, but they do not clearly associate with MS. Furthermore, we present epigenetic biomarkers for current interferon-beta treatment, and extensive validation shows that the ZBTB16 DMP is a signature for prior glucocorticoid treatment. Taken together, this study represents an important reference for epigenomic MS studies, identifies new candidate epigenetic markers, and highlights treatment effects and genetic background as major confounders.

Monozygotic (MZ) twins are ideal to study the influence of non-genetic factors on complex phenotypes. Here, Souren et al. perform an EWAS in peripheral blood mononuclear cells from 45 MZ twins discordant for multiple sclerosis and identify disease and treatment-associated epigenetic markers.

Whole-Genome Linkage Scan Combined With Exome Sequencing Identifies Novel Candidate Genes for Carotid Intima-Media Thickness

Carotid intima-media thickness (cIMT) is an established heritable marker for subclinical atherosclerosis. In this study, we aim to identify rare variants with large effects driving differences in cIMT by performing genome-wide linkage analysis of individuals in the extremes of cIMT trait distribution (>90th percentile) in a large family-based study from a genetically isolated population in the Netherlands. Linked regions were subsequently explored by fine-mapping using exome sequencing. We observed significant evidence of linkage on chromosomes 2p16.3 [rs1017418, heterogeneity LOD (HLOD) = 3.35], 19q13.43 (rs3499, HLOD = 9.09), 20p13 (rs1434789, HLOD = 4.10), and 21q22.12 (rs2834949, HLOD = 3.59). Fine-mapping using exome sequencing data identified a non-coding variant (rs62165235) in PNPT1 gene under the linkage peak at chromosome 2 that is likely to have a regulatory function. The variant was associated with quantitative cIMT in the family-based study population (effect = 0.27, p-value = 0.013). Furthermore, we identified several genes under the linkage peak at chromosome 21 highly expressed in tissues relevant for atherosclerosis. To conclude, our linkage analysis identified four genomic regions significantly linked to cIMT. Further analyses are needed to demonstrate involvement of identified candidate genes in development of atherosclerosis.

A Genome-Wide Linkage Study for Chronic Obstructive Pulmonary Disease in a Dutch Genetic Isolate Identifies Novel Rare Candidate Variants

Chronic obstructive pulmonary disease (COPD) is a complex and heritable disease, associated with multiple genetic variants. Specific familial types of COPD may be explained by rare variants, which have not been widely studied. We aimed to discover rare genetic variants underlying COPD through a genome-wide linkage scan. Affected-only analysis was performed using the 6K Illumina Linkage IV Panel in 142 cases clustered in 27 families from a genetic isolate, the Erasmus Rucphen Family (ERF) study. Potential causal variants were identified by searching for shared rare variants in the exome-sequence data of the affected members of the families contributing most to the linkage peak. The identified rare variants were then tested for association with COPD in a large meta-analysis of several cohorts. Significant evidence for linkage was observed on chromosomes 15q14-15q25 [logarithm of the odds (LOD) score = 5.52], 11p15.4-11q14.1 (LOD = 3.71) and 5q14.3-5q33.2 (LOD = 3.49). In the chromosome 15 peak, that harbors the known COPD locus for nicotinic receptors, and in the chromosome 5 peak we could not identify shared variants. In the chromosome 11 locus, we identified four rare (minor allele frequency (MAF) <0.02), predicted pathogenic, missense variants. These were shared among the affected family members. The identified variants localize to genes including neuroblast differentiation-associated protein (AHNAK), previously associated with blood biomarkers in COPD, phospholipase C Beta 3 (PLCB3), shown to increase airway hyper-responsiveness, solute carrier family 22-A11 (SLC22A11), involved in amino acid metabolism and ion transport, and metallothionein-like protein 5 (MTL5), involved in nicotinate and nicotinamide metabolism. Association of SLC22A11 and MTL5 variants were confirmed in the meta-analysis of 9,888 cases and 27,060 controls. In conclusion, we have identified novel rare variants in plausible genes related to COPD. Further studies utilizing large sample whole-genome sequencing should further confirm the associations at chromosome 11 and investigate the chromosome 15 and 5 linked regions.

Progress in Genetic Studies of Tourette's Syndrome

Tourette's Syndrome (TS) is a complex disorder characterized by repetitive, sudden, and involuntary movements or vocalizations, called tics. Tics usually appear in childhood, and their severity varies over time. In addition to frequent tics, people with TS are at risk for associated problems including attention deficit hyperactivity disorder (ADHD), obsessive-compulsive disorder (OCD), anxiety, depression, and problems with sleep. TS occurs in most populations and ethnic groups worldwide, and it is more common in males than in females. Previous family and twin studies have shown that the majority of cases of TS are inherited. TS was previously thought to have an autosomal dominant pattern of inheritance. However, several decades of research have shown that this is unlikely the case. Instead TS most likely results from a variety of genetic and environmental factors, not changes in a single gene. In the past decade, there has been a rapid development of innovative genetic technologies and methodologies, as well as significant progresses in genetic studies of psychiatric disorders. In this review, we will briefly summarize previous genetic epidemiological studies of TS and related disorders. We will also review previous genetic studies based on genome-wide linkage analyses and candidate gene association studies to comment on problems of previous methodological and strategic issues. Our main purpose for this review will be to summarize the new genetic discoveries of TS based on novel genetic methods and strategies, such as genome-wide association studies (GWASs), whole exome sequencing (WES) and whole genome sequencing (WGS). We will also compare the new genetic discoveries of TS with other major psychiatric disorders in order to understand the current status of TS genetics and its relationship with other psychiatric disorders.