| Literature DB >> 19776031 |
Glen M Borchert1, Brian L Gilmore, Ryan M Spengler, Yi Xing, William Lanier, Debashish Bhattacharya, Beverly L Davidson.
Abstract
Animals regulate gene expression at multiple levels, contributing to the complexity of the proteome. Among these regulatory events are post-transcriptional gene silencing, mediated by small non-coding RNAs (e.g. microRNAs), and adenosine-to-inosine (A-to-I) editing, generated by adenosine deaminases that act on double-stranded RNA (ADAR). Recent data suggest that these regulatory processes are connected at a fundamental level. A-to-I editing can affect Drosha processing or directly alter the microRNA (miRNA) sequences responsible for mRNA targeting. Here, we analyzed the previously reported adenosine deaminations occurring in human cDNAs, and asked if there was a relationship between A-to-I editing events in the mRNA 3' untranslated regions (UTRs) and mRNA:miRNA binding. We find significant correlations between A-to-I editing and changes in miRNA complementarities. In all, over 3000 of the 12 723 distinct adenosine deaminations assessed were found to form 7-mer complementarities (known as seed matches) to a subset of human miRNAs. In 200 of the ESTs, we also noted editing within a specific 13 nucleotide motif. Strikingly, deamination of this motif simultaneously creates seed matches to three (otherwise unrelated) miRNAs. Our results suggest the creation of miRNA regulatory sites as a novel function for ADAR activity. Consequently, many miRNA target sites may only be identifiable through examining expressed sequences.Entities:
Mesh:
Substances:
Year: 2009 PMID: 19776031 PMCID: PMC2778373 DOI: 10.1093/hmg/ddp443
Source DB: PubMed Journal: Hum Mol Genet ISSN: 0964-6906 Impact factor: 6.150
Figure 1.ADARs deaminate adenosine to inosine, potentially altering miRNA complementarities. (A) A cartoon depicting adenosine, deaminated adenosine (inosine) and guanine. (B) A representative deamination site occurring in the 3′-UTR of DNA fragmentation factor α (DFFA) is shown in both the unedited (left) and edited (right) state. For this transcript, the 7 nt sequence scan (red) of the 200 nts flanking the deamination site (blue) identified a miR-513 seed sequence (yellow). This was repeated for all miRNA seeds (miRBase 13.0) for all known edited transcripts.
miRNA seed matches markedly enriched by adenosine deamination
| miRNA | Seed | Totala/createdb/expectedc | |
|---|---|---|---|
| mir-513a-5pd | TCACAGG | 377/258/13.5 | 2.0E−226 |
| mir-769-3p/-450b-3pd | TGGGATC | 400/252/14.4 | 1.2E−212 |
| mir-140-3p | ACCACAG | 348/139/12.5 | 1.2E−92 |
| mir-129-5p | TTTTTGC | 322/105/11.6 | 4.0E−62 |
| mir-936 | CAGTAGA | 269/94/9.7 | 2.5E−58 |
| mir-412 | CTTCACC | 199/83/7.1 | 1.6E−57 |
| mir-326/-330 | CTCTGGG | 196/82/7.0 | 6.1E−57 |
| mir-646 | AGCAGCT | 311/91/11.2 | 2.7E−50 |
| mir-371-3p | AGTGCCG | 97/57/3.5 | 6.2E−48 |
| mir-763/-1207-3p | CAGCTGG | 413/99/14.8 | 3.9E−47 |
| mir-1229 | TCTCACC | 256/74/9.2 | 6.8E−41 |
| mir-630 | GTATTCT | 146/56/5.2 | 1.6E−37 |
| mir-325-5p | CTAGTAG | 161/57/5.8 | 2.3E−36 |
| mir-99ab/-100 | ACCCGTA | 46/27/1.7 | 1.4E−23 |
| mir-518a-dfe | AAAGCGC | 57/27/2.1 | 3.2E−21 |
| mir-643 | CTTGTAT | 66/28/2.4 | 1.0E−20 |
| mir-518ee | AAGCGCT | 61/27/2.2 | 1.7E−20 |
| mir-548c-3p | AAAAATC | 419/59/15.0 | 8.1E−18 |
| mir-578 | TTCTTGT | 140/32/5.0 | 7.8E−16 |
| mir-28-5p/-708 | AGGAGCT | 176/35/6.3 | 2.2E−15 |
| mir-149 | CTGGCTC | 160/30/5.7 | 8.8E−13 |
| mir-1306 | AGGAGCT | 176/28/6.3 | 2.0E−10 |
| mir-1289 | GGAGTCC | 92/20/3.3 | 4.3E−10 |
| mir-344-5p/-484 | CAGGCTC | 311/36/11.2 | 2.9E−9 |
| mir-324-3p | CTGCCCC | 134/22/4.8 | 9.3E−9 |
| mir-615-3p | CCGAGCC | 29/11/1.0 | 1.5E−8 |
| mir-1281 | CGCCTCC | 754/57/27.1 | 3.6E−7 |
| mir-1298 | TCATTCG | 5/5/0.2 | 1.3E−6 |
| mir-197 | TCACCAC | 440/38/15.8 | 1.5E−6 |
| mir-1247 | CCCGTCC | 15/7/0.5 | 1.6E−6 |
aTotal seed sites were tallied by counting the number of seed matches for each miRNA (miRBase 13.0) in every 7-mer window, in all 12 723 edited (central ‘G’) sequences. Results for remaining miRNAs and for unedited counts are provided in Supplementary Material, Table S1.
bBecause only 7 out of 195 7-mer windows (windows 95–101) are altered by editing, ‘created’ sites were determined by tallying the number of seed matches in these seven windows in the edited sequences.
cExpected values were calculated as seven times the number of miRNA seed matches found in any 7-mer window if the total seed sites were equally distributed [total sites × (7/195)].
dmiR-513 and miR-769-3p/-450b-3p are often (∼80%) complementary to the same deaminated sequences.
eThe two miR-518 family seeds ‘AAAGCGC’ and AAGCGCT’ are complementary to the same 27 deaminated sequences.
*P-values are calculated based on Poisson distribution and the cutoff (1.85E−5) determined by Bonferroni correction on a familywise error rate of 0.01.
Figure 2.A-to-I edits create miR-513 and miR-769-3p/-450b-3p complementarities. (A) 12 723 unique EST sequences (www.cgen.com), each consisting of a central A-to-I deamination and 100 nt flanks [i.e. n100 (A or I) n100], were screened for complementarity to human miRNAs. All human miRNA seed matches were identified within the individual 201 nt sequences originally identified as an A-to-I transition by Compugen (statistical significance is addressed in Table 1). The top two panels represent all miR-769-3p (and miR-450b-3p) seed matches occurring at each position in both the unedited (left) and edited (right) states. The lower panels represent all miR-513 seed matches occurring in unedited (left) and edited (right) states. (B) A cartoon of miR-513 and miR-769-3p complementarities to a miRNA Associating If Deaminated (MAID) site in both unedited (left) and edited (right) state is shown. The seed of miR-450b-3p (not shown) is identical to miR-769-3p. Perfect seed matches to miR-769-3p/-450b-3p (blue) and miR-513 (yellow) are significantly enriched in sequences containing characterized deaminations (red). Vertical lines indicate complementary base pairing. (C) Venn diagram depicting the overlap between miR-513 and miR-769-3p/-450b-3p target sites matching the full 13-bp MAID motif created by adjacent miR-513 and miR-769/-450-3p seeds and a G:U wobble. Nearly 100 additional sequences are identified by allowing for the G:U wobble immediately 3′ to the deamination.
Figure 3.MiR-513 and miR-769-3p target MAID sites but not the corresponding unedited sequence. (A) A diagram shows hairpin expression vectors and MAID reporter constructs. pAL-513 and pAL-769-3p reporter vectors have miR-513 and miR-769-3p hairpins downstream of the miR-517 Pol-III promoter (28). TAAT, TGAT and TGGT reporters contain three tandem copies of the 13 bp MAID sequence in the 3′-UTR of Renilla luciferase for testing activity in the unedited (TAAT) or edited (TGAT and TGGT) states. Guanines mimicking A-to-I edits are bolded and underscored. (B) Renilla luciferase activity (normalized to firefly luciferase and presented as percent mock transfected control) following co-transfection of miR-513, miR-769-3p, pooled miR-513 and miR-769-3p inhibitors and/or control miRNA inhibitor with the indicated reporters into HEK 293 cells (n = 3) is illustrated. *P < 0.005.
Figure 4.Endogenous MAIDs are targets for miR-513 and miR-769-3p repression. (A) A cartoon depicts the DFFA 3′-UTR and the localization of nine distinct MAIDs (lines above the 3′-UTR). (B) Alignment of the nine DFFA 3′-UTR MAIDs commonly deaminated in ESTs is represented. MAID sequences are shaded. Four MAIDs contain G:U wobbles from consensus (bold). (C) Alignment of DFFA_1 sequences from three independent DFFA clones isolated from HEK 293 cells and NB7 cells is shown. DFFA_1 was deaminated in NB7s (bold) but not in HEK 293 cells. RT reactions were performed using a thermostable reverse transcriptase. (D) A diagram of DFFA 3′-UTR reporter constructs is shown. In DFFA-Edited (-E) and DFFA-Unedited (-U), the Renilla luciferase 3′-UTRs are the cloned DFFA 3′-UTRs from NB7 and HEK 293 cells, respectively. DFFA-E nucleotides differing from DFFA-U are bolded and underscored (compare NB7_2 and 293_1 detailed in panel (C). (E) Luciferase assays performed identically to those in Fig. 3B except for the reporter constructs illustrated (n = 3). *P < 0.005.
Figure 5.MiR-769-3p selectively represses DFFA protein. (A) Relative mature miR-769-3p levels in HEK 293, A549, HT1080 and NB7 cell lines are shown as determined by quantitative PCR. MiR-513 was not detected in these cell lines. (B) Western blot analysis of endogenous DFFA in HEK 293 and NB7 cell lysates following transfection of miR-769-3p as indicated. Representative blots for DFFA and β-catenin (loading control) are shown. (C) Relative DFFA levels calculated as band intensity ratios of DFFA to β-catenin and normalized to mock (left most bar in each graph). 400, 200 and 100 refer to ng of miR-769-3p expression vector.