Jubao Duan

Common variants on chromosome 6p22.1 are associated with schizophrenia

Schizophrenia, a devastating psychiatric disorder, has a prevalence of 0.5–1%, with high heritability (80–85%) and complex transmission. Recent studies implicate rare, large, high-penetrance copy number variants in some cases, but the genes or biological mechanisms that underlie susceptibility are not known. Here we show that schizophrenia is significantly associated with single nucleotide polymorphisms (SNPs) in the extended major histocompatibility complex region on chromosome 6. We carried out a genome-wide association study of common SNPs in the Molecular Genetics of Schizophrenia (MGS) case-control sample, and then a meta-analysis of data from the MGS, International Schizophrenia Consortium and SGENE data sets. No MGS finding achieved genome-wide statistical significance. In the meta-analysis of European-ancestry subjects (8,008 cases, 19,077 controls), significant association with schizophrenia was observed in a region of linkage disequilibrium on chromosome 6p22.1 (P = 9.54 × 10-9). This region includes a histone gene cluster and several immunity-related genes—possibly implicating aetiological mechanisms involving chromatin modification, transcriptional regulation, autoimmunity and/or infection. These results demonstrate that common schizophrenia susceptibility alleles can be detected. The characterization of these signals will suggest important directions for research on susceptibility mechanisms.

Identification of loci associated with schizophrenia by genome-wide association and follow-up

We carried out a genome-wide association study of schizophrenia (479 cases, 2,937 controls) and tested loci with P < 10−5 in up to 16,726 additional subjects. Of 12 loci followed up, 3 had strong independent support (P < 5 × 10−4), and the overall pattern of replication was unlikely to occur by chance (P = 9 × 10−8). Meta-analysis provided strongest evidence for association around ZNF804A (P = 1.61 × 10−7) and this strengthened when the affected phenotype included bipolar disorder (P = 9.96 × 10−9).

Narcolepsy is strongly associated with the T-cell receptor alpha locus

Narcolepsy with cataplexy, characterized by sleepiness and rapid onset into REM sleep, affects 1 in 2,000 individuals. Narcolepsy was first shown to be tightly associated with HLA-DR2 (ref. 3) and later sublocalized to DQB1*0602 (ref. 4). Following studies in dogs and mice, a 95% loss of hypocretin-producing cells in postmortem hypothalami from narcoleptic individuals was reported. Using genome-wide association (GWA) in Caucasians with replication in three ethnic groups, we found association between narcolepsy and polymorphisms in the TRA@ (T-cell receptor alpha) locus, with highest significance at rs1154155 (average allelic odds ratio 1.69, genotypic odds ratios 1.94 and 2.55, P < 10−21, 1,830 cases, 2,164 controls). This is the first documented genetic involvement of the TRA@ locus, encoding the major receptor for HLA-peptide presentation, in any disease. It is still unclear how specific HLA alleles confer susceptibility to over 100 HLA-associated disorders; thus, narcolepsy will provide new insights on how HLA–TCR interactions contribute to organ-specific autoimmune targeting and may serve as a model for over 100 other HLA-associated disorders.

Copy Number Variants in Schizophrenia: Confirmation of Five Previous Findings and New Evidence for 3q29 Microdeletions and VIPR2 Duplications

OBJECTIVE To evaluate previously reported associations of copy number variants (CNVs) with schizophrenia and to identify additional associations, the authors analyzed CNVs in the Molecular Genetics of Schizophrenia study (MGS) and additional available data. METHOD After quality control, MGS data for 3,945 subjects with schizophrenia or schizoaffective disorder and 3,611 screened comparison subjects were available for analysis of rare CNVs (<1% frequency). CNV detection thresholds were chosen that maximized concordance in 151 duplicate assays. Pointwise and genewise analyses were carried out, as well as analyses of previously reported regions. Selected regions were visually inspected and confirmed with quantitative polymerase chain reaction. RESULTS In analyses of MGS data combined with other available data sets, odds ratios of 7.5 or greater were observed for previously reported deletions in chromosomes 1q21.1, 15q13.3, and 22q11.21, duplications in 16p11.2, and exon-disrupting deletions in NRXN1. The most consistently supported candidate associations across data sets included a 1.6-Mb deletion in chromosome 3q29 (21 genes, TFRC to BDH1) that was previously described in a mild-moderate mental retardation syndrome, exonic duplications in the gene for vasoactive intestinal peptide receptor 2 (VIPR2), and exonic duplications in C16orf72. The case subjects had a modestly higher genome-wide number of gene-containing deletions (>100 kb and >1 Mb) but not duplications. CONCLUSIONS The data strongly confirm the association of schizophrenia with 1q21.1, 15q13.3, and 22q11.21 deletions, 16p11.2 duplications, and exonic NRXN1 deletions. These CNVs, as well as 3q29 deletions, are also associated with mental retardation, autism spectrum disorders, and epilepsy. Additional candidate genes and regions, including VIPR2, were identified. Study of the mechanisms underlying these associations should shed light on the pathophysiology of schizophrenia.

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

A genome scan meta-a nalysis (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 (PSR) 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.

Genomewide Linkage Scan of 409 European-Ancestry and African American Families with Schizophrenia: Suggestive Evidence of Linkage at 8p23.3-p21.2 and 11p13.1-q14.1 in the Combined Sample

We report the clinical characteristics of a schizophrenia sample of 409 pedigrees--263 of European ancestry (EA) and 146 of African American ancestry (AA)--together with the results of a genome scan (with a simple tandem repeat polymorphism interval of 9 cM) and follow-up fine mapping. A family was required to have a proband with schizophrenia (SZ) and one or more siblings of the proband with SZ or schizoaffective disorder. Linkage analyses included 403 independent full-sibling affected sibling pairs (ASPs) (279 EA and 124 AA) and 100 all-possible half-sibling ASPs (15 EA and 85 AA). Nonparametric multipoint linkage analysis of all families detected two regions with suggestive evidence of linkage at 8p23.3-q12 and 11p11.2-q22.3 (empirical Z likelihood-ratio score [Z(lr)] threshold >/=2.65) and, in exploratory analyses, two other regions at 4p16.1-p15.32 in AA families and at 5p14.3-q11.2 in EA families. The most significant linkage peak was in chromosome 8p; its signal was mainly driven by the EA families. Z(lr) scores >2.0 in 8p were observed from 30.7 cM to 61.7 cM (Center for Inherited Disease Research map locations). The maximum evidence in the full sample was a multipoint Z(lr) of 3.25 (equivalent Kong-Cox LOD of 2.30) near D8S1771 (at 52 cM); there appeared to be two peaks, both telomeric to neuregulin 1 (NRG1). There is a paracentric inversion common in EA individuals within this region, the effect of which on the linkage evidence remains unknown in this and in other previously analyzed samples. Fine mapping of 8p did not significantly alter the significance or length of the peak. We also performed fine mapping of 4p16.3-p15.2, 5p15.2-q13.3, 10p15.3-p14, 10q25.3-q26.3, and 11p13-q23.3. The highest increase in Z(lr) scores was observed for 5p14.1-q12.1, where the maximum Z(lr) increased from 2.77 initially to 3.80 after fine mapping in the EA families.

The Role of Genetics in the Etiology of Schizophrenia

Schizophrenia is a complex genetic disorder manifesting combined environmental and genetic causation. Recently, genome-wide association experiments yielded remarkable new experimental evidence that is leading to a better understanding of the genetic models and the biological risk factors involved in schizophrenia. These studies have discovered uncommon copy number variations (mainly deletions) and common single nucleotide polymorphisms with alleles associated with schizophrenia. The aggregate data provide support for polygenic inheritance and for genetic overlap of schizophrenia with autism and with bipolar disorder. It is anticipated that the application of a myriad of tools from systems biology, in combination with biological functional experiments, will lead to a delineation of biological pathways involved in the pathophysiology of schizophrenia, and eventually to new therapies.

Mammalian Mutation Pressure, Synonymous Codon Choice, and mRNA Degradation

The usage of synonymous codons (SCs) in mammalian genes is highly correlated with local base composition and is therefore thought to be determined by mutation pressure. The usage is nonetheless structured. For instance, mammals share with Saccharomyces and Drosophila most preferences for the C-ending over the G-ending codon (or vice versa) within each fourfold-degenerate SC family and the fact that their SCs are placed along coding regions in ways that minimize the number of T|A and C|G dinucleotides (“|” being the codon boundary). TA and CG underrepresentations are observed everywhere in the mammalian genome affecting the SC usage, the amino acid composition of proteins, and the primary structure of introns and noncoding DNA. While the rarity of CG is ascribed to the high mutability of this dinucleotide, the rarity of TA in coding regions is considered adaptive because UA dinucleotides are cleaved by endoribonucleases. Here we present in vivo experimental evidence indicating that the number of T|A and/or C|G dinucleotides of a human gene can affect strongly the expression level and degradation of its mRNA. Our results are consistent with indirect evidence produced by other workers and with the detailed work that has been devoted to characterize UA cleavage in vitro and in vivo. We conclude that SC choice can influence strongly mRNA function and gene expression through effects not directly related to the codon–anticodon interaction. These effects should constrain heavily the nucleotide motif composition of the most abundant mRNAs in the transcriptome, in particular, their SC usage, a usage that must be reflected by cellular tRNA concentrations and thus defines for all other genes which SCs are translated fastest and most accurately. Furthermore, the need to avoid such effects genome-wide appears serious enough to have favored the evolution of biases in context-dependent mutation that reduce the occurrence of intrinsically unfavorable motifs, and/or, when possible, to have induced the molecular machinery mediating such effects to rely opportunistically on already existing motif rarities and abundances. This may explain why nucleotide motif preferences are very similar in transcribed and nontranscribed mammalian DNA even though the preferences appear to be adaptive only in transcribed DNA.

Polymorphisms in the Trace Amine Receptor 4 (TRAR4) Gene on Chromosome 6q23.2 Are Associated with Susceptibility to Schizophrenia

Several linkage studies across multiple population groups provide convergent support for a susceptibility locus for schizophrenia--and, more recently, for bipolar disorder--on chromosome 6q13-q26. We genotyped 192 European-ancestry and African American (AA) pedigrees with schizophrenia from samples that previously showed linkage evidence to 6q13-q26, focusing on the MOXD1-STX7-TRARs gene cluster at 6q23.2, which contains a number of prime candidate genes for schizophrenia. Thirty-one screening single-nucleotide polymorphisms (SNPs) were selected, providing a minimum coverage of at least 1 SNP/20 kb. The association observed with rs4305745 (P=.0014) within the TRAR4 (trace amine receptor 4) gene remained significant after correction for multiple testing. Evidence for association was proportionally stronger in the smaller AA sample. We performed database searches and sequenced genomic DNA in a 30-proband subsample to obtain a high-density map of 23 SNPs spanning 21.6 kb of this gene. Single-SNP analyses and also haplotype analyses revealed that rs4305745 and/or two other polymorphisms in perfect linkage disequilibrium (LD) with rs4305745 appear to be the most likely variants underlying the association of the TRAR4 region with schizophrenia. Comparative genomic analyses further revealed that rs4305745 and/or the associated polymorphisms in complete LD with rs4305745 could potentially affect gene expression. Moreover, RT-PCR studies of various human tissues, including brain, confirm that TRAR4 is preferentially expressed in those brain regions that have been implicated in the pathophysiology of schizophrenia. These data provide strong preliminary evidence that TRAR4 is a candidate gene for schizophrenia; replication is currently being attempted in additional clinical samples.

Polymorphisms in the 5'-untranslated region of the human serotonin receptor 1B (HTR1B) gene affect gene expression.

We present evidence of complex balancing regulation of HTR1B transcription by common polymorphisms in its promoter. Computational analysis of the HTR1B gene predicted that a 5′ segment, spanning common DNA sequence variations, T−261G, A−161T, and −182INS/DEL−181, contained a putative functional promoter. Using a secreted alkaline phosphatase (SEAP) reporter gene system, we found that the haplotype −261G_−182INS−181_A−161 enhanced transcriptional activity 2.3-fold compared with the haplotype T−261_−182INS−181_A−161. Conversely, −161T reversed this, and the net effect when −261G and −161T were in the same haplotype (−261G_−182INS−181_−161T) was equivalent to the major haplotype (T−261_−182INS−181_A−161). Electrophoretic mobility shift experiments showed that −261G and −161T modify the binding of transcription factors (TFs): −261G generates a new AP2 binding site, while alleles A−161 and −161T exhibit different binding characteristics to AP1. T−261G and A−161T were found to be in linkage disequilibrium (LD) with G861C in a European ancestry population. Interestingly, G861C has been reported to be associated with several psychiatric disorders. Our results indicate that HTR1B is the target of substantial transcriptional genetic regulation by common haplotypes, which are in LD with the HTR1B single-nucleotide polymorphism (SNP) most commonly used in association studies.