Chenxing Liu

MirSNP, a database of polymorphisms altering miRNA target sites, identifies miRNA-related SNPs in GWAS SNPs and eQTLs

BackgroundNumerous single nucleotide polymorphisms (SNPs) associated with complex diseases have been identified by genome-wide association studies (GWAS) and expression quantitative trait loci (eQTLs) studies. However, few of these SNPs have explicit biological functions. Recent studies indicated that the SNPs within the 3’UTR regions of susceptibility genes could affect complex traits/diseases by affecting the function of miRNAs. These 3’UTR SNPs are functional candidates and therefore of interest to GWAS and eQTL researchers.DescriptionWe developed a publicly available online database, MirSNP (http://cmbi.bjmu.edu.cn/mirsnp), which is a collection of human SNPs in predicted miRNA-mRNA binding sites. We identified 414,510 SNPs that might affect miRNA-mRNA binding. Annotations were added to these SNPs to predict whether a SNP within the target site would decrease/break or enhance/create an miRNA-mRNA binding site. By applying MirSNP database to three brain eQTL data sets, we identified four unreported SNPs (rs3087822, rs13042, rs1058381, and rs1058398), which might affect miRNA binding and thus affect the expression of their host genes in the brain. We also applied the MirSNP database to our GWAS for schizophrenia: seven predicted miRNA-related SNPs (p < 0.0001) were found in the schizophrenia GWAS. Our findings identified the possible functions of these SNP loci, and provide the basis for subsequent functional research.ConclusionMirSNP could identify the putative miRNA-related SNPs from GWAS and eQTLs researches and provide the direction for subsequent functional researches.

Further evidence for genetic association of CACNA1C and schizophrenia: New risk loci in a Han Chinese population and a meta-analysis

CACNA1C (12p13.3) has been implicated as a susceptibility gene for schizophrenia by several replicated genome wide association studies. While these results have been consistent among studies in European populations, the findings in East Asian populations have varied. To test whether CACNA1C is a risk gene for schizophrenia, we conducted a case-control study in 5897 schizophrenic patients and 6323 healthy control subjects selected from Han Chinese population. Our study replicated the positive associations of rs1006737 (P=0.0108, OR=1.16, 95% CI: 1.03-1.29) and rs1024582 (P=0.0062, OR=1.18, 95% CI: 1.05-1.33), and identified a novel risk locus, rs2007044 (P=0.0053, OR=1.08, 95% CI: 1.02-1.14). A meta-analysis of rs1006737 combining our study and previous studies was conducted in a total of 8222 schizophrenia cases and 24,661 healthy controls. In the meta-analysis, the association between rs1006737 and schizophrenia remained significant (OR=1.14, 95% CI: 1.07-1.22, P=0.0001). Stratified analysis showed no heterogeneity between East Asian and European ancestries (χ(2)[1]=0.07, P=0.795), and the difference in pooled ORs between ancestries was not significant (Z=0.25, P=0.801). Our results provide further support for associations of rs1006737 and rs1024582 with schizophrenia, identify a new risk locus rs2007044 in a Han Chinese population, and further establish CACNA1C as an important susceptibility gene for the disease across world populations.

Sequencing ASMT Identifies Rare Mutations in Chinese Han Patients with Autism

Melatonin is involved in the regulation of circadian and seasonal rhythms and immune function. Prior research reported low melatonin levels in autism spectrum disorders (ASD). ASMT located in pseudo-autosomal region 1 encodes the last enzyme of the melatonin biosynthesis pathway. A previous study reported an association between ASD and single nucleotide polymorphisms (SNPs) rs4446909 and rs5989681 located in the promoter of ASMT. Furthermore, rare deleterious mutations were identified in a subset of patients. To investigate the association between ASMT and autism, we sequenced all ASMT exons and its neighboring region in 398 Chinese Han individuals with autism and 437 healthy controls. Although our study did not detect significant differences of genotypic distribution and allele frequencies of the common SNPs in ASMT between patients with autism and healthy controls, we identified new rare coding mutations of ASMT. Among these rare variants, 4 were exclusively detected in patients with autism including a stop mutation (p.R115W, p.V166I, p.V179G, and p.W257X). These four coding variants were observed in 6 of 398 (1.51%) patients with autism and none in 437 controls (Chi-Square test, Continuity Correction p = 0.032, two-sided). Functional prediction of impact of amino acid showed that p.R115W might affect protein function. These results indicate that ASMT might be a susceptibility gene for autism. Further studies in larger samples are needed to better understand the degree of variation in this gene as well as to understand the biochemical and clinical impacts of ASMT/melatonin deficiency.

Association study of NRXN3 polymorphisms with schizophrenia and risperidone-induced bodyweight gain in Chinese Han population

Recent researches have implicated that mutations in the neurexin-3 (NRXN3) gene on chromosome 14q24.3-q31.1 might play a role in addiction, autism, and obesity. In order to explore the association of NRXN3 polymorphisms with schizophrenia, we examined seven single nucleotide polymorphisms (SNPs) in NRXN3 spanning 1.33 Mb of this gene, in a Chinese Han sample of 1214 schizophrenic patients and 1517 healthy control subjects. Our results showed that three SNPs were associated with schizophrenia (rs7157669: A>C, p=0.006; rs724373: C>T, p=0.014; rs7154021: C>T, p=0.018). After being corrected for multiple tests, the association of rs7157669 remained significant but those for two others were modest. According to the linkage disequilibrium pattern, the 7 SNPs may construct 3 haplotype blocks. Several haplotypes were significantly associated with schizophrenia, constructed by rs11624704-rs7157669-rs724373 (AAC, p=0.003; ACT, p=0.007, both remained significant after permutation tests), rs7154021-rs7142344 (TT, p=0.024; CT, p=0.012), respectively. Among the patients, 326 ones at first onset have received 6-week monotherapy of risperidone. Further analyses showed that two SNPs were associated with percentage of bodyweight gain following a 6-week therapy of risperidone (rs11624704: p=0.03; rs7154021: p=0.008) and rs7154021 remained significant after permutation test. Our findings suggested that NRXN3 might represent a major susceptibility gene for schizophrenia and have a role in bodyweight gain related to therapy of risperidone in Chinese Han population.

Protein-interaction-network-based analysis for genome-wide association analysis of schizophrenia in Han Chinese population

Schizophrenia is a severe neuropsychiatric disorder with a strong and complex genetic background. Recent genome-wide association studies (GWAS) have successfully identified several susceptibility loci of schizophrenia. In order to interpret the functional role of the genetic variants and detect the combined effects of some of these genes on schizophrenia, protein-interaction-network-based analysis (PINBA) has emerged as an effective approach. In the current study, we conducted a PINBA of our previous GWAS data taken from the Han Chinese population. In order to do so, we used dense module search (DMS), a method that locates densely connected modules for complex diseases by integrating the association signal from GWAS datasets into the human protein-protein interaction (PPI) network. As a result, we identified one gene set with a joint effect significantly associated with schizophrenia and gene expression profiling analysis suggested that they were mainly neuro- and immune-related genes, such as glutamatergic gene (GRM5), GABAergic genes (GABRB1, GABARAP) and genes located in the MHC region (HLA-C, TAP2, HIST1H1B). Further pathway enrichment analysis suggested that these genes are involved in processes related to neuronal and immune systems, such as the Adherens junction pathway, the Neurotrophin signaling pathway and the Toll-like receptor signaling pathway. In our study, we identified a set of susceptibility genes that had been missed in single-marker GWAS, and our findings could promote the study of the genetic mechanisms in schizophrenia.

No Association of Catechol-O-Methyltransferase Polymorphisms with Schizophrenia in the Han Chinese Population

AIMS Genetics play a major role in the etiology of schizophrenia (SZ). Catechol-O-methyltransferase (COMT) is one of the promising candidate genes for SZ. A nonsynonymous single-nucleotide polymorphism (SNP), rs4680, causing a Valine (Val) to Methionine (Met) substitution, has been widely studied in relation to psychiatric phenotypes, including SZ, but with conflicting results. We conducted a two-stage study to examine the association of COMT polymorphisms with SZ in the Han Chinese population. RESULTS Association analysis of nine SNPs in 768 patients and 1348 controls failed to detect any positive markers or haplotypes. Then, we tested rs4680 in a validation sample of 963 patients and 992 controls, and no significant association was observed, but the cases significantly deviated from Hardy-Weinberg equilibrium (p=5.7e-4). There was no association of rs4680 with SZ in the combined sample (n=4071, p=0.110, odds ratio=1.08). CONCLUSIONS Our results do not support the association of COMT with SZ in the Han Chinese population.

Systematic association analysis of microRNA machinery genes with schizophrenia informs further study

microRNAs (miRNAs) play a vital role in development via the post-transcriptional regulation of most genes. Variation in the miRNA machinery pathway proteins which mediate the biogenesis, maturation, transportation, and functioning of miRNAs might be relevant to human traits. In this work, we explored the role of 59 miRNA machinery genes in schizophrenia (SZ). Association analysis of 967 single nucleotide polymorphisms within these genes detected that an intronic polymorphism of EIF4ENIF1, rs7289941, was significantly associated with SZ (P=4.10E-5). We failed to replicate this result in a validation sample comprising 1027 healthy controls and 1012 SZ cases, and the combined data yielded nominal significance (P=0.013). We conducted a gene-based association analysis using VEGAS and SKAT, and found seven associated genes in total, including EIF4ENIF1, PIWIL2, and DGCR8, but none survived correction for multiple testing. Taken together, our data do not provide strong support for the association of common variants within miRNA machinery genes with SZ in the Han Chinese population, but implicate several promising candidate genes for further research.

Association analysis of a functional variant in ATXN2 with schizophrenia

Schizophrenia (SZ) is a severe mental disorder characterized by multiple neurodevelopmental dysfunctions including a breakdown of thinking process and a deficit of typical emotional responses. Ataxin-2 (ATXN2) plays vital roles in cell proliferation and growth, and functional mutations of ATXN2 cause neurodegenerative phenotypes, including spinocerebellar ataxia type 2 (SCA2) and amyotrophic lateral sclerosis (ALS). To explore the possible role of ATXN2 in SZ, we conducted a two-stage study to examine the association of ATXN2 polymorphisms with SZ in the Han Chinese population. Association analysis of seven SNPs in 768 patients and 1348 controls revealed two associated SNPs, including rs630511 (P=1.76E-4) and rs7969300 (P=5.08E-4). We examined these two SNPs in a validation sample of 1957 patients and 1509 controls, and observed an association of rs7969300 with SZ (P=5.03E-3). The SNP rs7969300 is a non-synonymous SNP causing a Ser to Asn substitution, which is predicted to increase the protein stability of ATXN2. Our data suggest that the ATXN2 gene may confer vulnerability for SZ, adding further evidence for the genetic variants within the developmental pathway in the illness.

Myosin Vb gene is associated with schizophrenia in Chinese Han population

Myosin Vb (MYO5B) has recently been implicated in the etiology of bipolar disorder in a genome-wide association study (GWAS). This gene is involved in amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) receptor subunit glutamate receptor 1 (GluR1) recycling and plays an important role in the primary excitatory neurotransmission. Dysfunction of the brain glutamate system has been postulated to be involved in the pathophysiology in schizophrenia. To further investigate the association between MYO5B polymorphisms and schizophrenia, we genotyped nine single nucleotide polymorphisms (SNPs) in an independent sample of 1463 individuals with schizophrenia and 1563 healthy control subjects, and detected three SNPs and two haplotype blocks which displayed significant association with schizophrenia. This association was further strengthened by the results of meta-analysis. Our data strongly supported that the MYO5B gene might be associated with schizophrenia in the Chinese Han population and they have implications for understanding the glutamate hypothesis of schizophrenia.

A Two-Stage Association Study Suggests BRAP as a Susceptibility Gene for Schizophrenia

Schizophrenia (SZ) is a neurodevelopmental disorder in which altered immune function typically plays an important role in mediating the effect of environmental insults and regulation of inflammation. The breast cancer suppressor protein associated protein (BRAP) is suggested to exert vital effects in neurodevelopment by modulating the mitogen-activated protein kinase cascade and inflammation signaling. To explore the possible role of BRAP in SZ, we conducted a two-stage study to examine the association of BRAP polymorphisms with SZ in the Han Chinese population. In stage one, we screened SNPs in BRAP from our GWAS data, which detected three associated SNPs, with rs3782886 being the most significant one (P  =  2.31E-6, OR  =  0.67). In stage two, we validated these three SNPs in an independently collected population including 1957 patients and 1509 controls, supporting the association of rs3782886 with SZ (P  =  1.43E-6, OR  =  0.73). Furthermore, cis-eQTL analysis indicates that rs3782886 genotypes are associated with mRNA levels of aldehyde dehydrogenase 2 family (ALDH2) (P  =  0.0039) and myosin regulatory light chain 2 (MYL2) (P < 1.0E-4). Our data suggest that the BRAP gene may confer vulnerability for SZ in Han Chinese population, adding further evidence for the involvement of developmental and/or neuroinflammatory cascades in the illness.