Genomics Proteomics Bioinformatics 12 (2014) 92–104
Genomics Proteomics Bioinformatics
www.elsevier.com/locate/gpb www.sciencedirect.com
ORIGINAL RESEARCH Mining the 30UTR of Autism-implicated Genes for SNPs Perturbing MicroRNA Regulation
Varadharajan Vaishnavi #, Mayakannan Manikandan #, Arasambattu Kannan Munirajan *
Department of Genetics, Dr. ALM PG Institute of Basic Medical Sciences, University of Madras, Taramani Campus, Chennai 600113, India
Received 13 September 2013; revised 9 December 2013; accepted 11 January 2014 Available online 18 April 2014
Handled by Zhongming Zhao
KEYWORDS Abstract Autism spectrum disorder (ASD) refers to a group of childhood neurodevelopmental dis- Autism spectrum disorder; orders with polygenic etiology. The expression of many genes implicated in ASD is tightly regulated by MicroRNA; various factors including microRNAs (miRNAs), a class of noncoding RNAs 22 nucleotides in Single nucleotide polymor- length that function to suppress translation by pairing with ‘miRNA recognition elements’ (MREs) phism; present in the 30untranslated region (30UTR) of target mRNAs. This emphasizes the role played by 30UTR-SNP; miRNAs in regulating neurogenesis, brain development and differentiation and hence any perturba- Gene regulation tions in this regulatory mechanism might affect these processes as well. Recently, single nucleotide polymorphisms (SNPs) present within 30UTRs of mRNAs have been shown to modulate existing MREs or even create new MREs. Therefore, we hypothesized that SNPs perturbing miRNA-medi- ated gene regulation might lead to aberrant expression of autism-implicated genes, thus resulting in disease predisposition or pathogenesis in at least a subpopulation of ASD individuals. We developed a systematic computational pipeline that integrates data from well-established databases. By following a stringent selection criterion, we identified 9 MRE-modulating SNPs and another 12 MRE-creating SNPs in the 30UTR of autism-implicated genes. These high-confidence candidate SNPs may play roles in ASD and hence would be valuable for further functional validation.
Introduction * Corresponding author. E-mail: [email protected], [email protected] (Munirajan AK). Autism spectrum disorder (ASD) represents a series of neuro- # Equal contribution. developmental disorders characterized by altered social inter- Peer review under responsibility of Beijing Institute of Genomics, ests, communicative deficits, and restricted and repetitive Chinese Academy of Sciences and Genetics Society of China. behaviors, usually with an onset before the age of three [1]. Although ASD is defined by a triad of symptoms, the levels of severity and presentation may vary among individuals, Production and hosting by Elsevier demonstrating the tremendous heterogeneity of the condition.
1672-0229/$ - see front matter ª 2014 Beijing Institute of Genomics, Chinese Academy of Sciences and Genetics Society of China. Production and hosting by Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.gpb.2014.01.003 Vaishnavi V et al / MicroRNA Binding Site Polymorphisms in Autism 93
The burden of ASD is increasing worldwide with recent prev- At this point, computational approaches are still the driving alence estimates of 62 per 10,000 children [2]. Research in the force in identifying truly functional genetic variants within last several years has substantiated the strong genetic basis and MREs, and numerous databases have been developed to assist has further established the polygenic etiology of ASD [3,4]. in this process. Aided by a rigorous bioinformatics approach, For instance, a comprehensive transcriptomic analysis was we identified a panel of high-confidence SNPs that either mod- performed to measure mRNA levels in post-mortem brains ulate existing MREs or create new MREs in the 30UTR of aut- from autistic individuals and controls using Illumina micro- ism-implicated genes. We propose that alterations in the arrays. As a result, 444 genes were found to be differentially miRNA:MRE interaction caused by these MRE-SNPs could expressed [5]. This disease-specific, aberrant expression of be an alternate mechanism underlying aberrant gene expres- genes may be caused by various factors, such as the presence sion, which may contribute to predisposition, pathogenesis, of submicroscopic structural chromosomal alterations in the inter-individual variation in gene expression, and heterogene- genome denoted as copy number variants (CNVs), or by per- ity of ASD (Figure 1). Future experimental validations may turbations in many regulatory mechanisms that govern gene further clarify the effect and consequences of these candidate expression [6–8]. Very recently, methylomic analysis of SNPs. monozygotic twins discordant for ASD identified epigenetic variations that mediate disease susceptibility via altered gene Results dosage [9]. Besides gene dysfunction, studies reveal that the known etiological factors of ASD eventually converge toward the unifying theme of aberrant gene expression [5–9]. Study design and preliminary analyses MicroRNAs (miRNAs) have emerged as key regulators for nearly 50% of protein-coding genes in the human genome, We followed a stepwise, integrative and stringent computa- participating in all vital cellular processes [10–12]. miRNAs tional pipeline to identify and characterize MRE-SNPs in the represent a class of short ( 22 nt) endogenous noncoding 30UTR of autism-implicated genes (Figure 2). A total of 484 RNAs that pairs with complementary sequences called miR- genes were retrieved from the ‘‘Human Gene Module’’ of NA recognition elements (MREs), which are located mostly ‘SFARI Gene’, a comprehensive resource for obtaining reli- in the 30untranslated region (30UTR) of target mRNAs. This able data on the genetics of ASD [30]. SNPs that modulate miRNA:MRE pairing leads to translational inhibition or or create MREs across the 30UTR of these genes were identi- mRNA destabilization by recruiting the RNA-induced silenc- fied using MirSNP database [31]. The resultant list contained ing complex (RISC). In mammals, complete complementarity 13,945 SNPs in 468 out of the 484 autism-implicated genes. between miRNA and MRE is rare, since a minimal 6-bp match We chose MirSNP due to its high sensitivity in covering a of nucleotides 2–8 at the 50 end of miRNA (known as the seed majority of the validated miRNA-related SNPs when com- region) is sufficient for functionality [13]. Thus a single miR- pared to other databases [31]. MirSNP also provides informa- NA can target approximately 200 transcripts, and more than tion on the plausible effect of SNPs on miRNA binding and one miRNA can act upon a single mRNA target [14]. Impor- categorizes them into one of the four functional classes. These tantly, almost 70% of experimentally-detectable miRNAs are include (i) break –– SNP completely disrupts the binding of expressed in the human nervous system [15]. Several miRNAs miRNA to the MRE, (ii) decrease –– SNP reduces the binding orchestrate a myriad of diverse neurodevelopmental processes efficacy of the miRNA to the MRE, (iii) enhance –– SNP including neuronal cell specification and differentiation, synap- increases the binding affinity of the miRNA to the MRE and tic plasticity and memory formation [16]. Accumulating evi- (iv) create –– SNP creates new MREs for miRNAs. We dence suggests that alterations in miRNA expression or grouped the first three functional classes as MRE modulating function are associated with the cognitive deficits and neurode- SNPs (MREm-SNPs) since the existing miRNA recognition velopmental abnormalities observed in autism, schizophrenia sites were modified. The fourth category that created new and other forms of intellectual dysfunction [17–20]. According MREs (MREc-SNPs) was considered as a separate group. to a recent bioinformatics study, co-expression of neural miR- Next, we checked the minor allele frequencies (MAFs) of these NAs with their target mRNAs is a common feature [21], sug- SNPs in Exome Variant Server (EVS) maintained by the gesting that any interference in miRNA:MRE interactions National Heart, Lung and Blood Institute (NHLBI) Exome might affect the physiological control of gene expression. Sequencing Project. The EVS was chosen for several reasons. Additionally, single nucleotide polymorphisms (SNPs) in the First, the database serves as a large repository for more than 30UTR of genes are recognized for their ability to perturb miR- 10 million human SNPs identified by sequencing 15,336 genes NA binding by modulating existing MREs or by creating new in 6515 individuals of European American (EA) and African MREs [22–25]. The pathological significance of this unique American (AA) ancestry [32]. Second, even SNPs with MAF class of functional polymorphisms is gaining recognition and as low as 0.1% can be identified owing to the larger sample their biological relevance has been explored in various diseases size. Third, EVS also indicates the level of conservation and [26,27], including neurodegenerative diseases [28,29]. To our evolutionary constraints of the SNP loci by providing phyloge- knowledge, MRE-SNPs have not formed the basis of any pre- netic analysis with space/time model for conservation (phast- vious case-control studies in autism, and their systematic Cons) and genomic evolutionary rate profiling (GERP) assessment might extend our understanding of the complex scores. Finally, for many SNPs, the MAF is unavailable in genetic architecture of this highly heterogeneous disorder. NCBI’s dbSNP. Of the 13,945 SNPs we identified, it was inter- However, the identification of ASD-associated MRE-SNPs esting to note that only 387 spanning the 30UTR of 181 autism- poses a considerable challenge due to the lack of an established implicated genes were documented with a MAF > 0.1% in at model that accounts for the heterogeneity of ASD and the least one population (EA or AA). This suggests that the vast ever-expanding modes of miRNA-mediated gene regulation. majority of MRE-SNPs are either very rare or not detected 94 Genomics Proteomics Bioinformatics 12 (2014) 92–104
3′
Figure 1 SNPs in MREs of autism-implicated genes may underlie ASD etiology SNPs in the 30UTR of autism-implicated genes may either modulate existing MREs or create new MREs thereby leading to aberrant expression of these genes, which may underlie predisposition, pathogenesis or genetic heterogeneity of autism. The ‘‘DNA–brain’’ image is adopted from the public domain http://neurosciencenews.com/files/2013/06/dna-brain-familial-alziheimers-public.jpg, which is credited to the National Institute of Neurological Disorders and Stroke (NINDS)/National Institutes of Health (NIH) and is appropriately modified. MRE, miRNA recognition element. in samples of the specified ancestry. Alternatively, it may indi- a functional duplex. We computed the free energy of binding cate the action of purifying selection on MRE-SNPs [33]. of the miRNA:MRE duplex by analyzing their sequences in RNAhybrid. A maximum hybridization energy of 15 kcal/ Minimizing false positive MREs mol was used as the scoring criteria to select stable duplexes. All of these sequential filtering steps aided the identification The first step in the identification of MREm-SNPs relies on the of high-confidence functional MREs and thus, MREm-SNPs. correct prediction of functional MREs in the 30UTR of aut- However, the biological relevance of these SNPs depends on ism-implicated genes. MirSNP database employs the miRanda whether the miRNA and its target mRNAs are co-expressed miRNA target prediction algorithm [34]. Despite the high sen- in a particular tissue or cell type. Considering the neurobio- sitivity of miRanda [35], the reliability of a single prediction logical basis of ASD, we shortlisted only those MREs that method is modest. To account for this, the MREs were sub- served as binding sites for the brain expressed human miR- jected to consistent cross-prediction by TargetScanHuman NAs. This brain-specific miRNA dataset consisted of overlap- 6.2 [36,37] and RNAhybrid 2.1 [38,39]. Furthermore, we ping miRNAs from two independent studies that profiled the improved the stringency of the prediction by considering cer- expression of miRNAs in different regions of post-mortem tain quantitative scores. As an added feature, the MirSNP human brain tissue samples [41,42]. In total, 9 MREm-SNPs database provides an miRSVR score, which was developed were identified across the 30UTR of 17 transcripts (Table 1). by supervised training on mRNA expression data from a panel The effect on miRNA binding was verified by computing the of miRNA transfection experiments [40]. We considered only difference in hybridization energy (DDG) caused by that par- the top miRSVR predictions with a score of <