Genome scan of schizophrenia families in a large Veterans Affairs Cooperative Study sample: evidence for linkage to 18p11.32 and for racial heterogeneity on chromosomes 6 and 14

High resolution chromosomal microarray analysis in paediatric obsessive-compulsive disorder

Background

Obsessive-Compulsive Disorder (OCD) is a common and chronic disorder in which a person has uncontrollable, reoccurring thoughts and behaviours. It is a complex genetic condition and, in case of early onset (EO), the patients manifest a more severe phenotype, and an increased heritability. Large (>500 kb) copy number variations (CNVs) previously associated with autism and schizophrenia have been reported in OCD. Recently, rare CNVs smaller than 500 kb overlapping risk loci for other neurodevelopmental conditions have also been reported in OCD, stressing the importance of examining CNVs of any size range. The aim of this study was to further investigate the role of rare and small CNVs in the aetiology of EO-OCD.

Methods

We performed high-resolution chromosomal microarray analysis in 121 paediatric OCD patients and in 124 random controls to identify rare CNVs (>50 kb) which might contribute to EO-OCD.

Results

The frequencies and the size of the observed rare CNVs in the patients did not differ from the controls. However, we observed a significantly higher frequency of rare CNVs affecting brain related genes, especially deletions, in the patients (OR = 1.98, 95% CI 1.02–3.84; OR = 3.61, 95% CI 1.14–11.41, respectively). Similarly, enrichment-analysis of CNVs gene content, performed with three independent methods, confirmed significant clustering of predefined genes involved in synaptic/brain related functional pathways in the patients but not in the controls. In two patients we detected de-novo CNVs encompassing genes previously associated with different neurodevelopmental disorders (NRXN1, ANKS1B, UHRF1BP1).

Conclusions

Our results further strengthen the role of small rare CNVs, particularly deletions, as susceptibility factors for paediatric OCD.

Electronic supplementary material

The online version of this article (10.1186/s12920-017-0299-5) contains supplementary material, which is available to authorized users.

Neuropeptides shaping the central nervous system development: Spatiotemporal actions of VIP and PACAP through complementary signaling pathways

Pituitary adenylate cyclase-activating polypeptide (PACAP) and vasoactive intestinal peptide (VIP) are neuropeptides with wide, complementary, and overlapping distributions in the central and peripheral nervous systems, where they exert important regulatory roles in many physiological processes. VIP and PACAP display a large range of biological cellular targets and functions in the adult nervous system including regulation of neurotransmission and neuroendocrine secretion and neuroprotective and neuroimmune responses. As the main focus of the present review, VIP and PACAP also have been long implicated in nervous system development and maturation through their interaction with the seven transmembrane domain G protein-coupled receptors, PAC1, VPAC1, and VPAC2, initiating multiple signaling pathways. Compared with PAC1, which solely binds PACAP with very high affinity, VPACs exhibit high affinities for both VIP and PACAP but differ from each other because of their pharmacological profile for both natural accessory peptides and synthetic or chimeric molecules, with agonistic and antagonistic properties. Complementary to initial pharmacological studies, transgenic animals lacking these neuropeptides or their receptors have been used to further characterize the neuroanatomical, electrophysiological, and behavioral roles of PACAP and VIP in the developing central nervous system. In this review, we recapitulate the critical steps and processes guiding/driving neurodevelopment in vertebrates and superimposing the potential contribution of PACAP and VIP receptors on the given timeline. We also describe how alterations in VIP/PACAP signaling may contribute to both (neuro)developmental and adult pathologies and suggest that tuning of VIP/PACAP signaling in a spatiotemporal manner may represent a novel avenue for preventive therapies of neurological and psychiatric disorders. © 2016 Wiley Periodicals, Inc.

Animal models of gene-environment interaction in schizophrenia: A dimensional perspective

Schizophrenia has long been considered as a disorder with multifactorial origins. Recent discoveries have advanced our understanding of the genetic architecture of the disease. However, even with the increase of identified risk variants, heritability estimates suggest an important contribution of non-genetic factors. Various environmental risk factors have been proposed to play a role in the etiopathogenesis of schizophrenia. These include season of birth, maternal infections, obstetric complications, adverse events at early childhood, and drug abuse. Despite the progress in identification of genetic and environmental risk factors, we still have a limited understanding of the mechanisms whereby gene-environment interactions (G × E) operate in schizophrenia and psychoses at large. In this review we provide a critical analysis of current animal models of G × E relevant to psychotic disorders and propose that dimensional perspective will advance our understanding of the complex mechanisms of these disorders.

PACAP and PAC1 receptor in brain development and behavior

Pituitary adenylate cyclase-activating polypeptide (PACAP) and vasoactive intestinal peptide (VIP) act through three class B G-protein coupled receptors, PAC1, VPAC1 and VPAC2, initiating multiple signaling pathways. In addition to natural peptides ligands, a number of synthetic peptides and a small molecular antagonist have been generated. Genetically modified animals have been produced for the neuropeptides and receptors. Neuroanatomical, electrophysiological, behavioral and pharmacological characterization of the mutants and transgenic mice uncovered diverse roles of PACAP-PAC1-VAPC2 signaling in peripheral tissues and in the central nervous system. Human genetic studies suggest that the PACAP-PAC1-VPAC2 signaling can be associated with psychiatric illness via mechanisms of not only loss-of-function, but also gain-of-function. For example, a duplication of chromosome 7q36.3 (encoding the VPAC2 receptor) was shown to be associated with schizophrenia, and high levels of PACAP-PAC1 signaling are associated with posttraumatic stress disorder. Whereas knockout animals are appropriate to address loss-of-function of human genetics, transgenic mice overexpressing human transgenes in native environment using artificial chromosomes are particularly valuable and essential to address the consequences of gain-of-function. This review focuses on role of PACAP and PAC1 receptor in brain development, behavior of animals and potential implication in human neurodevelopmental disorders. It also encourages keeping an open mind that alterations of VIP/PACAP signaling may associate with psychiatric illness without overt neuroanatomic changes, and that tuning of VIP/PACAP signaling may represent a novel avenue for the treatment of the psychiatric illness.

Schizophrenia: from epidemiology to rehabilitation

Purpose/Objective:

We discuss recent evidences about schizophrenia (frequency, onset, course, risk factors and genetics) and their influences to some epidemiological myths about schizophrenia diffuse between psychiatric and psychopathology clinicians. The scope is to evaluate if the new acquisitions may change the rehabilitation approaches to schizophrenia modifying the balance about the neurodevelopmental hypothesis of schizophrenia accepting that the cognitive deficits are produced by errors during the normal development of the brain (neurodevelopmental hypothesis) that remains stable in the course of illness and the neurodegenerative hypothesis according of which they derived from a degenerative process that goes on inexorably.

Research Method/Design:

A review of the literature about epidemiology of schizophrenia has been performed and the contributions of some of these evidence to neurodevelopmental hypothesis and to rehabilitation has been described.

Results:

It cannot be definitively concluded for or against the neurodevelopmental or degenerative hypothesis, but efforts in understanding basis of schizophrenia must go on. Until now, rehabilitation programs are based on the vulnerability-stress model: supposing an early deficit that go on stable during the life under favorable circumstances. So, rehabilitation approaches (as neuro-cognitive approaches, social skill training, cognitive-emotional training) are focused on the individual and micro-group coping skills, aiming to help people with schizophrenia to cope with environmental stress factors.

Conclusions/Implications:

Coping of cognitive deficits in schizophrenia may represents the starting-point for further research on schizophrenia, cohort studies and randomized trials are necessary to defined the range of effectiveness and the outcome of the treatments.

Meta-analysis of 32 genome-wide linkage studies of schizophrenia

A genome scan meta-analysis (GSMA) was carried out on 32 independent genome-wide linkage scan analyses that included 3255 pedigrees with 7413 genotyped cases affected with schizophrenia (SCZ) or related disorders. The primary GSMA divided the autosomes into 120 bins, rank-ordered the bins within each study according to the most positive linkage result in each bin, summed these ranks (weighted for study size) for each bin across studies and determined the empirical probability of a given summed rank (P(SR)) by simulation. Suggestive evidence for linkage was observed in two single bins, on chromosomes 5q (142-168 Mb) and 2q (103-134 Mb). Genome-wide evidence for linkage was detected on chromosome 2q (119-152 Mb) when bin boundaries were shifted to the middle of the previous bins. The primary analysis met empirical criteria for 'aggregate' genome-wide significance, indicating that some or all of 10 bins are likely to contain loci linked to SCZ, including regions of chromosomes 1, 2q, 3q, 4q, 5q, 8p and 10q. In a secondary analysis of 22 studies of European-ancestry samples, suggestive evidence for linkage was observed on chromosome 8p (16-33 Mb). Although the newer genome-wide association methodology has greater power to detect weak associations to single common DNA sequence variants, linkage analysis can detect diverse genetic effects that segregate in families, including multiple rare variants within one locus or several weakly associated loci in the same region. Therefore, the regions supported by this meta-analysis deserve close attention in future studies.

Linkage analysis of schizophrenia in African-American families

While many studies have sought a window into the genetics of schizophrenia, few have focused on African-American families. An exception is the Project among African-Americans to Explore Risks for Schizophrenia (PAARTNERS), which seeks to identify novel and known risk variation for schizophrenia by genetic analyses of African-American families. We report a linkage study of diagnostic status in 217 African-American families using the Illumina Linkage Panel. Due to assumed incomplete and time-dependent penetrance, we performed linkage analysis using two different treatments of diagnosis: (1) treating both affected and unaffected individuals as informative for linkage (using the program SIBPAL) and (2) treating only affected individuals as informative (using the program MERLIN). We also explore three definitions of affected status: narrowly defined schizophrenia; one broadened to include schizoaffective disorder; and another including all diagnoses indicating psychosis. Several regions show a decrease in the evidence for linkage as the definition broadens 8q22.1 (rs911, 99.26 cM; SIBPAL p-value [p] goes from 0.006 to 0.02), 16q24.3 (rs1006547, 130.48 cM; p from 0.00095 to 0.0085), and 20q13.2 (rs1022689, 81.73 cM; p from 0.00015 to 0.032). One region shows a substantial increase in evidence for linkage, 11p15.2 (rs722317, 24.27 cM; p from 0.0022 to 0.0000003); MERLIN results support the significance of the SIBPAL results (p=0.00001). Our linkage results overlap two broad, previously-reported linkage regions: 8p23.3-p12 found in studies sampling largely families of European ancestry; and 11p11.2-q22.3 reported by a study of African-American families. These results should prove quite useful for uncovering loci affecting risk for schizophrenia.

Copy number variations and risk for schizophrenia in 22q11.2 deletion syndrome

22q11.2 Deletion Syndrome (22q11.2DS) is a common microdeletion syndrome with congenital and late-onset features. Testing for the genomic content of copy number variations (CNVs) may help elucidate the 22q11.2 deletion mechanism and the variable clinical expression of the syndrome including the high (25%) risk for schizophrenia. We used genome-wide microarrays to assess CNV content and the parental origin of 22q11.2 deletions in a cohort of 100 adults with 22q11.2DS (44 with schizophrenia) and controls. 22q11.2DS subjects with schizophrenia failed to exhibit de novo CNVs or any excess of novel inherited CNVs outside the 22q11.2 region. There were no significant effects of parental origin of the 22q11.2 deletion, deletion length, parental age or family history on expression of schizophrenia. There was no evidence for a general increase of de novo CNVs in 22q11.2DS. A novel finding was the relative paucity of males with de novo 22q11.2 deletions of paternal origin (P = 0.019). The Y chromosome may play a mediating role in the mechanism of 22q11.2 deletion events during spermatogenesis, resulting in the previously observed excess of maternal de novo 22q11.2 deletions. Hemizygosity of the 22q11.2 region appears to be the major CNV-related risk factor for schizophrenia in 22q11.2DS. The results reinforce the need for further efforts to identify specific molecular mechanisms underlying this expression and to identify the 1% of patients with schizophrenia who carry 22q11.2 deletions.

Association of SNPs and haplotypes in APOL1, 2 and 4 with schizophrenia

Prior work found the APOL1, 2 and 4 genes, located on chromosome 22q12.3-q13.1, to be upregulated in brains of schizophrenic patients. We performed a family-based association study using 130 SNPs tagging the APOL gene family (APOL1-6). The subjects were 112 African-American (AA), 114 European-American (EA), 109 Chinese (Ch) and 42 Japanese (Jp) families with schizophrenia (377 families, 1161 genotyped members and 647 genotyped affected in total). Seven SNPs had p-values<0.05 in the APOL1, 2 and 4 regions for the AA, EA and combined (AA and EA) samples. In the AA sample, two SNPs, rs9610449 and rs6000200 showed low p-values; and a haplotype which comprised these two SNPs yielded a p-value of 0.00029 using the global test (GT) and the allele specific test (AST). The two SNPs and the haplotype were associated with risk for schizophrenia in African-Americans. In the combined (AA and EA) sample, two SNPs, rs2003813 and rs2157249 showed low p-values; and a three SNP haplotype including these two SNPs was significant using the GT (p=0.0013) and the AST (p=0.000090). The association of this haplotype with schizophrenia was significant for the entire (AA, EA, Ch and Jp) sample using the GT (p=0.00054) and the AST (p=0.00011). Although our study is not definitive, it suggests that the APOL genes should be more extensively studied in schizophrenia.

The myelin-pathogenesis puzzle in schizophrenia: a literature review

Schizophrenia is a serious and disabling mental disorder with symptoms such as auditory hallucinations, disordered thinking and delusions, avolition, anhedonia, blunted affect and apathy. In this review article we seek to present the current scientific findings from linkage studies and susceptible genes and the pathophysiology of white matter in schizophrenia. The article has been reviewed in two parts. The first part deals with the linkage studies and susceptible genes in schizophrenia in order to have a clear-cut picture of the involvement of chromosomes and their genes in schizophrenia. The genetic linkage results seem to be replicated in some cases but in others are not. From these results, we cannot draw a fine map to a single locus or gene, leading to the conclusion that schizophrenia is not caused by a single factor/gene. In the second part of the article we present the oligodendrocyte-related genes that are associated with schizophrenia, as we hypothesize a potential role of oligodendrocyte-related genes in the pathology of the disorder.

TGFB-induced factor (TGIF): a candidate gene for psychosis on chromosome 18p

Schizophrenia (SC) and bipolar disorder (BP) share many clinical features, among them psychosis. We previously identified a putative gene locus for psychosis on chromosome 18p in a sample from the Central Valley of Costa Rica (CVCR) population. The present study replicated the association to a specific allele of microsatellite marker D18S63 on 18p11.3, using a newly collected sample from the CVCR. A combined analysis of both samples, plus additional subjects, showed that this specific allele on D18S63, which lies within an intron on the TGFB-induced factor (TGIF) gene, is strongly associated (P-value=0.0005) with psychosis. Eleven additional SNP markers, spanning five genes in the region, were analyzed in the combined sample from the CVCR. Only the four SNPs within the TGIF gene were in strong linkage disequilibrium with D18S63 (D'=1.00). A specific haplotype for all five markers within the TGIF gene showed evidence of association (P-value=0.011) to psychosis. A second, distinct haplotype, containing a newly identified nonsynonymous polymorphism in exon 5 of the TGIF gene, showed a nonsignificant trend towards association to psychosis (P-value=0.077). TGIF is involved in neurodevelopment, neuron survival and controls the expression of dopamine receptors. Altogether, our results point to the possible involvement of TGIF in the pathophysiology of psychotic disorders in the CVCR population.

Cerebellar cortical-layer-specific control of neuronal migration by pituitary adenylate cyclase-activating polypeptide

Migration of immature neurons is essential for forming the cortical layers and nuclei. Impairment of migration results in aberrant neuronal cytoarchitecture, which leads to various neurological disorders. Neurons alter the mode, tempo and rate of migration when they translocate through different cortical layers, but little is known about the mechanisms underlying this process. Here we show that endogenous pituitary adenylate cyclase-activating polypeptide (PACAP) has short-term and cortical-layer-specific effects on granule cell migration in the early postnatal mouse cerebellum. Application of exogenous PACAP significantly slowed the migration of isolated granule cells and shortened the leading process in the microexplant cultures of the postnatal day (P)0-3 cerebella. Interestingly, in the cerebellar slices of P10 mice, application of exogenous PACAP significantly inhibited granule cell migration in the external granular layer (EGL) and molecular layer (ML), but failed to alter the movement in the Purkinje cell layer (PCL) and internal granular layer (IGL). In contrast, application of PACAP antagonist accelerated granule cell migration in the PCL, but did not change the movement in the EGL, ML and IGL. Inhibition of the cAMP signaling and the activity of phospholipase C significantly reduced the effects of exogenous PACAP on granule cell migration. The PACAP action on granule cell migration was transient, and lasted for approximately 2 h. The duration of PACAP action on granule cell migration was determined by the desensitization of its receptors and prolonged by inhibiting the protein kinase C. Endogenous PACAP was present sporadically in the bottom of the ML, intensively in the PCL, and throughout the IGL. Collectively, these results indicated that PACAP acts on granule cell migration as "a brake (stop signal) for cell movement." Furthermore, these results suggest that endogenous PACAP slows granule cell migration when the cells enter the PACAP-rich PCL, and 2 h later the desensitization of PACAP receptors allows the cells to accelerate the rate of migration and to actively move within the PACAP-rich IGL. Therefore, endogenous PACAP may provide a cue that regulates granule cell migration in a cerebellar cortical-layer-specific manner.

A genome-wide scan for schizophrenia and psychosis susceptibility loci in families of Mexican and Central American ancestry

Schizophrenia is a complex psychiatric disorder, likely to be caused in part by multiple genes. In this study, linkage analyses were performed to identify chromosomal regions most likely to be associated with schizophrenia and psychosis in multiplex families of Mexican and Central American origin. Four hundred and fifty-nine individuals from 99 families, containing at least two siblings with hospital diagnoses of schizophrenia or schizoaffective disorder, were genotyped. Four hundred and four microsatellite markers were genotyped for all individuals and multipoint non-parametric linkage analyses were performed using broad (any psychosis) and narrow (schizophrenia and schizoaffective disorder) models. Under the broad model, three chromosomal regions (1pter-p36, 5q35, and 18p11) exhibited evidence of linkage with non-parametric lod (NPL) scores greater than 2.7 (equivalent to empirical P values of less than 0.001) with the peak multipoint NPL = 3.42 (empirical P value = 0.00003), meeting genomewide evidence for significant linkage in the 1pter-p36 region. Under the narrow model, the same three loci showed (non-significant) evidence of linkage. These linkage findings (1pter-p36, 18p11, and 5q35) highlight where genes for psychosis and schizophrenia are most likely to be found in persons of Mexican and Central American ancestry, and correspond to recent linkages of schizophrenia or psychosis in other populations which were formed in part from emigrants from the Spanish empire of the 15th and 16th centuries.

Expression of the human PAC1 receptor leads to dose-dependent hydrocephalus-related abnormalities in mice

Hydrocephalus is a common and potentially devastating birth defect affecting the CNS, and its relationship with G protein-coupled receptors (GPCRs) is unknown. We have expressed 2, 4, or 6 copies of a GPCR--the human PAC1 receptor with a 130-kb transgene in the mouse nervous system in a pattern closely resembling that of the endogenous gene. Consistent with PAC1 actions, PKA and PKC activity were elevated in the brains of Tg mice. Remarkably, Tg mice developed dose-dependent hydrocephalus-like characteristics, including enlarged third and lateral ventricles and reduced cerebral cortex, corpus callosum, and subcommissural organ (SCO). Neuronal proliferation and apoptosis were implicated in hydrocephalus, and we observed significantly reduced neuronal proliferation and massively increased neuronal apoptosis in the developing cortex and SCO of Tg embryos, while neurite outgrowth and neuronal migration in vitro remain uncompromised. Ventricular ependymal cilia are crucial for directing cerebrospinal fluid flow, and ependyma of Tg mice exhibited disrupted cilia with increased phospho-CREB immunoreactivity. These data demonstrate that altered neuronal proliferation/apoptosis and disrupted ependymal cilia are the main factors contributing to hydrocephalus in PAC1-overexpressing mice. This is the first report to our knowledge demonstrating that misregulation of GPCRs can be involved in hydrocephalus-related neurodevelopmental disorders.

An evaluation of the replicate pool method: quick estimation of genome-wide linkage peak p-values

The calculation of empirical p-values for genome-wide non-parametric linkage tests continues to present significant computational challenges for many complex disease mapping studies. The gold standard approach is to use gene dropping to simulate null genome scans. Unfortunately, this approach is too computationally expensive for many data sets of interest. An alternative, more efficient method for sampling null genome scans is to pre-calculate pools of family-specific statistics and then resample from these replicate pools to generate "pseudo-replicate" genome scans. In this study, we use simulations to explore properties of the replicate pool p-value estimator pRP and show that it provides an excellent approximation to the traditional gene-dropping estimator for significantly less computational effort. While the computational efficiency of the replicate pool estimator is noticeable in almost all data sets, by applying the replicate pool method to several previously characterized data sets we show that savings in computational effort can be especially significant (on the order of 10,000-fold or more) when one or more large families are analyzed. We also estimate replicate pool p-values for the schizophrenia data described by Abecasis et al. and show that pRP closely approximates gene-drop p-values for all linkage peaks reported for this study. Lastly, we expand upon Song et al.'s previous work by deriving a conservative estimator of the variance for PRP that can easily be computed in practical settings. We have implemented the replicate pool method along with our variance estimator in a new program called Pseudo, which is the first widely available automated implementation of the replicate pool method.

An algorithm to construct genetically similar subsets of families with the use of self-reported ethnicity information

We present a simple algorithm that uses self-reported ethnicity information, pedigree structure, and affection status to group families into genetically more homogeneous subsets. This algorithm should prove useful to researchers who wish to perform genetic analyses on more-homogeneous subsets when they suspect that ignoring heterogeneity could lead to false-positive results or loss of power. We applied our algorithm to the self-reported ethnicity information of 159 families from the Veterans Affairs Cooperative Study of schizophrenia. We compared these estimates of population membership with those obtained using the program structure in an analysis of 378 microsatellite markers. We found excellent concordance between family classifications determined using self-reported ethnicity information and our algorithm and those determined using genetic marker data and structure; 158 of the 159 families had concordant classifications. In addition, the degree of admixture estimated using our algorithm and self-reported ethnicity information correlated well with that predicted using the genotype information.