MicroRNAs targeting Nicastrin regulate Aβ production and are affected by target site polymorphisms

Introduction

The small, non-coding miRNAs function as important regulators of protein expression (Ambros, 2004). They act by imperfect base-pairing to target messenger RNAs, leading to translational repression, degradation, or both (Pillai, 2005). Important for this action is the seed sequence, located at positions 2–8 within the miRNA sequence. It is now well documented that specific miRNAs are altered in a number of neurodegenerative disorders including Alzheimer’s disease (AD; reviewed in Delay et al., 2012). These observations evoke the hypothesis that miRNA dysregulation could directly be implicated in the disease processes, either detrimental or as an attempt of neurons to survive. Moreover, increasing evidence suggests that single-nucleotide polymorphisms (SNP) within miRNA target sites might contribute to disease risk. For instance, the Parkinson’s disease-associated rs12720208 SNP is located in the miR-433 target site of Fibroblast growth factor 20, and affects miRNA-mediated repression (Wang et al., 2008). The rs5848 SNP, associated with TDP43-positive frontotemporal dementias, is located in the miR-659 target site of the Progranulin gene, also affecting the miRNA-mediated repression (Rademakers et al., 2008). Previously, we have shown that AD-associated SNPs present in the untranslated region (3′UTR) of APP could directly affect miRNA function (Delay et al., 2011). In theory, loss of miRNA function could contribute significantly to risk for AD by increasing APP expression and therefore Aβ production. Aβ peptides are the major constituents of amyloid (senile) plaques that accumulate in AD brain (Glenner and Wong, 1984).

Another essential player in Aβ production is the γ-secretase complex, composed of four subunits (i.e., Presenilin, Nicastrin, Aph-1, and Pen2; De Strooper, 2003). Recent studies have documented rare Nicastrin (NCSTN) variants associated with AD (Lupton et al., 2011), whereas studies in mice have shown that increased expression of NCSTN is sufficient to promote AD pathology (Goo et al., 2013). In this study, we evaluated whether miRNAs could modulate NCSTN expression, and whether polymiRTS within the 3′UTR of the Nicastrin gene could abrogate miRNA function.

Materials and Methods

Patients

Cases of AD and controls were participants from two large Canadian cohorts: the Canadian Study of Health and Aging (CSHA) and A Canadian Collaborative Cohort of Cognitive Impairment and Related Dementia (ACCORD). Ethical approval of these studies was obtained from the ethics review board in each of the study centers. Details of each cohort have been published elsewhere (Anonymous, 1994; Feldman et al., 2003; Hsiung et al., 2004). The clinical examination was developed in collaboration with the US Consortium to Establish a Registry for Alzheimer’s Disease. AD cases were diagnosed by NINCDS – ADRDA clinical criteria (McKhann et al., 1984), and controls were subjects who remained not cognitively impaired until the end of the study or censorship.

Cell Culture

Human HeLa and HEK293-APPSw cells were cultured in DMEM medium (Invitrogen, Carlsbad, CA, USA) supplemented with 10% heat-inactivated fetal bovine serum. One day before transfection, HeLa cells were plated at 20% confluence in six-well or 24-well plates and 250,000 HEK-APPSw cells were plated into six-well plates. Transfection was performed using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions.

cDNA Constructs

The full-length hNCSTN 3′UTR was extracted by PCR from SH-SY5Y cells using the following primers: forward: 5′-TTTTTCTAGAGGAGGACCCCAGCTTTTC, reverse 5′-AAAAGGATCCCGTGTGGGATAATCTATTTT, and were cloned using XbaI (New England Biolabs) and BamHI (New England Biolabs) enzymes into the pGL3-promotor vector (Promega, USA). Mutagenesis was performed by TOPgene technologies (Montreal, Quebec, Canada) and validated by sequencing.

Luciferase Assay and Protein Analysis

HeLa cells were transfected with 25 or 50 nM pre-miRs (Applied Biosystems, USA), 2.5 ng/cm²; pRL control vector, and 50 ng/cm²; pGL3_HSV TK_3′UTR hNCSTN WT or C460T (rs1043329), T623G (rs113810300) or delCA515–516 (rs141849450) mutant variant plasmids. Twenty-four hours post-transfection, cells were lyzed, and luciferase activity was measured according to the manufacturer’s instructions (Promega, USA). For western blots, cells were lyzed in RIPA buffer [50 mM Tris HCl, 1% NP40, 0.9% NaCl, 0.25% Na-deoxycholate, 1 mM EDTA, 1× proteinase inhibitors (Roche, Basel, Switzerland), 1 mM PMSF, 1 mM Na₃VO₄, and 1 mM NaF], mixed with LDS sample buffer (Invitrogen, Carlsbad, CA, USA) containing 5% β-mercapto-ethanol and boiled at 95°C for 8 min. Crude extracts (10 μg) were immunoblotted with the Aph-1a [clone B80.2, kind gift from B. De Strooper (Nyabi et al., 2003)], Pen2 (Cell Signaling, clone D2G6), Nicastrin (Cell Signaling, clone D38F9), Presenilin 1 (Millipore, clone PS1-loop), and Gapdh (Millipore, clone 6c5) antibodies. Membranes were detected using the ECL detection kit (Millipore, Billerica, MA, USA). Quantifications were performed using the Multi Gauge software (FUJIFILM, Minato-ku, Tokyo, Japan).

Elisa

HEK293-APPSw cells were transfected with 50 nM pre-miRs (Applied Biosystems, USA). Twenty-four hours post-transfection, cell lysates were collected, spun at 1000 × g for 10 min to remove debris, and supernatants were kept on ice. Human (soluble) Aβ1–40 and Aβ1–42 peptides were measured using Aβ40 and Aβ42 Human ELISA Kits, following the manufacturer’s conditions (Life Technologies, Carlsbad, CA, USA).

Statistics

Statistical significance of western blots and luminescence quantifications was determined using one-way ANOVA, two-way ANOVA, or Student’s paired t-test as indicated in the text. Calculations were made using the GraphPad Prism 5 software.

Genotyping

DNA was extracted from peripheral blood (AutogenFlex STAR, Holliston MA, USA), and SNP genotyping was performed with TaqMan assays in an optimized ABI 7300 (Applied Biosystems). The two SNPs chosen in this study were rs113810300 and rs141849450 in human Nicastrin.

Results

We first used the National Center for Biotechnology Information (NCBI) database to generate a list of all known SNPs located within the 3′UTR of human NCSTN (Table 1). In order to identify SNPs that could affect miRNA binding, we used the on-line algorithms TargetScan and microRNA.org (Betel et al., 2008; Garcia et al., 2011). We included both canonical (i.e., perfect seed match) and non-canonical (i.e., containing G:U wobbles or seed mismatches) miRNA binding sites. This analysis resulted in a list of 22 miRNAs (Table 1), which we narrowed down to six (i.e., miR-24, miR-186, miR-340, miR-455, miR-656, and miR-1301) based on our previous expression profiling studies in the human cerebral neocortex (45 raw reads cut-off; Hébert et al., 2013). Whether excluded miRNAs are expressed at higher frequency in other brain regions and/or tissues remains to be determined.

TABLE 1

TABLE 1. Polymorphisms located in or near miRNA target sites located in the 3′UTR of hNCSTN.

To validate our bioinformatic analyses, we cloned the full-length (∼700 bp) 3′UTR of human NCSTN into a luciferase reporter vector (Figure 1A), and co-transfected this construct with candidate miRNAs into human HeLa cells. We observed a significant reduction of luciferase expression upon miR-24, miR-186, and miR-455 expression (Figure 1A) compared to a scrambled miRNA (SCR) negative control. We also confirmed the downregulation of endogenous NCSTN after transfection of candidate miRNAs into HeLa cells (Figure 1B). Notably, only miR-186 and miR-455 decreased both the mature and immature forms of NCSTN (Edbauer et al., 2002). While both miR-340 and miR-656 also affected luciferase and endogenous NCSTN expression, this observation did not reach significance. Overexpression of identified NCSTN targeting miRNAs did not significantly affect endogenous Presenilin 1 (PS1), APH1A, or PEN2 expression (Figure 1C). We thus identified miR-24, miR-186, and miR-455 as endogenous regulators of human NCSTN.

FIGURE 1

FIGURE 1. In vitro functional analysis of NCSTN targeting miRNAs. (A) Schematic representation (not to scale) of the luciferase reporter construct used in this study (upper panel). TK; thymidine kinase promoter, AAAA(n); PolyA site, HeLa cells were transfected with 50 nM pre-miRs (as indicated) as well as a reporter construct containing the untranslated region (3′UTR) of hNCSTN. The cells were lyzed 24 h post-transfection and luciferase signal was measured. Signals were normalization for transfection efficiency and graph represents the luciferase signals compared to the scrambled control (SCR). Statistical significance was assessed by one-way ANOVA with Bonferroni post-hoc test. *p < 0.05. (B) HeLa cells were transfected with 50 nM pre-miRNAs (as indicated). The cells were lyzed 48 h post-transfection and western blotting was performed. Representative (n = 3, in triplicate) western blots are shown. The ratios of the NCSTN/GAPDH signals are presented (both mature and immature NSCTN). Measurements were normalized to the scrambled control (SCR). Statistical significance was assessed by one-way ANOVA with Bonferroni post-hoc test. *p < 0.05, **p < 0.01, ***p < 0.001. (C) HeLa cells were transfected with 50 nM pre-miRNAs (as indicated). The cells were lyzed 48 h post-transfection and western blotting was performed. Representative (n = 3, in triplicate) western blots and quantifications are shown. Statistical significance was assessed by one-way ANOVA with Bonferroni post-hoc test. ns, not significant.

In order to determine the functional consequences of miR-24, miR-186, and miR-455 expression on Aβ production, we performed ELISA using HEK293-APPSwe cells. Both miR-186 and miR-455 decreased (soluble) Aβ40 and Aβ42 levels, while miR-24 had a small, nonetheless significant effect on Aβ42 (Figure 2). The observations correlate with our western blot data (Figure 1B), showing that highest NCSTN repression leads to the lowest Aβ secretion.

FIGURE 2

FIGURE 2. miR-24, miR-186, and miR-455 expression results in decreased Aβ secretion. Relative quantities of human Aβ40 and Aβ42 peptides measured by ELISA in the culture medium of HEK-APPSwe cells after transfection with 50 nM pre-miRNAs (as indicated). Signals were normalized to the basal Aβ secretion found in cells transfected with SCR miRNAs. Statistical significance was assessed by one-way ANOVA with Bonferroni post-hoc test. *p < 0.05, **p < 0.01, ***p < 0.001.

Next, we tested whether these miRNAs were functionally affected by 3′UTR polymorphisms. To address this, we generated human NCSTN 3′UTR luciferase constructs containing the C460T (rs1043329), T623G (rs113810300), or delCA515–516 (rs141849450) variants. As before, these mutant vectors were co-transfected with candidate miRNAs and luciferase activity was measured. These screens indicated that the polymorphism C460T, which is located in the 5′ compensatory region of the miRNA binding site, did not significantly affect miR-24 function (Figure 3A). On the other hand, both seed region SNPs T623G and delCA515–516 reduced miR-186 and miR-455-mediated repression, respectively (Figures 3B–D). It is interesting to note that both miR-186 and miR-455 have previously been shown to have an altered expression level in AD cerebrospinal fluid (CSF) samples, which further underscores the possible association between these miRNAs and AD development. Taken together, we identified miR-186 and miR-455, which are functionally affected by polymorphisms.

FIGURE 3

FIGURE 3. SNP delCA 515–516 and SNP T623G reduce miR-455- and miR-186-mediated NCSTN repression. (A–C) HeLa cells were transfected with 25 nM pre-miRs (as indicated) as well as a reporter construct containing the WT, delCA 515–516, T623G and C460T mutated 3′UTR of hNCSTN. The cells were lyzed 24 h post-transfection and luciferase signal was measured. After normalization for transfection efficiency, the signals were compared to the SCR control. Representative results (n= 3, performed in triplicate) are shown. Statistical significance was assessed by Student’s paired t-test (*p < 0.05, ***p < 0.001, ns, not significant). SCR; scrambled. (D) Schematic representation of base pair matching between miRNAs and the 3′UTR of hNCSTN. The seed region of the miRNAs is indicated. The red bases represent the SNPs C460T (upper panel), T623G (middle panel), and delCA 515–516 (lower panel).

We finally asked whether identified SNPs could be associated with AD. According to the NCBI dbSNP database (February 2014), the rs113810300 SNP has no reported frequency, while the rs141849450 SNP has a 1/1000 frequency in the European reference panel (1000 Genomes Phase I V3 SHAPEIT2, December 2013). Using imputation from this latter panel, rs141849450 has an estimated frequency of 5.6/1000 in an independent French cohort of more than 2000 individuals (data not shown). To validate and extend these observations, we performed TaqMan SNP genotyping assays in a Canadian patient population of AD cases and non-demented controls. These screens identified 1 out of 511 AD individuals harboring the rs141849450 deletion mutation. This individual was aged 65 years and was a carrier of the APOE4 (e3/e4) allele. Of notice, the rs141849450 SNP was not present in 631 control. No rs113810300 polymorphism was identified in either groups, as everyone was homozygous T/T. Given the low frequencies of these SNPs in the general population, much larger cohorts are required to determine whether candidate polymorphisms are indeed associated with disease. Nevertheless, we identified at least one SNP co-segregating with an AD patient, which suggests that rare polymiRTS could contribute to risk for disease in a subset of patients. In theory, the rs141849450 SNP would lead to increased NCSTN expression and Aβ load.

Discussion

In the current study, we identified a number of brain-expressed miRNAs that could regulate endogenous NCSTN expression and Aβ peptide production. Furthermore, we identified rare 3′UTR polymorphisms that could affect miRNA repressional activity toward NCSTN, at least in our experimental settings. It has been demonstrated recently that, besides expression variation in miRNAs and genetic variants in miRNA loci, polymiRTS (including both SNPs and small insertions and deletions) are important sources of expression variation of miRNA target genes in humans (Gamazon et al., 2012). Moreover, several SNPs that are in linkage disequilibrium with those that might affect miRNA binding have been identified and are associated with the risk of developing disorders including breast cancer, asthma, and Parkinson’s disease (Richardson et al., 2011; Thomas et al., 2011). These observations provide significant evidence for the role of genetic-variation-mediated regulation by miRNAs in the etiology of a broad spectrum of human complex phenotypes, including neurodegenerative diseases, metabolic traits, and autoimmune disorders. However, whether the SNPs identified in this study could also confer risk for AD remains an interesting possibility but will require extensive follow-up studies in larger cohorts of patients. Nevertheless, these observations strengthen the potential role of polymiRTS in neurodegenerative disorders such as AD.

There are some caveats associated with our study. For instance, the luciferase reporter system, although a robust, quantitative method of analysis, has limited clinical relevance. Indeed, the miRNome of immortalized cell lines and a human brain cell are definitely different, making it difficult to account for compensatory or cumulative effects by other miRNAs on the same miRNA target site. Secondly, the basal levels of the miRNA of interest as well as the epigenetics of the region at study might be different between brain and cell lines, possibly resulting in different responses. Validation in vivo would be more accurate but would require extensive screening (gene sequencing) studies in patients with AD. Unfortunately, our Canadian patient samples were limited to DNA, so we could not perform in vivo Nicastrin expression studies on mRNA. Another alternative would be to generate multiple knock-in cell lines or animals (e.g., mice) harboring candidate polymorphisms within the endogenous NCSTN 3′UTR, which extends the current scope of this study. It should be noticed that current mouse models expressing human NCSTN do not contain the (full) 3′UTR. Moreover, it remains uncertain if the identified miRNAs, although expressed in human brain, are similarly expressed in the mouse brain. Nevertheless, we do provide functional data in human cells which is in line with our main hypothesis that specific miRNAs can regulate Aβ production.

What are the clinical implications of these findings? As mentioned above, this is a notoriously difficult question to answer given the low frequencies of identified SNPs (e.g., 5/1000 for rs141849450) in the general population. These results are somewhat expected and consistent with the notion that functional 3′UTR polymorphisms are uncommon, particularly in the seed region (Chen and Rajewsky, 2006; Saunders et al., 2007; Yu et al., 2007), likely because of evolutionary pressure. They are nevertheless important for human disease, as polymiRTS have been associated with a wide range of diseases, including neurodegenerative disorders such as Parkinson’s disease (Wang et al., 2008). Thus, larger cohorts are required to determine whether candidate polymorphisms might be associated with AD. Of notice, we also explored the frequency and possible correlation between Aβ load and the rs141849450 Nicastrin variant in a large French (European) cohort. Unfortunately, we could not conclude anything due to the lack of SNP-bearing subjects (unpublished observations). We do, however, provide experimental validation that the rs141849450 Nicastrin variant exists in at least one AD patient (Canadian cohort). This is important since not based only on prediction algorithms and imputation values. Whether this has clinical relevance is interesting, but it requires extensive follow-up studies in much larger cohorts of patients.

In conclusion, we identified NCSTN-targeting miRNAs (miR-24, miR-186, and miR-455) that could decrease Aβ secretion. MiR-186 and miR-455 were previously shown to be altered in AD CSF samples, and we demonstrated that the miRNA-mediated repression of NSCTN mRNA is altered by the presence of SNPs rs113810300 and rs141849450. Our data provide additional proof of principle that polymiRTS could contribute to neurodegenerative disorders, and that miR-186 and miR-455 could be important players in AD pathology. We also stress the importance of including non-translated regions in genome-wide association studies.

Author Contributions

Charlotte Delay, G.-Y. Hsiung, Sébastien S. Hébert: design of the study. Charlotte Delay, Véronique Dorval, Alice Fok: experiments. Benjamin Grenier-Boley, Jean-Charles Lambert, Charlotte Delay: statistical analyses. Charlotte Delay, Sébastien S. Hébert: manuscript writing. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank Claudia Goupil for technical expertise. This work was supported by the CIHR and the Alzheimer Society of Canada, including a fellowship for Charlotte Delay.

References

Ambros, V. (2004). The functions of animal microRNAs. Nature 431, 350–355. doi: 10.1038/nature02871

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Anonymous. (1994). Canadian study of health and aging: study methods and prevalence of dementia. CMAJ 150, 899–913.

Betel, D., Wilson, M., Gabow, A., Marks, D. S., and Sander, C. (2008). The microRNA.org resource: targets and expression. Nucleic Acids Res. 36, D149–D153. doi: 10.1093/nar/gkm995

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chen, K., and Rajewsky, N. (2006). Natural selection on human microRNA binding sites inferred from SNP data. Nat. Genet. 38, 1452–1456. doi: 10.1038/ng1910

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Delay, C., Calon, F., Mathews, P., and Hébert, S. S. (2011). Alzheimer-specific variants in the 3′UTR of Amyloid precursor protein affect microRNA function. Mol. Neurodegener. 6, 70. doi: 10.1186/1750-1326-6-70

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Delay, C., Mandemakers, W., and Hébert, S. S. (2012). MicroRNAs in Alzheimer’s disease. Neurobiol. Dis. 46, 285–290. doi: 10.1016/j.nbd.2012.01.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

De Strooper, B. (2003). Aph-1, Pen-2, and Nicastrin with Presenilin generate an active gamma-Secretase complex. Neuron 38, 9–12. doi: 10.1016/S0896-6273(03)00205-8

CrossRef Full Text

Edbauer, D., Winkler, E., Haass, C., and Steiner, H. (2002). Presenilin and nicastrin regulate each other and determine amyloid beta-peptide production via complex formation. Proc. Natl. Acad. Sci. U.S.A. 99, 8666–8671. doi: 10.1073/pnas.132277899

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Feldman, H., Levy, A. R., Hsiung, G.-Y., Peters, K. R., Donald, A., Black, S. E.,et al. (2003). A Canadian cohort study of cognitive impairment and related dementias (ACCORD): study methods and baseline results. Neuroepidemiology 22, 265–274. doi: 10.1159/000071189

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gamazon, E. R., Ziliak, D., Im, H. K., LaCroix, B., Park, D. S., Cox, N. J.,et al. (2012). Genetic architecture of microRNA expression: implications for the transcriptome and complex traits. Am. J. Hum. Genet. 90, 1046–1063. doi: 10.1016/j.ajhg.2012.04.023

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Garcia, D. M., Baek, D., Shin, C., Bell, G. W., Grimson, A., and Bartel, D. P. (2011). Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs. Nat. Struct. Mol. Biol. 18, 1139–1146. doi: 10.1038/nsmb.2115

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Glenner, G. G., and Wong, C. W. (1984). Alzheimer’s disease: initial report of the purification and characterization of a novel cerebrovascular amyloid protein. Biochem. Biophys. Res. Commun. 120, 885–890. doi: 10.1016/S0006-291X(84)80190-4

CrossRef Full Text

Goo, J. S., Kim, Y. B., Shim, S. B., Jee, S. W., Lee, S. H., Kim, J. E.,et al. (2013). Nicastrin overexpression in transgenic mice lnduces aberrant behavior and APP processing. Mol. Neurobiol. 48, 232–243. doi: 10.1007/s12035-013-8453-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hébert, S. S., Wang, W.-X., Zhu, Q., and Nelson, P. T. (2013). A study of small RNAs from cerebral neocortex of pathology-verified Alzheimer’s disease, dementia with lewy bodies, hippocampal sclerosis, frontotemporal lobar dementia, and non-demented human controls. J. Alzheimers Dis. 35, 335–348. doi: 10.3233/JAD-122350

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hsiung, G.-Y., Sadovnick, A. D., and Feldman, H. (2004). Apolipoprotein E epsilon4 genotype as a risk factor for cognitive decline and dementia: data from the Canadian Study of Health and Aging. CMAJ 171, 863–867. doi: 10.1503/cmaj.1031789

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lupton, M. K., Proitsi, P., Danillidou, M., Tsolaki, M., Hamilton, G., Wroe, R.,et al. (2011). Deep sequencing of the Nicastrin gene in pooled DNA, the identification of genetic variants that affect risk of Alzheimer’s disease. PLoS ONE 6:e17298. doi: 10.1371/journal.pone.0017298

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

McKhann, G., Drachman, D., Folstein, M., Katzman, R., Price, D., and Stadlan, E. M. (1984). Clinical diagnosis of Alzheimer’s disease: report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer’s Disease. Neurology 34, 939–944. doi: 10.1212/WNL.34.7.939

CrossRef Full Text

Nyabi, O., Bentahir, M., Horré, K., Herreman, A., Gottardi-Littell, N., Van Broeckhoven, C.,et al. (2003). Presenilins mutated at Asp-257 or Asp-385 restore Pen-2 expression and Nicastrin glycosylation but remain catalytically inactive in the absence of wild type Presenilin. J. Biol. Chem. 278, 43430–43436. doi: 10.1074/jbc.M306957200

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pillai, R. S. (2005). MicroRNA function: multiple mechanisms for a tiny RNA? RNA 11, 1753–1761. doi: 10.1261/rna.2248605

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Rademakers, R., Eriksen, J. L., Baker, M., Robinson, T., Ahmed, Z., Lincoln, S. J.,et al. (2008). Common variation in the miR-659 binding-site of GRN is a major risk factor for TDP43-positive frontotemporal dementia. Hum. Mol. Genet. 17, 3631–3642. doi: 10.1093/hmg/ddn257

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Richardson, K., Lai, C.-Q., Parnell, L. D., Lee, Y.-C., and Ordovas, J. M. (2011). A genome-wide survey for SNPs altering microRNA seed sites identifies functional candidates in GWAS. BMC Genomics 12:504. doi: 10.1186/1471-2164-12-504

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Saunders, M. A., Liang, H., and Li, W.-H. (2007). Human polymorphism at microRNAs and microRNA target sites. Proc. Natl. Acad. Sci. U.S.A. 104, 3300–3305. doi: 10.1073/pnas.0611347104

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Thomas, L. F., Saito, T., and Sætrom, P. (2011). Inferring causative variants in microRNA target sites. Nucleic Acids Res. 39:e109. doi: 10.1093/nar/gkr414

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wang, G.,van Der Walt, J. M., Mayhew, G., Li, Y.-J., Züchner, S., Scott, W. K.,et al. (2008). Variation in the miRNA-433 binding site of FGF20 confers risk for Parkinson disease by overexpression of alpha-synuclein. Am. J. Hum. Genet. 82, 283–289. doi: 10.1016/j.ajhg.2007.09.021

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yu, Z., Li, Z., Jolicoeur, N., Zhang, L., Fortin, Y., Wang, E.,et al. (2007). Aberrant allele frequencies of the SNPs located in microRNA target sites are potentially associated with human cancers. Nucleic Acids Res. 35, 4535–4541. doi: 10.1093/nar/gkm480

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text