Differential Regulation of Schizophrenia-associated microRNA Gene Function by Variable Number Tandem Repeats (VNTR) polymorphism

Schizophrenia is a complex neuropsychiatric disorder with high heritability (80–90%) and phenotypic variability common to other psychiatric spectrum disorders (International Schizophrenia et al., 2009). In 2011, the Schizophrenia Psychiatric GWAS Consortium (PGC) published the largest schizophrenia GWAS to date, with an initial discovery sample of 21 856 individuals and a replication sample of 29 839 individuals (Schizophrenia Psychiatric Genome-Wide Association Study, 2011). Of the seven genome-wide significant hits, the strongest finding was in an intron of the primary transcript for human microRNA 137 (miR-137). Since then, multiple studies have identified associations between miR-137 and brain function in patients with schizophrenia strongly supporting miR-137 involvement in the disease etiology (Green et al., 2013; Whalley et al., 2012). MicroRNAs (miRNAs), small non-coding RNA molecules predominantly negatively regulating gene expression, are involved in normal neurodevelopmental and physiological processes, and are implicated in the pathogenesis of several neuropsychiatric disorders (Amiel et al., 2012). MiR-137 is a brain-enriched miRNA shown to modulate proliferation and differentiation of adult neurons and affect neuronal maturation, dendritic spine development and morphogenesis (Schizophrenia Psychiatric Genome-Wide Association Study, 2011). In addition to miR-137, four of the seven genome-wide significant genes from PGC were predicted and experimentally validated as miR-137 gene targets (Kwon et al., 2011). These results warrant further investigation into a potential mechanism of susceptibility to schizophrenia mediated by dysregulation of miR-137 and its downstream targets. Variable Nucleotide Tandem Repeats (VNTRs) have been previously associated with neuropsychiatric disorders such as attention deficit hyperactivity disorder, addiction, and schizophrenia (McGeary, 2009; Mill et al., 2005; Prata et al., 2009); however, the molecular effects of VNTRs on disease development have been largely understudied and their role in the pathogenesis of neuropsychiatric disorders remains unclear. A 15bp VNTR at the 5’ end of the miR-137 precursor was first discovered in a melanoma cell line in 2008 (Bemis et al., 2008), and was identified in our postmortem brain sample of 54 schizophrenia cases and 54 controls donated by the Mount Sinai School of Medicine using Sanger sequencing. In this sample the VNTR consisted of 9 alleles (3–10 and 13 repeats) with varying frequency; the major VNTR allele is 3 (72% frequency), followed by 4 (9%), then 9 (6%) and all others were under 5%. We hypothesize that this VNTR is involved in disease susceptibility by affecting the expression of a schizophrenia-associated miRNA and its schizophrenia-associated downstream targets and we tested this hypothesis using in vitro approaches. Due to low allele frequency of the VNTR alleles we were not able to reliably estimate their effects on miR-137 expression in the postmortem sample (data not shown). Therefore, we attempted to assess the effect of the VNTR alleles individually on miR-137 expression in HEK293 cells. These cells were selected based on their low endogenous miR-137 expression, high transfection efficiency and gene expression profile similar to that of neurons (Shaw et al., 2002). Cells were transfected with expression vectors carrying each VNTR allele (3–13 repeats) in the precursor of miR-137 and cells carrying the 9 and 13 VNTR alleles exhibited significantly reduced miR-137 expression as measured by qPCR (Figure 1a; Bonferroni-corrected p<0.038 and p<0.025, respectively). Figure 1 (a) MiRNA expression was assessed by qPCR and relative expression values calculated using 2−ΔΔCt. Data were analyzed using the Kruskal-Wallis test, followed by Dunn’s post-hoc test correcting for multiple testing (GraphPad ... To further clarify the mechanism by which the VNTR alleles impact miR-137 function, we performed a series of luciferase assays to measure the strength of interaction between miR-137 and its previously validated gene targets: TCF4, CACNA1C, WBP1L (formerly: C10orf26), CSMD1 and ZNF804A (Kim et al., 2012; Schizophrenia Psychiatric Genome-Wide Association Study, 2011). Considering that only 9 and 13 VNTR alleles exhibited a significant impact on miR-137 expression, we assessed those alleles’ effects on miR-137 function compared to the major VNTR allele, 3 (as the control). Cells transfected with 9 and 13 VNTR alleles showed elevated luciferase activity compared to 3 VNTR allele and these effects were not observed with the mutated target sequences, confirming that these interactions are miRNA- and gene target-specific (Figures 1b, 1c and 1d). In conclusion, we report a novel mechanism by which a schizophrenia-implicated miRNA’s function is controlled by a VNTR that modulates miR-137 expression and its schizophrenia-associated target genes. Given that the VNTR is present within the miR-137 precursor sequence we hypothesized that 9 and 13 VNTR alleles affect mature miR-137 maturation, expression and function by altering the precursor’s secondary and tertiary structures. Using three dimensional (3D) RNA folding algorithms we observed differences in the tertiary structures caused by the 9 and 13 VNTRs, which we hypothesize is disrupting the mature miRNA’s processing from its precursor sequence (Figure 1e). Although additional studies are needed to further clarify the mechanisms underlying the VNTR’s effects on miR-137 function, our results offer intriguing insights into a potentially novel mechanism of susceptibility to schizophrenia that aligns major GWAS findings within a unified genetic and molecular mechanism in order to provide greater understanding of otherwise disparate findings.