Vladimir I. Vladimirov

Replication of association between schizophrenia and ZNF804A in the Irish Case-Control Study of Schizophrenia sample

A recent genome-wide association study reported association between schizophrenia and the ZNF804A gene on chromosome 2q32.1. We attempted to replicate these findings in our Irish Case–Control Study of Schizophrenia (ICCSS) sample (N=1021 cases, 626 controls). Following consultation with the original investigators, we genotyped three of the most promising single-nucleotide polymorphisms (SNPs) from the Cardiff study. We replicate association with rs1344706 (trend test one-tailed P=0.0113 with the previously associated A allele) in ZNF804A. We detect no evidence of association with rs6490121 in NOS1 (one-tailed P=0.21), and only a trend with rs9922369 in RGRIP1L (one-tailed P=0.0515). On the basis of these results, we completed genotyping of 11 additional linkage disequilibrium-tagging SNPs in ZNF804A. Of 12 SNPs genotyped, 11 pass quality control criteria and 4 are nominally associated, with our most significant evidence of association at rs7597593 (P=0.0013) followed by rs1344706. We observe no evidence of differential association in ZNF804A on the basis of family history or sex of case. The associated SNP rs1344706 lies in ∼30 bp of conserved mammalian sequence, and the associated A allele is predicted to maintain binding sites for the brain-expressed transcription factors MYT1l and POU3F1/OCT-6. In controls, expression is significantly increased from the A allele of rs1344706 compared with the C allele. Expression is increased in schizophrenic cases compared with controls, but this difference does not achieve statistical significance. This study replicates the original reported association of ZNF804A with schizophrenia and suggests that there is a consistent link between the A allele of rs1344706, increased expression of ZNF804A and risk for schizophrenia.

Genome-wide pharmacogenomic analysis of response to treatment with antipsychotics

Schizophrenia is an often devastating neuropsychiatric illness. Understanding the genetic variation affecting response to antipsychotics is important to develop novel diagnostic tests to match individual schizophrenia patients to the most effective and safe medication. In this study, we use a genome-wide approach to detect genetic variation underlying individual differences in response to treatment with the antipsychotics olanzapine, quetiapine, risperidone, ziprasidone and perphenazine. Our sample consisted of 738 subjects with DSM-IV schizophrenia who took part in the Clinical Antipsychotic Trials of Intervention Effectiveness. Subjects were genotyped using the Affymetrix 500 K genotyping platform plus a custom 164 K chip to improve genome-wide coverage. Treatment outcome was measured using the Positive and Negative Syndrome Scale. Our criterion for genome-wide significance was a prespecified threshold that ensures that, on an average, only 10% of the significant findings are false discoveries. The top statistical result reached significance at our prespecified threshold and involved a single-nucleotide polymorphism (SNP) in an intergenic region on chromosome 4p15. In addition, SNPs in Ankyrin Repeat and Sterile Alpha Motif Domain-Containing Protein 1B (ANKS1B) and in the Contactin-Associated Protein-Like 5 gene (CNTNAP5), which mediated the effects of olanzapine and risperidone on Negative symptoms, were very close to our threshold for declaring significance. The most significant SNP in CNTNAP5 is nonsynonymous, giving rise to an amino-acid substitution. In addition to highlighting our top results, we provide all P-values for download as a resource for investigators with the requisite samples to carry out replication. This study demonstrates the potential of genome-wide association studies to discover novel genes that mediate the effects of antipsychotics, which could eventually help to tailor drug treatment to schizophrenic patients.

A meta-analysis of gene expression quantitative trait loci in brain

Current catalogs of brain expression quantitative trait loci (eQTL) are incomplete and the findings do not replicate well across studies. All existing cortical eQTL studies are small and emphasize the need for a meta-analysis. We performed a meta-analysis of 424 brain samples across five studies to identify regulatory variants influencing gene expression in human cortex. We identified 3584 genes in autosomes and chromosome X with false discovery rate q<0.05 whose expression was significantly associated with DNA sequence variation. Consistent with previous eQTL studies, local regulatory variants tended to occur symmetrically around transcription start sites and the effect was more evident in studies with large sample sizes. In contrast to random SNPs, we observed that significant eQTLs were more likely to be near 5’-untranslated regions and intersect with regulatory features. Permutation-based enrichment analysis revealed that SNPs associated with schizophrenia and bipolar disorder were enriched among brain eQTLs. Genes with significant eQTL evidence were also strongly associated with diseases from OMIM (Online Mendelian Inheritance in Man) and the NHGRI (National Human Genome Research Institute) genome-wide association study catalog. Surprisingly, we found that a large proportion (28%) of ~1000 autosomal genes encoding proteins needed for mitochondrial structure or function were eQTLs (enrichment P-value=1.3 × 10−9), suggesting a potential role for common genetic variation influencing the robustness of energy supply in brain and a possible role in the etiology of some psychiatric disorders. These systematically generated eQTL information should be a valuable resource in determining the functional mechanisms of brain gene expression and the underlying biology of associations with psychiatric disorders.

The AVPR1A Gene and Substance Use Disorders: Association, Replication, and Functional Evidence

BACKGROUND The liability to addiction has been shown to be highly genetically correlated across drug classes, suggesting nondrug-specific mechanisms. METHODS In 757 subjects, we performed association analysis between 1536 single nucleotide polymorphisms (SNPs) in 106 candidate genes and a drug use disorder diagnosis (DUD). RESULTS Associations (p ≤ .0008) were detected with three SNPs in the arginine vasopressin 1A receptor gene, AVPR1A, with a gene-wise p value of 3 × 10(-5). Bioinformatic evidence points to a role for rs11174811 (microRNA binding site disruption) in AVPR1A function. Based on literature implicating AVPR1A in social bonding, we tested spousal satisfaction as a mediator of the association of rs11174811 with the DUD. Spousal satisfaction was significantly associated with DUD in males (p < .0001). The functional AVPR1A SNP, rs11174811, was associated with spousal satisfaction in males (p = .007). Spousal satisfaction was a significant mediator of the relationship between rs11174811 and DUD. We also present replication of the association in males between rs11174811 and substance use in one clinically ascertained (n = 1399) and one epidemiologic sample (n = 2231). The direction of the association is consistent across the clinically-ascertained samples but reversed in the epidemiologic sample. Lastly, we found a significant impact of rs11174811 genotype on AVPR1A expression in a postmortem brain sample. CONCLUSIONS The findings of this study call for expansion of research into the role of the arginine vasopressin and other neuropeptide system variation in DUD liability.

Detecting miRNAs in deep-sequencing data: a software performance comparison and evaluation

Deep sequencing has become a popular tool for novel miRNA detection but its data must be viewed carefully as the state of the field is still undeveloped. Using three different programs, miRDeep (v1, 2), miRanalyzer and DSAP, we have analyzed seven data sets (six biological and one simulated) to provide a critical evaluation of the programs performance. We selected these software based on their popularity and overall approach toward the detection of novel and known miRNAs using deep-sequencing data. The program comparisons suggest that, despite differing stringency levels they all identify a similar set of known and novel predictions. Comparisons between the first and second version of miRDeep suggest that the stringency level of each of these programs may, in fact, be a result of the algorithm used to map the reads to the target. Different stringency levels are likely to affect the number of possible novel candidates for functional verification, causing undue strain on resources and time. With that in mind, we propose that an intersection across multiple programs be taken, especially if considering novel candidates that will be targeted for additional analysis. Using this approach, we identify and performed initial validation of 12 novel predictions in our in-house data with real-time PCR, six of which have been previously unreported.

Association of CRHR1 variants and posttraumatic stress symptoms in hurricane exposed adults

Posttraumatic stress disorder (PTSD) is a moderately heritable anxiety disorder that may develop after exposure to trauma. However, only few genetic variants that relate to PTSD have been studied. This study examined the relationship between 12 single nucleotide polymorphisms (SNPs) in the corticotropin-releasing hormone receptor 1 gene (CRHR1) and post-disaster PTSD symptoms and diagnosis in adults exposed to 2004 Florida hurricanes. CRHR1 regulates the hypothalamic-pituitary-adrenal (HPA) axis; dysregulation of the HPA axis is characteristic of stress phenotypes. Final analyses were conducted in the European-American (EA) subsample (n=564) due to population stratification. After correction for multiple testing, rs12938031 and rs4792887 remained associated with post-hurricane PTSD symptoms. Additionally, rs12938031 was associated with post-hurricane diagnosis of PTSD. This study is the first to examine CRHR1 in relation to PTSD in adults, and provides evidence for the importance of CRHR1 variation in the etiology of PTSD. Although results are preliminary and require replication, they justify follow-up efforts to characterize how this gene relates to PTSD.

Metal transcription factor-1 regulation via MREs in the transcribed regions of selenoprotein H and other metal- responsive genes

Selenoprotein H is a redox-sensing DNA binding protein that upregulates genes involved in antioxidant responses. Given the known links between oxidative stress and heavy metals, we investigated the potential for regulation of selenoprotein H by metals. In silico analysis of the selenoprotein H genes from nine species reveals multiple predicted metal response elements (MREs). To validate MRE function, we investigated the effects of zinc or cadmium addition and metal-responsive transcription factor 1 (MTF-1) knockout on selenoprotein H mRNA levels. Chromatin immunoprecipitation was used to directly assess physical binding of the transcription factor to MREs in the human and mouse selenoprotein H genes. The results reported herein show that selenoprotein H is a newly identified target for MTF-1. Further, whereas nearly all prior studies of MREs focused on those located in promoters, we demonstrate binding of MTF-1 to MREs located downstream of the transcription start sites in the human and murine selenoprotein H genes. Finally, we identified MREs in downstream sequences in 15 additional MTF-1 regulated genes lacking promoter MREs, and demonstrated MTF-1 binding in three of these genes. This regulation via sequences downstream of promoters highlights a new direction for identifying previously unrecognized target genes for MTF-1.

Cellular models to study bipolar disorder: a systematic review

BACKGROUND There is an emerging interest in the use of cellular models to study psychiatric disorders. We have systematically reviewed the application of cellular models to understand the biological basis of bipolar disorder (BD). METHOD Published scientific literature in MEDLINE, PsychINFO and SCOPUS databases were identified with the following search strategy: [(Lymphoblastoid OR Lymphoblast OR Fibroblast OR Pluripotent OR Olfactory epithelium OR Olfactory mucosa) AND (Bipolar disorder OR Lithium OR Valproate OR Mania)]. Studies were included if they had used cell cultures derived from BD patients. RESULTS There were 65 articles on lymphoblastoid cell lines, 14 articles on fibroblasts, 4 articles on olfactory neuronal epithelium (ONE) and 2 articles on neurons reprogrammed from induced pluripotent stem cell lines (IPSC). Several parameters have been studied, and the most replicated findings are abnormalities in calcium signaling, endoplasmic reticulum (ER) stress response, mitochondrial oxidative pathway, membrane ion channels, circadian system and apoptosis related genes. These, although present in basal state, seem to be accentuated in the presence of cellular stressors (e.g. oxidative stress--rotenone; ER stress--thapsigargin), and are often reversed with in-vitro lithium. CONCLUSION Cellular modeling has proven useful in BD, and potential pathways, especially in cellular resilience related mechanisms have been identified. These findings show consistency with other study designs (genome-wide association, brain-imaging, and post-mortem brain expression). ONE cells and IPSC reprogrammed neurons represent the next generation of cell models in BD. Future studies should focus on family-based study designs and combine cell models with deep sequencing and genetic manipulations.

Integrating mRNA and miRNA Weighted Gene Co-Expression Networks with eQTLs in the Nucleus Accumbens of Subjects with Alcohol Dependence

Alcohol consumption is known to lead to gene expression changes in the brain. After performing weighted gene co-expression network analyses (WGCNA) on genome-wide mRNA and microRNA (miRNA) expression in Nucleus Accumbens (NAc) of subjects with alcohol dependence (AD; N = 18) and of matched controls (N = 18), six mRNA and three miRNA modules significantly correlated with AD were identified (Bonferoni-adj. p≤ 0.05). Cell-type-specific transcriptome analyses revealed two of the mRNA modules to be enriched for neuronal specific marker genes and downregulated in AD, whereas the remaining four mRNA modules were enriched for astrocyte and microglial specific marker genes and upregulated in AD. Gene set enrichment analysis demonstrated that neuronal specific modules were enriched for genes involved in oxidative phosphorylation, mitochondrial dysfunction and MAPK signaling. Glial-specific modules were predominantly enriched for genes involved in processes related to immune functions, i.e. cytokine signaling (all adj. p≤ 0.05). In mRNA and miRNA modules, 461 and 25 candidate hub genes were identified, respectively. In contrast to the expected biological functions of miRNAs, correlation analyses between mRNA and miRNA hub genes revealed a higher number of positive than negative correlations (χ2 test p≤ 0.0001). Integration of hub gene expression with genome-wide genotypic data resulted in 591 mRNA cis-eQTLs and 62 miRNA cis-eQTLs. mRNA cis-eQTLs were significantly enriched for AD diagnosis and AD symptom counts (adj. p = 0.014 and p = 0.024, respectively) in AD GWAS signals in a large, independent genetic sample from the Collaborative Study on Genetics of Alcohol (COGA). In conclusion, our study identified putative gene network hubs coordinating mRNA and miRNA co-expression changes in the NAc of AD subjects, and our genetic (cis-eQTL) analysis provides novel insights into the etiological mechanisms of AD.

A region of 35 kb containing the trace amine associate receptor 6 (TAAR6) gene is associated with schizophrenia in the Irish study of high-density schizophrenia families.

The TAAR6 gene has been previously associated with schizophrenia in 192 pedigrees of European and African ancestry. To replicate these findings we performed an association study of TAAR6 in 265 pedigrees of the Irish Study of High-Density Schizophrenia Families (ISHDSF). Of the 24 genotyped single-nucleotide polymorphisms only rs12189813 and rs9389011 provided single-marker evidence for association (0.0094⩽P⩽0.03). Two-marker haplotypes (rs7772821 and rs12189813; 0.0071⩽P⩽0.0023) and four-marker haplotypes (rs8192622, rs7772821, rs12189813 and rs9389011; 0.0047⩽P⩽0.018) gave strongest evidence for association. The associated high-risk (HR) haplotype in the ISHDSF is defined by the major alleles at rs7772821 and rs12189813 (0.00097⩽P⩽0.023). The associated HR remains positive in a case only test of association by Operational Criteria score analysis in which significant association was observed only with the highest threshold for delusions (P<0.009). When analysis was restricted to affected individuals under the core schizophrenia (D2) diagnosis, the observed associations for HR were most significant in the highest threshold for delusions (P<0.004) and hallucinations (P<0.0004), supporting the family-based association with schizophrenia.