Common genetic risk factors in ASD and ADHD co-occurring families.

Autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) are two major neurodevelopmental disorders that frequently co-occur. However, the genetic mechanism of the co-occurrence remains unclear. The New Jersey Language and Autism Genetics Study (NJLAGS) collected more than 100 families with at least one member affected by ASD. NJLAGS families show a high prevalence of ADHD and provide a good opportunity to study shared genetic risk factors for ASD and ADHD. The linkage study of the NJLAGS families revealed regions on chromosomes 12 and 17 that are significantly associated with ADHD. Using whole-genome sequencing data on 272 samples from 73 NJLAGS families, we identified potential risk genes for ASD and ADHD. Within the linkage regions, we identified 36 genes that are associated with ADHD using a pedigree-based gene prioritization approach. KDM6B (Lysine Demethylase 6B) is the highest-ranking gene, which is a known risk gene for neurodevelopmental disorders, including ASD and ADHD. At the whole-genome level, we identified 207 candidate genes from the analysis of both small variants and structure variants, including both known and novel genes. Using enrichment and protein-protein interaction network analyses, we identified gene ontology terms and pathways enriched for ASD and ADHD candidate genes, such as cilia function and cation channel activity. Candidate genes and pathways identified in our study improve the understanding of the genetic etiology of ASD and ADHD and will lead to new diagnostic or therapeutic interventions for ASD and ADHD in the future.

Structural Variations Contribute to the Genetic Etiology of Autism Spectrum Disorder and Language Impairments.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by restrictive interests and/or repetitive behaviors and deficits in social interaction and communication. ASD is a multifactorial disease with a complex polygenic genetic architecture. Its genetic contributing factors are not yet fully understood, especially large structural variations (SVs). In this study, we aimed to assess the contribution of SVs, including copy number variants (CNVs), insertions, deletions, duplications, and mobile element insertions, to ASD and related language impairments in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Within the cohort, ~77% of the families contain SVs that followed expected segregation or de novo patterns and passed our filtering criteria. These SVs affected 344 brain-expressed genes and can potentially contribute to the genetic etiology of the disorders. Gene Ontology and protein-protein interaction network analysis suggested several clusters of genes in different functional categories, such as neuronal development and histone modification machinery. Genes and biological processes identified in this study contribute to the understanding of ASD and related neurodevelopment disorders.

Characterization hiPSC-derived neural progenitor cells and neurons to investigate the role of NOS1AP isoforms in human neuron dendritogenesis.

Abnormal dendritic arbor development has been implicated in a number of neurodevelopmental disorders, such as autism and Rett syndrome, and the neuropsychiatric disorder schizophrenia. Postmortem brain samples from subjects with schizophrenia show elevated levels of NOS1AP in the dorsolateral prefrontal cortex, a region of the brain associated with cognitive function. We previously reported that the long isoform of NOS1AP (NOS1AP-L), but not the short isoform (NOS1AP-S), negatively regulates dendrite branching in rat hippocampal neurons. To investigate the role that NOS1AP isoforms play in human dendritic arbor development, we adapted methods to generate human neural progenitor cells and neurons using induced pluripotent stem cell (iPSC) technology. We found that increased protein levels of either NOS1AP-L or NOS1AP-S decrease dendrite branching in human neurons at the developmental time point when primary and secondary branching actively occurs. Next, we tested whether pharmacological agents can decrease the expression of NOS1AP isoforms. Treatment of human iPSC-derived neurons with d-serine, but not clozapine, haloperidol, fluphenazine, or GLYX-13, results in a reduction in endogenous NOS1AP-L, but not NOS1AP-S, protein expression; however, d-serine treatment does not reverse decreases in dendrite number mediated by overexpression of NOS1AP isoforms. In summary, we demonstrate how an in vitro model of human neuronal development can help in understanding the etiology of schizophrenia and can also be used as a platform to screen drugs for patients.

d-Serine administration affects nitric oxide synthase 1 adaptor protein and DISC1 expression in sex-specific manner.

Antipsychotic medications are inefficient at treating symptoms of schizophrenia (SCZ), and N-methyl d-aspartate receptor (NMDAR) agonists are potential therapeutic alternatives. As such, these agonists may act on different pathways and proteins altered in the brains of patients with SCZ than do antipsychotic medications. Here, we investigate the effects of administration of the antipsychotic haloperidol and NMDAR agonist d-serine on function and expression of three proteins that play significant roles in SCZ: nitric oxide synthase 1 adaptor protein (NOS1AP), dopamine D2 (D2) receptor, and disrupted in schizophrenia 1 (DISC1). We administered haloperidol or d-serine to male and female Sprague Dawley rats via intraperitoneal injection for 12 days and subsequently examined cortical expression of NOS1AP, D2 receptor, and DISC1. We found sex-specific effects of haloperidol and d-serine treatment on the expression of these proteins. Haloperidol significantly reduced expression of D2 receptor in male, but not female, rats. Conversely, d-serine reduced expression of NOS1AP in male rats and did not affect D2 receptor expression. d-serine treatment also reduced expression of DISC1 in male rats and increased DISC1 expression in female rats. As NOS1AP is overexpressed in the cortex of patients with SCZ and negatively regulates NMDAR signaling, we subsequently examined whether treatment with antipsychotics or NMDAR agonists can reverse the detrimental effects of NOS1AP overexpression in vitro as previously reported by our group. NOS1AP overexpression promotes reduced dendrite branching in vitro, and as such, we treated cortical neurons overexpressing NOS1AP with different antipsychotics (haloperidol, clozapine, fluphenazine) or d-serine for 24 h and determined the effects of these drugs on NOS1AP expression and dendrite branching. While antipsychotics did not affect NOS1AP protein expression or dendrite branching in vitro, d-serine reduced NOS1AP expression and rescued NOS1AP-mediated reductions in dendrite branching. Taken together, our data suggest that d-serine influences the function and expression of NOS1AP, D2 receptor, and DISC1 in a sex-specific manner and reverses the effects of NOS1AP overexpression on dendrite morphology.

Behavioral and Molecular Genetics of Reading-Related AM and FM Detection Thresholds.

Auditory detection thresholds for certain frequencies of both amplitude modulated (AM) and frequency modulated (FM) dynamic auditory stimuli are associated with reading in typically developing and dyslexic readers. We present the first behavioral and molecular genetic characterization of these two auditory traits. Two extant extended family datasets were given reading tasks and psychoacoustic tasks to determine FM 2 Hz and AM 20 Hz sensitivity thresholds. Univariate heritabilities were significant for both AM (h 2  = 0.20) and FM (h 2  = 0.29). Bayesian posterior probability of linkage (PPL) analysis found loci for AM (12q, PPL = 81 %) and FM (10p, PPL = 32 %; 20q, PPL = 65 %). Bivariate heritability analyses revealed that FM is genetically correlated with reading, while AM was not. Bivariate PPL analysis indicates that FM loci (10p, 20q) are not also associated with reading.

Rare exonic deletions implicate the synaptic organizer Gephyrin (GPHN) in risk for autism, schizophrenia and seizures.

The GPHN gene codes for gephyrin, a key scaffolding protein in the neuronal postsynaptic membrane, responsible for the clustering and localization of glycine and GABA receptors at inhibitory synapses. Gephyrin has well-established functional links with several synaptic proteins that have been implicated in genetic risk for neurodevelopmental disorders such as autism spectrum disorder (ASD), schizophrenia and epilepsy including the neuroligins (NLGN2, NLGN4), the neurexins (NRXN1, NRXN2, NRXN3) and collybistin (ARHGEF9). Moreover, temporal lobe epilepsy has been linked to abnormally spliced GPHN mRNA lacking exons encoding the G-domain of the gephyrin protein, potentially arising due to cellular stress associated with epileptogenesis such as temperature and alkalosis. Here, we present clinical and genomic characterization of six unrelated subjects, with a range of neurodevelopmental diagnoses including ASD, schizophrenia or seizures, who possess rare de novo or inherited hemizygous microdeletions overlapping exons of GPHN at chromosome 14q23.3. The region of common overlap across the deletions encompasses exons 3-5, corresponding to the G-domain of the gephyrin protein. These findings, together with previous reports of homozygous GPHN mutations in connection with autosomal recessive molybdenum cofactor deficiency, will aid in clinical genetic interpretation of the GPHN mutation spectrum. Our data also add to the accumulating evidence implicating neuronal synaptic gene products as key molecular factors underlying the etiologies of a diverse range of neurodevelopmental conditions.

The value of regenotyping older linkage data sets with denser marker panels.

OBJECTIVES: Linkage analysis can help determine regions of interest in whole-genome sequence studies. However, many linkage studies rely on older microsatellite (MSAT) panels. We set out to determine whether results would change if we regenotyped families using a dense map of SNPs. METHODS: We selected 47 Hispanic-American families from the NIMH Repository and Genomics Resource (NRGR) schizophrenia data repository. We regenotyped all individuals with DNA available from the NRGR on the Affymetrix Lat Array. After optimizing SNP selection for inclusion on the linkage map, we compared information content (IC) and linkage results using MSAT, SNP and MSAT+SNP maps. RESULTS: As expected, SNP provided a higher average IC (0.78, SD 0.03) than MSAT (0.51, SD 0.10) in a direct 'apples-to-apples' comparison using only individuals genotyped on both platforms; while MSAT+SNP provided only a slightly higher IC (0.82, SD 0.03). However, when utilizing all available individuals, including those who had genotypes available on only one platform, the IC was substantially increased using MSAT+SNP (0.76, SD 0.05) compared to SNP (0.61, SD 0.02). Linkage results changed appreciably between MSAT and MSAT+SNP in terms of magnitude, rank ordering and localization of peaks. CONCLUSIONS: Regenotyping older family data can substantially alter the conclusions of linkage analyses.

The Rosetta Phenotype Harmonization Method Facilitates Finding a Relationship Quantitative Trait Locus for a Complex Cognitive Trait.

Genetics researchers increasingly combine data across many sources to increase power and to conduct analyses that cross multiple individual studies. However, there is often a lack of alignment on outcome measures when the same constructs are examined across studies. This inhibits comparison across individual studies and may impact the findings from meta-analysis. Using a well-characterized genotypic (brain-derived neurotrophic factor: BDNF) and phenotypic constructs (working memory and reading comprehension), we employ an approach called Rosetta, which allows for the simultaneous examination of primary studies that employ related but incompletely overlapping data. We examined four studies of BDNF, working memory, and reading comprehension with a combined sample size of 1711 participants. Although the correlation between working memory and reading comprehension over all participants was high, as expected (ρ = 0.45), the correlation between working memory and reading comprehension was attenuated in the BDNF Met/Met genotype group (ρ = 0.18, n.s.) but not in the Val/Val (ρ = 0.44) or Val/Met (ρ = 0.41) groups. These findings indicate that Met/Met carriers may be a unique and robustly defined subgroup in terms of memory and reading comprehension. This study demonstrates the utility of the Rosetta method when examining complex phenotypes across multiple studies, including psychiatric genetic studies, as shown here, and also for the mega-analysis of cohorts generally.

MicroRNA and MicroRNA-Target Variants Associated with Autism Spectrum Disorder and Related Disorders.

Autism spectrum disorder (ASD) is a childhood neurodevelopmental disorder with a complex and heterogeneous genetic etiology. MicroRNA (miRNA), a class of small non-coding RNAs, could regulate ASD risk genes post-transcriptionally and affect broad molecular pathways related to ASD and associated disorders. Using whole-genome sequencing, we analyzed 272 samples in 73 families in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Families with at least one ASD patient were recruited and were further assessed for language impairment, reading impairment, and other associated phenotypes. A total of 5104 miRNA variants and 1,181,148 3' untranslated region (3' UTR) variants were identified in the dataset. After applying several filtering criteria, including population allele frequency, brain expression, miRNA functional regions, and inheritance patterns, we identified high-confidence variants in five brain-expressed miRNAs (targeting 326 genes) and 3' UTR miRNA target regions of 152 genes. Some genes, such as SCP2 and UCGC, were identified in multiple families. Using Gene Ontology overrepresentation analysis and protein-protein interaction network analysis, we identified clusters of genes and pathways that are important for neurodevelopment. The miRNAs and miRNA target genes identified in this study are potentially involved in neurodevelopmental disorders and should be considered for further functional studies.

Autism NPCs from both idiopathic and CNV 16p11.2 deletion patients exhibit dysregulation of proliferation and mitogenic responses.

Neural precursor cell (NPC) dysfunction has been consistently implicated in autism. Induced pluripotent stem cell (iPSC)-derived NPCs from two autism groups (three idiopathic [I-ASD] and two 16p11.2 deletion [16pDel]) were used to investigate if proliferation is commonly disrupted. All five individuals display defects, with all three macrocephalic individuals (two 16pDel, one I-ASD) exhibiting hyperproliferation and the other two I-ASD subjects displaying hypoproliferation. NPCs were challenged with bFGF, and all hyperproliferative NPCs displayed blunted responses, while responses were increased in hypoproliferative cells. mRNA expression studies suggest that different pathways can result in similar proliferation phenotypes. Since 16pDel deletes MAPK3, P-ERK was measured. P-ERK is decreased in hyperproliferative but increased in hypoproliferative NPCs. While these P-ERK changes are not responsible for the phenotypes, P-ERK and bFGF response are inversely correlated with the defects. Finally, we analyzed iPSCs and discovered that 16pDel displays hyperproliferation, while idiopathic iPSCs were normal. These data suggest that NPC proliferation defects are common in ASD.

Genetic covariation underlying reading, language and related measures in a sample selected for specific language impairment.

Specific language impairment is a developmental language disorder characterized by failure to develop language normally in the absence of a specific cause. Previous twin studies have documented the heritability of reading and language measures as well as the genetic correlation between those measures. This paper presents results from an alternative to the classical twin designs by estimating heritability from extended pedigrees. These pedigrees were previously studied as part of series of molecular genetic studies of specific language impairment where the strongest genetic findings were with reading phenotypes rather than language despite selecting pedigrees based on language impairments. To explore the relationship between reading and language in these pedigrees, variance components estimates of heritability of reading and language measures were conducted showing general agreement with the twin literature, as were genetics correlations between reading and language. Phonological short-term memory, phonological awareness and auditory processing were evaluated as candidate mediators of the reading-language genetic correlations. Only phonological awareness showed significant genetic correlations with all reading measures and several language measures while phonological short-term memory and auditory processing did not.

Increasing genotype-phenotype model determinism: application to bivariate reading/language traits and epistatic interactions in language-impaired families.

While advances in network and pathway analysis have flourished in the era of genome-wide association analysis, understanding the genetic mechanism of individual loci on phenotypes is still readily accomplished using genetic modeling approaches. Here, we demonstrate two novel genotype-phenotype models implemented in a flexible genetic modeling platform. The examples come from analysis of families with specific language impairment (SLI), a failure to develop normal language without explanatory factors such as low IQ or inadequate environment. In previous genome-wide studies, we observed strong evidence for linkage to 13q21 with a reading phenotype in language-impaired families. First, we elucidate the genetic architecture of reading impairment and quantitative language variation in our samples using a bivariate analysis of reading impairment in affected individuals jointly with language quantitative phenotypes in unaffected individuals. This analysis largely recapitulates the baseline analysis using the categorical trait data (posterior probability of linkage (PPL) = 80%), indicating that our reading impairment phenotype captured poor readers who also have low language ability. Second, we performed epistasis analysis using a functional coding variant in the brain-derived neurotrophic factor (BDNF) gene previously associated with reduced performance on working memory tasks. Modeling epistasis doubled the evidence on 13q21 and raised the PPL to 99.9%, indicating that BDNF and 13q21 susceptibility alleles are jointly part of the genetic architecture of SLI. These analyses provide possible mechanistic insights for further cognitive neuroscience studies based on the models developed herein.

NOS1AP regulates dendrite patterning of hippocampal neurons through a carboxypeptidase E-mediated pathway.

During neuronal development, neurons form elaborate dendritic arbors that receive signals from axons. Additional studies are needed to elucidate the factors regulating the establishment of dendritic patterns. Our work explored possible roles played by nitric oxide synthase 1 adaptor protein (NOS1AP; also known as C-terminal PDZ ligand of neuronal nitric oxide synthase or CAPON) in dendritic patterning of cultured hippocampal neurons. Here we report that the long isoform of NOS1AP (NOS1AP-L) plays a novel role in regulating dendrite outgrowth and branching. NOS1AP-L decreases dendrite number when overexpressed at any interval between day in vitro (DIV) 0 and DIV 12, and knockdown of NOS1AP-L results in increased dendrite number. In contrast, the short isoform of NOS1AP (NOS1AP-S) decreases dendrite number only when overexpressed during DIV 5-7. Using mutants of NOS1AP-L, we show that neither the PDZ-binding domain nor the PTB domain is necessary for the effects of NOS1AP-L. We have functionally narrowed the region of NOS1AP-L that mediates this effect to the middle amino acids 181-307, a region that is not present in NOS1AP-S. Furthermore, we performed a yeast two-hybrid screen and identified carboxypeptidase E (CPE) as a binding partner for the middle region of NOS1AP-L. Biochemical and cellular studies reveal that CPE mediates the effects of NOS1AP on dendrite morphology. Together, our results suggest that NOS1AP-L plays an important role in the initiation, outgrowth, and maintenance of dendrites during development.

Evidence for a role of the NOS1AP (CAPON) gene in schizophrenia and its clinical dimensions: an association study in a South American population isolate.

BACKGROUND/AIMS: Recent studies have implicated a region on chromosome 1q21-23, including the NOS1AP gene, in susceptibility to schizophrenia. However, replication studies have been inconsistent, a fact that could partly relate to the marked psychopathological heterogeneity of schizophrenia. The aim of this study is to evaluate association of polymorphisms in the NOS1AP gene region to schizophrenia, in patients from a South American population isolate, and to assess if these variants are associated with specific clinical dimensions of the disorder. METHODS: We genotyped 24 densely spaced SNPs in the NOS1AP gene region in a schizophrenia trio sample. The transmission disequilibrium test (TDT) was applied to single marker and haplotype data. Association to clinical dimensions (identified by factor analysis) was evaluated using a quantitative transmission disequilibrium test (QTDT). RESULTS: We found significant association between eight SNPs in the NOS1AP gene region to schizophrenia (minimum p value = 0.004). The QTDT analysis of clinical dimensions revealed an association to a dimension consisting mainly of negative symptoms (minimum p value 0.001). CONCLUSIONS: Our findings are consistent with a role for NOS1AP in susceptibility to schizophrenia, especially for the 'negative syndrome' of the disorder.

Childhood maltreatment, serotonin transporter gene, and risk for callous and unemotional traits: A prospective investigation.

Previous studies have reported associations between the serotonin transporter 5-HTTLPR genotype and antisocial and aggressive traits and between child maltreatment and antisocial traits. However, few studies have examined whether 5-HTTLPR moderates the influence of childhood maltreatment on callous and unemotional traits, a hallmark of psychopathy. Using a prospective cohort design, children with documented cases of maltreatment and matched controls were followed up and interviewed in adulthood. DNA was extracted from blood and saliva (N = 414) and callous-unemotional (CU) traits were assessed. Childhood maltreatment predicted higher CU scores in adulthood, whereas the effect of 5-HTTLPR was not significant. The effect of child maltreatment on CU traits did not differ by genetic risk (high or low activity 5-HTTLPR), whereas controls with the LL genotype had higher CU scores than controls with the SS genotype. Similar results were found for females and White, non-Hispanics, but not for males and Blacks. Variations in 5-HTTLPR did not affect the impact of child maltreatment on CU traits in adulthood. Genetic risk had a stronger effect on adults with lower environmental risk (controls). Having a history of child maltreatment or the LL genotype placed participants at risk for higher levels of callous and unemotional trait scores.

Social (Pragmatic) Communication Disorder: Another name for the Broad Autism Phenotype?

The Diagnostic and Statistical Manual of Mental Disorders' (5th ed.) Social (Pragmatic) Communication Disorder is meant to capture the social elements of communication dysfunction in children who do not meet autism spectrum disorder criteria. It is unclear whether Social (Pragmatic) Communication Disorder captures these elements without overlapping with Autism Spectrum Disorder or the Diagnostic and Statistical Manual of Mental Disorders' (5th ed.) Language Disorder. Standardized behavioral assessments administered during a family genetics study were used to evaluate the social communication impairment and the restricted interests and repetitive behaviors in persons with autism spectrum disorder, language impairment, or neither. Social communication impairment and restricted interests and repetitive behavior were significantly correlated in all family members regardless of affection status. Rates of social communication impairment and restricted interests and repetitive behavior were highest in individuals with autism spectrum disorder. One-third of family members with language impairment presented with at least mild/moderate levels of social communication impairment (36.6%) and restricted interests and repetitive behavior (43.3%). A subset of unaffected members also presented with mild/moderate levels of social communication impairment (parents = 10.1%, siblings 11.6%) and restricted interests and repetitive behavior (parents = 14.0%, siblings = 22.1%). The majority of child family members with mild/moderate levels of social communication impairment had similar restricted interest and repetitive behavior levels reflecting criteria representing the Broad Autism Phenotype. These data suggest that social pragmatic communication disorder does not capture the profiles of children who have both social communication impairment and restricted interests and repetitive behavior but are in need of clinical services.

Within-task variability on standardized language tests predicts autism spectrum disorder: a pilot study of the Response Dispersion Index.

BACKGROUND: Qualitatively atypical language development characterized by non-sequential skill acquisition within a developmental domain, which has been called developmental deviance or difference, is a common characteristic of autism spectrum disorder (ASD). We developed the Response Dispersion Index (RDI), a measure of this phenomenon based on intra-subtest scatter of item responses on standardized psychometric assessments, to assess the within-task variability among individuals with language impairment (LI) and/or ASD. METHODS: Standard clinical assessments of language were administered to 502 individuals from the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Participants were divided into four diagnostic groups: unaffected, ASD-only, LI-only, and ASD + LI. For each language measure, RDI was defined as the product of the total number of test items and the sum of the weight (based on item difficulty) of test items missed. Group differences in RDI were assessed, and the relationship between RDI and ASD diagnosis among individuals with LI was investigated for each language assessment. RESULTS: Although standard scores were unable to distinguish the LI-only and ASD/ASD + LI groups, the ASD/ASD + LI groups had higher RDI scores compared to LI-only group across all measures of expressive, pragmatic, and metalinguistic language. RDI was positively correlated with quantitative ASD traits across all subgroups and was an effective predictor of ASD diagnosis among individuals with LI. CONCLUSIONS: The RDI is an effective quantitative metric of developmental deviance/difference that correlates with ASD traits, supporting previous associations between ASD and non-sequential skill acquisition. The RDI can be adapted to other clinical measures to investigate the degree of difference that is not captured by standard performance summary scores.

Improvements to bead-based oligonucleotide ligation SNP genotyping assays.

We describe a bead-based, multiplexed, oligonucleotide ligation assay (OLA) performed on the Luminex flow cytometer. Differences between this method and those previously reported include the use of far fewer beads and the use of a universal oligonucleotide for signal detection. These innovations serve to significantly reduce the cost of the assay, while maintaining robustness and accuracy. Comparisons are made between the Luminex OLA and both pyrosequencing and direct sequencing. Experiments to assess conversion rates, call rates, and concordance across technical replicates are also presented.

NOS1AP in schizophrenia.

NOS1AP is an attractive candidate gene for schizophrenia susceptibility. Linkage and association studies from multiple samples drawn from different populations indicate that a schizophrenia susceptibility gene is located in the region of chromosome 1 containing NOS1AP. Increased NOS1AP expression is observed in postmortem samples from individuals with schizophrenia. NOS1AP binds to neuronal nitric oxide synthase and synapsin, other candidate genes for schizophrenia, and may disrupt signal transduction through the N-methyl d-aspartate receptor complex, leading to hypofunctioning of that system. In this review, I present the evidence supporting NOS1AP as a schizophrenia susceptibility gene, with a focus on explaining the strengths and weaknesses of the evidence obtained from each type of study that has been conducted.

Validation of a microRNA target site polymorphism in H3F3B that is potentially associated with a broad schizophrenia phenotype.

Despite much progress, few genetic findings for schizophrenia have been assessed by functional validation experiments at the molecular level. We previously reported evidence for genetic linkage of broadly defined schizophrenia to chromosome 17q25 in a sample of 24 multiplex families. 2,002 SNPs under this linkage peak were analyzed for evidence of linkage disequilibrium using the posterior probability of linkage (PPL) framework. SNP rs1060120 produced the strongest evidence for association, with a PPLD|L score of 0.21. This SNP is located within the 3'UTR of the histone gene H3F3B and colocalizes with potential gene target miR-616. A custom miRNA target prediction program predicted that the binding of miR-616 to H3F3B transcripts would be altered by the allelic variants of rs1060120. We used dual luciferase assays to experimentally validate this interaction. The rs1060120 A allele significantly reduced luciferase expression, indicating a stronger interaction with miR-616 than the G allele (p = 0.000412). These results provide functional validation that this SNP could alter schizophrenia epigenetic mechanisms thereby contributing to schizophrenia-related disease risk.