<title>Introduction</title> <p>Alternative splicing enhances transcriptomic diversity and presumably leads to speciation and higher organism complexity, especially in mammals ##REF##11173120##[1]##–##REF##15510168##[3]##. There are four major types of alternative splicing: exon skipping, which is the most prevalent form in higher vertebrates; alternative 5′ and 3′ splice site (5′ss and 3′ss) selection; and intron retention, which is the rarest form in both vertebrates and invertebrates ##REF##17158149##[4]##,##REF##14992493##[5]##. At least 74% (and probably much more) of human genes that contain introns produce more than one type of mRNA transcript through alternative splicing; however, it is unclear which of these products are biologically functional and which are non-functional products of inaccurate splicing ##REF##15510168##[3]##, ##REF##14684825##[6]##–##REF##15851065##[8]##. Thus, understanding the changes in the genome that dictate fixation of beneficial alternative splicing events or deleterious events (e.g., mutations leading to genetic disorders or cancer), or aberrant splicing events (noise in the system) is of great interest.</p> <p>There are three known origins of alternatively spliced exons: 1) exon shuffling, which is a form of gene duplication ##REF##622185##[9]##–##REF##12045209##[11]##; 2) exonization of intronic sequences ##REF##12097342##[12]##–##REF##17594509##[16]##; and 3) change in the mode of splicing from constitutive to alternative splicing during evolution ##REF##18020709##[17]##,##REF##17530917##[18]##. One mechanism responsible for the shift from constitutive to alternative splicing is accumulation of mutations in the 5′ splice site region. Here we set out to examine additional mechanisms involved in the transition from constitutive to alternative splicing.</p> <p>The primate-specific retrotransposons called <italic>Alu</italic> are ∼280 nucleotides long. These are the most abundant retrotransposed elements in the human genome with about 1.1 million copies ##REF##17594509##[16]##,##REF##11988762##[19]##,##REF##18046749##[20]##. A large fraction of these <italic>Alu</italic> elements are located within intronic sequences, in both the sense and the antisense orientation relative to the mRNA, and can potentially form long regions of double-stranded RNA (dsRNA) ##REF##17594509##[16]##, ##REF##15534692##[21]##–##REF##12048240##[25]##. There are indications that extensive secondary structure occurs between <italic>Alu</italic> elements. The evidence is embedded in analyses of the RNA editing mechanism: The human transcriptome undergoes extensive adenosine to inosine RNA editing ##REF##15342557##[23]##,##REF##16940548##[26]##,##REF##15661352##[27]##. RNA editing is directed by an adenosine deamination mechanism catalyzed by specific adenosine deaminases, termed dsRADs (double-stranded RNA adenosine deaminases) or ADARs ##REF##16940548##[26]##,##REF##9292492##[28]##,##UREF##0##[29]##. ADARs are required for the formation of the dsRNA molecules that serve as substrates for the deamination process ##REF##8269514##[30]##,##REF##7937939##[31]##. Hence, dsRNA regions formed between two <italic>Alu</italic>s in opposite orientation within 2000 nucleotides of each other may serve as substrates for ADAR ##REF##15534692##[21]##–##REF##12048240##[25]##,##REF##15661352##[27]##,##REF##9292492##[28]##,##REF##15833131##[32]##.</p> <p>More than 90% of known editing sites are found in <italic>Alu</italic> elements and editing occurs in sense and antisense pairs of <italic>Alu</italic>s but not in flanking non-<italic>Alu</italic> sequences ##REF##15342557##[23]##,##REF##15258596##[24]##,##REF##16574103##[33]##. Recently, it was shown that a pair of inverted <italic>Alu</italic>s located within the 3′UTR of <italic>EGFP</italic> mRNA serves as a substrate for A-to-I RNA editing that stabilizes the binding of the p54 protein to the mRNA. This causes nuclear retention of the mRNA and the silencing of <italic>EGFP</italic> expression ##REF##18497743##[34]##. Another example comes from the <italic>NARF</italic> gene, where formation of <italic>Alu</italic>-<italic>Alu</italic> dsRNA and its subsequent editing generates a functional 3′ splice site that is essential for exonization of that intronic <italic>Alu</italic>; moreover, editing in that <italic>Alu</italic> eliminates a stop codon and modulates the strength of exonic splicing regulatory sequences (ESRs). Interestingly, the nucleotides surrounding the editing site are important not only for editing of that particular site but also for editing at other sites located downstream in the same exon. It was also shown that the C nucleotide thought to pair with the edited site on the dsRNA is important for editing ##REF##17326827##[35]##. There is emerging evidence that the secondary structure of precursor mRNA plays a role in regulation of alternative splicing. However, in most studies the double-stranded structure was made up of only 10–40 base pairs and sequestered exonic or splice site sequences ##REF##17580311##[36]##–##REF##18020710##[47]##.</p> <p>In this study, we have bioinformatically and experimentally evaluated the effects of intronic <italic>Alu</italic> elements on splicing. We found that different regulatory constraints act on <italic>Alu</italic> insertions into introns that flank constitutively or alternatively spliced exons. We further demonstrated that two <italic>Alu</italic> elements which were inserted into introns in opposite orientation have the potential to undergo base-pairing, as evident by RNA editing, and affect the splicing patterns of a downstream exon, by shifting it from constitutive to alternative. Finally, as <italic>Alu</italic> elements are abundant in introns, the findings we present suggest that the effect of intronic <italic>Alu</italic> elements on the transcriptome could be substantial, and could result in transcriptomic novelties. The new isoforms could then be subjected to purifying selection which will determine their fixation.</p>
<title>Materials and Methods</title> <title>Dataset Compilation</title> <p>A dataset of 596 alternatively spliced exons, conserved between human and mouse, was derived from a previously compiled dataset ##REF##14992493##[5]##. In addition, 45,553 human-mouse conserved constitutively spliced exons were obtained from Carmel et al. ##REF##15100438##[74]##. Species-specific exons (354) were extracted from a dataset of 4,262 human-mouse orthologous exons that are suspected to splice differently in human and mouse based on initial EST analysis ##REF##15100438##[74]##. For details of how the datasets were built see ##REF##18020709##[17]##.</p> <title>Human and Mouse Intronic Dataset</title> <p>Introns and exons for human (<italic>Homo sapiens</italic>, Build 35.4) and mouse (<italic>Mus musculus</italic>, Build 34.1) were extracted from the Exon-Intron Database (<ext-link ext-link-type="uri" xlink:href="http://hsc.utoledo.edu/depts/bioinfo/database.html">http://hsc.utoledo.edu/depts/bioinfo/database.html</ext-link>) ##REF##10592221##[75]##. These intron sequences were analyzed with RepeatMasker software version 3.1.0 ##UREF##4##[76]## (<ext-link ext-link-type="uri" xlink:href="http://www.repatmasker.org">www.repatmasker.org</ext-link>) using Repbase update files ##REF##16093699##[77]##.</p> <title>Retrotransposed Elements Analysis</title> <p>Since <italic>Alu</italic>s are primate specific, the distribution was computed only from human flanking introns. The density of <italic>Alu</italic> elements was calculated per 1000-bp intron length according to the following equation:N = number of <italic>Alu</italic>s within the intron; L = the length of the intron; <italic>Alu</italic>
<sub>density</sub> = <italic>Alu</italic> density. For the detection of retrotransposed elements, we used the RepeatMasker (<ext-link ext-link-type="uri" xlink:href="http://www.repeatmasker.org">http://www.repeatmasker.org</ext-link>) software version 3.1.0 ##UREF##4##[76]## and Repbase update ##REF##16093699##[77]##.</p> <title>Statistical Analysis</title> <p>A T-test was used to calculate statistical differences between two populations; for χ<sup>2</sup> test with 2×2 contingency table, Fisher's exact test was used.</p> <title>The Potential Alternatively Skipped Exons under Intronic <italic>Alu</italic> Regulation</title> <p>The alternatively skipped exons in the human genome (build hg18) were extracted by downloading the knownAlt table from UCSC genome browser ##UREF##1##[48]##. The presence or absence of <italic>Alu</italic> within the upstream intron was determined using RepeatMasker tables downloaded from UCSC. The conservation of these introns was analyzed using MAF pairwise alignments between the human genome (build hg18) and the mouse genome (build mm9) downloaded from UCSC genome browser. The intersections between these tables were done using the Galaxy sever ##REF##16169926##[49]##.</p> <title>Bioinformatic Analysis of <italic>Alu</italic> Editing within Intronic Sequences</title> <p>We extracted data of on 10,113 nucleotides that undergo mRNA editing in the human genome from Levanon et al. ##REF##15258596##[24]##. The UCSC genome browser ##REF##18276646##[58]## was then used to extract data of all <italic>Alu</italic> elements in the human genome (build hg18) using the RepeatMasker ##UREF##4##[76]## annotations, and to extract human RefSeq intronic sequences. Intersection of these three datasets yielded 953 <italic>Alu</italic> elements that are embedded in intronic sequences and undergo mRNA editing.</p> <title>Plasmid Constructs</title> <p>The <italic>RABL5</italic> (RAB member RAS oncogene family-like 5) minigene was generated by amplifying a human genomic fragment using PCR reaction. Each primer contained an additional sequence encoding a restriction enzyme. The PCR product was restriction digested and inserted into the pEGFP-C3 plasmid (Clontech) and sequenced to confirm that the desired construct was obtained. The <italic>RABL5</italic> minigene, contains exons 2 through 4 (2.7 kb). The intron replacements with the <italic>RABL5 Alu</italic>1, <italic>Alu</italic>2, and <italic>Alu</italic>3 were done by PCR opening of the plasmid lacking the specific <italic>Alu</italic> (#1, #2 or #3) and ligation with a fragment of 270 intronic-nucleotides taken from a PCR amplification directed to the <italic>IKBKAP</italic> gene intron number 20 (primer forward, <named-content content-type="gene">
<sup>5′</sup>AGAATCGTGACACTCATCATATAAAGGAGG<sup>3′</sup>
</named-content>; and primer reverse, <named-content content-type="gene">
<sup>5′</sup>CAAAACATTAGTATAGATCTTTCCAATACA<sup>3′</sup>
</named-content>). The 800-nucleotide insertion1 was taken from PCR amplification directed to the <italic>IMP</italic> gene intron number 11 (primer forward, <named-content content-type="gene">
<sup>5′</sup>ATCACTCTGCACTTTCTCCCAT<sup>3′</sup>
</named-content>; primer reverse <named-content content-type="gene">
<sup>5′</sup>ACCATGTCCACTTCATCCAGTTC<sup>3′</sup>
</named-content>). Insertion2 is a 25-bp sequence, free of any known splicing regulatory elements, that was doubled or tripled into 50-bp and 75-bp sequences, (<named-content content-type="gene">
<sup> 5′</sup>CTATCTGATAAGCTGCGAGCAATT<sup>3′</sup>
</named-content>).</p> <title>cDNA Amplification</title> <p>Endogenous PCR amplification was done on a cDNA template originating from a neuroblastoma cell-line (SH-SY5Y). Amplification was performed for 30 cycles, consisting of denaturation for 30 seconds at 94°C, annealing for 45 seconds at 52°C or 56°C, and elongation for 1 minute at 72°C. The products were separated in a 1.5% agarose gel<bold>.</bold> The upper PCR product was Topo-ligated (Invitrogen) and sequenced. The primers used were: forward (exon 2), <named-content content-type="gene">
<sup>5′</sup>CAGAATCTTCTGACATCACTG<sup>3′</sup>
</named-content>; or forward (intron 2), <named-content content-type="gene">
<sup>5′</sup>GTGAGCCCTGACAAATCTGTGT<sup>3′</sup>
</named-content>; and reverse (exon 3) <named-content content-type="gene">
<sup>5′</sup>GTTGCTGGTAACATGCGGGTTC<sup>3′</sup>
</named-content>.</p> <title>Site-Directed Mutagenesis, Transfection, RNA Isolation, and RT-PCR Amplification</title> <p>For details see ##REF##18332115##[78]##.</p>
<title>Results</title> <title>Genome-Wide Analysis of <italic>Alu</italic>s within Intronic Sequences</title> <p>To examine potential effects of insertion of <italic>Alu</italic> elements into introns on splicing of the flanking exons, we downloaded data of human introns (hg18) and <italic>Alu</italic> elements and determined the intersected set using the UCSC genome browser and GALAXY ##UREF##1##[48]##,##REF##16169926##[49]##. Overall, 730,622 <italic>Alu</italic> elements that reside within introns and 185,534 introns were extracted. This analysis showed that there are 85,126 introns that contain at least one <italic>Alu</italic> element; of these, 5009 introns contained at least two <italic>Alu</italic> elements in opposite orientation. The median length of introns containing at least one <italic>Alu</italic> element is 3829 base pairs (bp), whereas the median length of introns that do not contain an <italic>Alu</italic> is 521 bp (for intron length distribution see ##SUPPL##0##Figure S1##). This suggests that double strand formation as a result of base pairing between two nearby <italic>Alu</italic> elements in opposite orientation might be common.</p> <title>Bioinformatic Analysis of <italic>Alu</italic>-Editing within Intronic Sequences</title> <p>An antisense <italic>Alu</italic> and a sense <italic>Alu</italic> that are within 2000 nucleotides of each other can form dsRNA and be subjected to mRNA editing ##REF##15534692##[21]##–##REF##15258596##[24]##,##REF##12045112##[50]##. To examine the configurations of possible <italic>Alu-Alu</italic> dsRNA we extracted data on 10,113 nucleotides that undergo mRNA editing in the human genome from Levanon et al. ##REF##15258596##[24]##. Intersection with data on all <italic>Alu</italic> elements in the human genome and human RefSeq intronic sequences yielded 953 <italic>Alu</italic> elements that are embedded in intronic sequences that undergo mRNA editing (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). For each of the 953 edited <italic>Alu</italic> elements, the nearest <italic>Alu</italic> element in the opposite orientation was identified. For the vast majority of edited <italic>Alu</italic> elements (880 out of the 953; 92%), the closest <italic>Alu</italic> element in the opposite orientation resided within the same intron, with an average distance of 682 bp from the edited <italic>Alu</italic>. However, we found 73 cases (8%) where the nearest <italic>Alu</italic> element in the opposite orientation was in a different intron; in 61 of these cases it is at least 500 bp closer to the edited <italic>Alu</italic> element than the nearest <italic>Alu</italic> element in the opposite orientation in the same intron (##SUPPL##1##Figure S2##). In these cases, the average distance between the edited <italic>Alu</italic> and the nearest <italic>Alu</italic> in the opposite orientation is 1993 bp. In fact, in 43 out of the 61 cases the distance was less than 2000 bp (averaging 1122 bp). This close proximity between the two <italic>Alu</italic> elements, along with the evidence that at least one of them undergoes editing, suggests that these regions may base pair.</p> <p>Since in 92% of edited <italic>Alu</italic> elements, the closest <italic>Alu</italic> element in the opposite orientation resided within the same intron, we decided to examine the splicing process in these cases. But first, we set to examine the distribution of <italic>Alu</italic> elements within datasets of exons conserved within human and mouse having different splicing patterns.</p> <title>
<italic>Alu</italic>s More Common in the Flanking Introns of Alternatively Spliced Exons that Changed Their Splicing Pattern after the Primate-Rodent Split</title> <p>We analyzed three datasets of human-mouse orthologous exons and their flanking introns and exons: 1) conserved constitutively spliced exons (constitutively spliced in both species, 45,553 exons), 2) conserved alternatively spliced exons (alternatively spliced in both species, 596 exons), and 3) exons that are alternatively spliced in human and constitutively spliced in mouse (species-specific alternative exons; 354 exons). Analysis of <italic>Alu</italic> insertions into introns flanking these exons revealed that species-specific alternative exons exhibited the highest level of <italic>Alu</italic> insertions, followed by conserved alternative exons; the group with the fewest intronic insertions were conserved constitutively spliced exons. We calculated the density of <italic>Alu</italic> insertions, namely the number of <italic>Alu</italic>s divided by the total intron length (and then multiplied by 1000 for convenience), in order to control for the fact that different intronic lengths might influence <italic>Alu</italic> insertion (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). On average, 0.42 <italic>Alu</italic> elements were found per 1000 bp within the upstream introns and 0.41 <italic>Alu</italic> elements were found per 1000 bp within the downstream introns of constitutively spliced exons; 0.49 and 0.44 <italic>Alu</italic>s per 1000 bp were found in the upstream and in the downstream introns of alternatively spliced exons, respectively; and 0.66 and 0.65 <italic>Alus</italic> were found per 1000 bp in the upstream and in the downstream introns of species-specific alternatively spliced exons, respectively. Thus, the number of <italic>Alu</italic> elements present in species-specific alternative exons differed significantly from that found in constitutively spliced exons (p-value = 7.16E-10, p-value = 5.22E-09, for upstream and downstream introns, respectively) and also differed from that in the alternatively spliced exons (p-value = 0.000503, p-value = 0.000014, for upstream and downstream introns, respectively).</p> <p>Furthermore, analysis of the distribution of antisense and sense <italic>Alu</italic>s upstream of conserved constitutively spliced exons revealed a selection against the presence of <italic>Alu</italic> elements adjacent to exons, specifically, against <italic>Alu</italic> elements in the antisense orientation (##FIG##0##Figure 1##). There are significantly fewer antisense <italic>Alu</italic>s compared to sense <italic>Alu</italic>s within 150 bp of constitutively spliced exons (p-value = 0.000012). Examination of the downstream intron did not reveal a significant bias (p-value = 0.056). This implies that a selective pressure exists against insertion of <italic>Alu</italic>s in close proximity upstream to constitutively spliced exons; this bias is stronger against <italic>Alu</italic>s in the antisense orientation than against the sense orientation.</p> <p>B1 is a rodent-specific retrotransposed element of ∼150 nucleotides that has the same ancestral origin as <italic>Alu</italic>: the 7SLRNA ##REF##17307271##[51]##,##REF##1774786##[52]##. Like <italic>Alu</italic>s, large numbers of B1 elements (a total of 331,015) reside within intronic sequences ##REF##17594509##[16]##. We found 236,036 B1 elements within 177,766 mouse introns that are found in GenBank (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). In mouse, 70,516 introns (39.6%) contain B1 elements. Overall, there are 1.32 B1 elements per intron. The median length of introns containing at least one B1 element is 3278 bp, whereas the median length of introns that do not contain B1 is 636 bp. These results indicate that <italic>Alu</italic> and B1 containing introns are substantially longer than other introns.</p> <p>We observed a significant correlation between the number of insertions of <italic>Alu</italic> and B1 elements within orthologous introns (Pearson correlation coefficient of 0.73 with p-value &lt;0.0001). Namely, orthologous introns show the same tendencies for <italic>Alu</italic> and B1 insertion, although these events happened independently after the split of the mouse and human lineages. We then set out to analyze whether insertion of B1 into rodent introns was biased in terms of location and orientation as was the case for <italic>Alu</italic> in primates. Analysis of B1 insertions within the flanking introns of conserved constitutively spliced exons, conserved alternatively spliced, and species-specific alternatively spliced exons yielded the same trend as that of <italic>Alu</italic> insertions in human. There was no statistical difference in the density of B1 between conserved constitutively spliced and conserved alternatively spliced exons; however, the upstream introns of the 258 species-specific exons (alternatively spliced in mouse, but constitutive in human) were significantly more enriched with B1 elements than were the upstream introns of conserved constitutively spliced exons (p-value = 0.0012) or constitutively spliced downstream introns (p-value = 0.014). This was also the case when the regions upstream of exons that are alternatively spliced in mouse and constitutively spliced in human were compared to the upstream introns of conserved alternatively spliced exons (p-value = 0.042) but not the downstream introns (p-value = 0.155). Therefore, insertion of retrotransposed elements into intronic sequences is correlated with the mode of splicing of the flanking exons.</p> <p>Five reports ##REF##12884004##[53]##–##REF##10737976##[57]## indicate that <italic>de novo Alu</italic> insertions into intronic sequences in antisense orientation and in close proximity to the affected exon (between 19–50 nucleotides) cause the downstream exon to shift from constitutive splicing to full exon skipping (three cases) or to alternative splicing (two cases) (##TAB##0##Table 1##). This effect of <italic>Alu</italic> elements on adjacent exons may be due to the <italic>Alu</italic> structure. <italic>Alu</italic> elements are comprised of two very similar segments, termed left and right arms. When an <italic>Alu</italic> is located in a gene in the antisense orientation and transcribed it contributes two poly-T stretches to the mRNA precursor. These poly-T regions might act as polypyrimidine tract (PPT) and, in combination with downstream 3′ and 5′ pseudo splice sites, might act as pseudo-exon ##REF##18276646##[58]##. Hence, such antisense <italic>Alu</italic>s that act as pseudo-exons might compete with nearby exons for the binding of splicing factors. These five cases of <italic>de novo Alu</italic> insertions imply that <italic>Alu</italic>s located in close proximity to exons might affect splicing of adjacent exons. This and the finding of <italic>de-novo Alu</italic> insertions that affect splicing imply that this is an on-going evolutionary process, which may result in novel transcripts that are deleterious and inflict genetic diseases. On the other hand, a shift in the splicing pattern from constitutive to alternative might be advantageous in some cases, and could enable testing new mRNA options without eliminating the old ones. Moreover, such a shift could introduce a premature termination codons enabling the expression of truncated proteins at certain needed times or in specific cell types and could be delicately regulated by the levels of splicing regulatory proteins ##REF##17361132##[59]##,##UREF##2##[60]##.</p> <p>In order to determine how many alternatively spliced exons are potentially regulated by the insertion of an antisense <italic>Alu</italic>, we used the alternative splicing track in the UCSC genome browser (##UREF##1##[48]## see also <xref ref-type="sec" rid="s4">Materials and Methods</xref>). In 269 events (∼1.5% out of 17,151 alternatively spliced cassette-exons), an antisense <italic>Alu</italic> was found within 100 bp upstream of an exon (150 have additional sense <italic>Alu</italic> within 2000 bp), 491 (∼2.8%) events in which an antisense <italic>Alu</italic> was found within 150 bp (273 have additional sense <italic>Alu</italic> within 2000 bp), and 689 (∼4%) events in which an antisense <italic>Alu</italic> was found within 200 bp (373 have additional sense <italic>Alu</italic> within 2000 bp). Out of these 689 alternative exons, 525 (76.1%) are conserved between human and mouse (23 events were recorded as alternatively spliced also in mouse alternative splicing track in version mm9). Within the human genome, almost 85% of alternative cassette exon skipping events are conserved in mouse, however only 76% of the cassette exon skipping events that have an adjacent <italic>Alu</italic> in opposite orientation are conserved within mouse genome. This is statistically significant, implying that there is a bias for <italic>Alu</italic> in antisense orientation in the regulation of alternative exons within non-conserved alternative splicing events (χ<sup>2</sup> test p-value = 1.6×10<sup>−8</sup>).</p> <title>Intronic <italic>Alu</italic>s Affect Splicing of Flanking Exons</title> <p>The above results suggest that stable insertion of <italic>Alu</italic>s into introns is associated with the mode of splicing of the flanking exons—especially the downstream exon—and that most <italic>Alu-Alu</italic> dsRNA is formed between sequences within the same intron. To test this hypothesis, exon 3 of the human <italic>RABL5</italic> gene was analyzed experimentally to examine the connection between intronic <italic>Alu</italic> and alternative splicing. A minigene containing exons 2 through 4 of the human <italic>RABL5</italic> gene (a conserved gene within all vertebrate genomes) was cloned. Exon 3 of <italic>RABL5</italic> is alternatively spliced in human and constitutively spliced in mouse, rat, dog, chicken, and zebrafish (see ##FIG##1##Figure 2## in ##REF##18020709##[17]##). Based on the phylogenetic relationships among the analyzed organisms, we conclude that the alternatively spliced variant is a derived form and the constitutively spliced variant is the ancestral one.</p> <p>Six <italic>Alu</italic>s have been inserted into the flanking introns of exon 3 since the last common ancestor of human and mouse: two in the upstream intron and four in the downstream intron (##FIG##1##Figure 2A##). <italic>Alu</italic>4, which is located in the downstream intron, resulted from insertion of an <italic>Alu</italic> within another <italic>Alu</italic> and will be regarded as one <italic>Alu</italic> (see also ##SUPPL##3##Text S1##). The minigene was transfected into 293T cells and the splicing products were examined following RNA extraction and RT-PCR analysis. Exon 3 in the <italic>RABL5</italic> minigene is alternatively spliced with approximately 40% inclusion (##FIG##1##Figure 2B##, lane 1). Removal of all intronic <italic>Alu</italic>s shifted splicing from alternative to constitutive (##FIG##1##Figure 2B##, compare lanes 1 and 2). This indicates that the insertion of <italic>Alu</italic>s into the flanking introns during primate evolution shifted exon 3 splicing from constitutive to alternative.</p> <p>Our experiments revealed that the orientation and position of the <italic>Alu</italic>s within the upstream intron affected splicing of exon 3. Deletion of the <italic>Alu</italic>s in the upstream intron, namely <italic>Alu</italic>1 and <italic>Alu</italic>2 (Δ1+2), had the same effect as deleting all <italic>Alu</italic>s (##FIG##1##Figure 2B##, lanes 7 and 2). The same effect was observed if one of the <italic>Alu</italic>s was deleted and the other was replaced with a 270-nucleotide non-<italic>Alu</italic> intronic sequence (##FIG##1##Figure 2B##, lanes 18 and 21, see also ##SUPPL##4##Text S2##). The replacement of each of the intronic <italic>Alu</italic>s with a non-<italic>Alu</italic> intronic sequence of a similar length eliminated the possibility that the effect observed after deletion related to shortening of the intron.</p> <p>The <italic>Alu</italic>s in the downstream intron, however, had a little or no effect on splicing (##FIG##1##Figure 2B##, lane 13). Taken together, it seems that the shift from constitutive to alternative splicing in the lineage leading to human is mediated mainly by <italic>Alu</italic>s 1 and 2. Deletion of <italic>Alu</italic>1 or replacement with a 270-nucleotide non-<italic>Alu</italic> intronic sequence resulted in almost complete exon skipping (##FIG##1##Figure 2B##, lanes 3 and 17, respectively). Replacement or deletion of <italic>Alu</italic>2 resulted in constitutive exon splicing (lanes 4 and 20).</p> <p>Interestingly, <italic>Alu</italic>1 and <italic>Alu</italic>2 have opposite effects on splicing. Deletion of both <italic>Alu</italic>s has the same effect as deleting only <italic>Alu</italic>2. Therefore, we concluded that <italic>Alu</italic>2 is dominant over <italic>Alu</italic>1. The dominance of <italic>Alu</italic>2 is also supported by two other observations. First, if all <italic>Alu</italic>s except <italic>Alu</italic>2 are removed, we observe almost total exon skipping (##FIG##1##Figure 2B##, lane 14). This indicates that <italic>Alu</italic>2 is a negative regulator of exon 3 recognition, unless <italic>Alu</italic>1, which is in opposite orientation to <italic>Alu</italic>2, is present (compare lanes 13 and 14). Furthermore, deletion or replacement with a 270-nucleotide non-<italic>Alu</italic> intronic sequence of <italic>Alu</italic>2 in combination with any additional intronic <italic>Alu</italic>s leads to constitutive splicing (lanes 7, 8, 10, 18–22). In the absence of <italic>Alu</italic>1 and the presence of <italic>Alu</italic>2, the dominance of <italic>Alu</italic>2 over the other <italic>Alu</italic>s is observed, leading to exon skipping (lanes 9, 12, 14, 15, and 16). As expected, in the presence of both <italic>Alu</italic>1 and <italic>Alu</italic>2, deletion of <italic>Alu</italic>s from the downstream intron had a marginal effect on splicing (##FIG##1##Figure 2B##, lanes 5, 6, 11, and 23).</p> <p>We demonstrated that the antisense orientation of <italic>Alu</italic>2 is essential for alternative splicing of exon 3. We first noted that the exact <italic>Alu</italic> family is not an important factor in determining splicing pattern: replacement of <italic>Alu</italic>1 of the Jo family with the sequence of <italic>Alu</italic>3 from the Sx family (both <italic>Alu</italic>s are in the sense orientation) did not affect the splicing pattern (##FIG##1##Figure 2B##, lane 24). Thus, the important factor is the presence of <italic>Alu</italic>1 in sense orientation. Our analysis showed that only when <italic>Alu</italic>2 is in the antisense orientation and <italic>Alu</italic>1 is in the sense orientation is alternative splicing of exon 3 observed (##FIG##1##Figure 2B##, lanes 25–26). These results indicate that the two <italic>Alu</italic>s in the upstream intron regulate alternative splicing of exon 3, whereas the three downstream intronic <italic>Alu</italic>s have no apparent effect on splicing of that exon. Moreover, <italic>Alu</italic>2 in the antisense orientation suppressed inclusion of exon 3, whereas <italic>Alu</italic>1 in the sense orientation antagonized the effect of <italic>Alu</italic>2.</p> <title>Formation of Double-Stranded RNA between the Two Intronic <italic>Alu</italic>s</title> <p>How do the two intronic <italic>Alu</italic>s regulate alternative splicing of exon 3? It is apparent that if <italic>Alu</italic>2 alone is present in the mRNA precursor, the exon is always skipped. We therefore postulated that in a population of the mRNA precursors that contain both <italic>Alu</italic>1 and <italic>Alu</italic>2, the two might form dsRNA formation (##FIG##2##Figure 3A##). This sequestration leads to exon inclusion; in the fraction of mRNA precursors with no base pairing between <italic>Alu</italic>1 and <italic>Alu</italic>2 the exon is skipped. To support this hypothesis, we examined whether RNA editing occurred in intron 2; editing would be indicative of formation of dsRNA. We searched the human EST/cDNA dataset and found five different mRNAs sequences containing intron 2 and comparison with genomic sequence indicated extensive RNA editing in both <italic>Alu</italic>1 and <italic>Alu</italic>2 (marked in red in ##FIG##2##Figure 3B##). The region of the editing was found to be in the middle of both <italic>Alu</italic>s. This suggests that these two <italic>Alu</italic> regions are in a double-stranded form.</p> <p>To further confirm the formation of dsRNA, we generated cDNA from a neuroblastoma cell-line using specific primers (##FIG##2##Figure 3C##). By using primers that hybridize in the exons flanking intron 2, we were able to observe a small population of mRNA molecules that contain intron 2; the majority of mRNAs are spliced (##FIG##2##Figure 3C## lane 1). We enriched the intron-containing fraction using primers that hybridize within the intron and within the downstream exon (##FIG##2##Figure 3C##, lane 2). Sequencing of the higher molecular weight PCR product using primer to <italic>Alu</italic>2 allowed us to identify four editing sites within it (##FIG##2##Figure 3D##).</p> <p>To confirm the importance of pairing between <italic>Alu</italic>1 and <italic>Alu</italic>2 on editing, we used the Δ<italic>Alu</italic>1 mutant (see ##FIG##1##Figure 2## lane 3) that led to a full exon skipping and examined the effect on editing within <italic>Alu</italic>2. There is one editing site within <italic>Alu</italic>2 that is dependent on the presence of <italic>Alu</italic>1; without <italic>Alu</italic>1 no editing at this site was observed (##FIG##2##Figure 3E##, the site is also highlighted in green in ##FIG##2##Figure 3B##). The other putative editing sites found in EST/cDNA show relatively low level of editing in the minigene.</p> <title>The Importance of the Distance between the Intronic <italic>Alu</italic>s and Exon 3</title> <p>We next set to examine if the distance between exon 3 and the intronic <italic>Alu</italic>s and the distance between <italic>Alu</italic>1 and <italic>Alu</italic>2 were important factors in splicing of exon 3. <italic>Alu</italic>2 is located 24 nucleotides upstream of exon 3. We identified the putative branch site of intron 2 and inserted an 800-nucleotide non-<italic>Alu</italic> intronic sequence upstream of the branch sequence and downstream of intronic <italic>Alu</italic>1 and <italic>Alu</italic>2 (marked B in ##FIG##3##Figure 4A##; see ##SUPPL##5##Text S3##). This insertion caused a shift from alternative to constitutive inclusion of exon 3 (##FIG##3##Figure 4B##, compare lane 1 and 2). Only when this insertion was shortened to less than 68 nucleotides did we begin to detect restoration of alternative splicing of exon 3; the level of skipping was further elevated when the inserted sequence was shortened to 56 or to 44 nucleotides (##FIG##3##Figure 4B##, lanes 3–9). To rule out the possibility that the sequence that was inserted contained intronic splicing regulatory sequences, we designed a fragment of 25 nucleotides free from known splicing regulatory sequences (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). We inserted this sequence into site B and also duplicated and triplicated this sequence to generate 50 and 75 nucleotides insertions. The longer is the inserted sequence, the higher is the inclusion level (##FIG##3##Figure 4B##, lanes 10–13). This indicates that the distance between the intronic <italic>Alu</italic>2 from exon 3 affects the mode of splicing.</p> <p>We also examined the effect of the distance between <italic>Alu</italic>1 and <italic>Alu</italic>2 on the splicing of exon 3. The insertion of the same fragment of 800 nucleotides between the two elements (marked A in ##FIG##3##Figure 4A##) led to substantial reduction in the inclusion level of exon 3, although alternative splicing was still observed (##FIG##3##Figure 4C##, compare lane 1 and lane 2). We then shortened this sequence, ultimately to 24 nucleotides; when the distance was shorter than 550 nucleotides, almost complete inclusion of exon 3 was observed (##FIG##3##Figure 4C## lanes 2–9). Our results indicate that the distance between <italic>Alu</italic>1 and <italic>Alu</italic>2 is important for maintaining the alternative splicing of exon 3; however, it is not as important as the distance between the intronic <italic>Alu</italic> elements and exon 3. We also note that increasing the distance between exon 3 and <italic>Alu</italic>2 leads to exon inclusion, whereas increasing the distance between <italic>Alu</italic>1 and <italic>Alu</italic>2 enhances exon skipping.</p> <title>
<italic>Alu</italic>2 Regulates Alternative Splicing of Exon 3</title> <p>
##FIG##1##Figure 2## demonstrates that <italic>Alu</italic>2 suppresses the inclusion of exon 3. We therefore analyzed the sequence of <italic>Alu</italic>2 to determine what regions might be critical for this effect. ##FIG##4##Figure 5A## shows the sequence of <italic>Alu</italic>2 and the mutations made. We deleted each of the two arms of <italic>Alu</italic>2 separately (##FIG##4##Figure 5B##, lanes 6 and 8) and mutated putative splicing signals (##FIG##4##Figure 5B##, lanes 2–5, 9–11). We found that the left arm of <italic>Alu</italic>2 is involved in the constitutive-to-alternative shift. Deletion of the left arm enhanced inclusion, whereas deletion of the right arm caused only a marginal effect (##FIG##4##Figure 5B##, compare lanes 1, 6, and 8). Mutations in the putative splicing signals in the right and left arms of <italic>Alu</italic>2 did not affect the splicing pattern (##FIG##4##Figure 5B##, lanes 2–5, 9–13). Our results do not support the possibility that the left arm functions as a pseudo-exon that abolishes or reduces selection of the exon 3 by competing with splicing factors (see also ##REF##18276646##[58]##,##REF##10958678##[61]##). However, analysis of this data is not straight forward, because deletion of the entire right arm together with the left-arm-polypyrimidine tract (LPPT), which produced a short Alu2 sequence, caused complete skipping of exon 3 (##FIG##4##Figure 5B##, lane 7). Insertion of a complementary sequence to the short <italic>Alu</italic>2 sequence along with its upstream intronic sequence (to complete an <italic>Alu</italic>-like length of 280 nucleotides), 105bp upstream to the original short <italic>Alu</italic>2 (mimicking the original distance of <italic>Alu</italic>2 from <italic>Alu</italic>1), either in the sense or antisense orientation, did not affect full skipping of exon 3. Deletion of the right arm alone or the LPPT alone had a marginal effect on splicing of exon 3 (##FIG##4##Figure 5B##, compare lanes 6 and 9 to 7). These results imply that multiple sequences along <italic>Alu</italic>2 combine to suppress the recognition of exon 3.</p> <p>Based on the location of the editing sites shown in ##FIG##2##figure 3B## we concluded that in a large part of the left arm of <italic>Alu</italic>2 there is no editing, suggesting that this part might not participate in a dsRNA structure. Within the left arm we identified a potential sequence, which is not part of the <italic>Alu</italic>1-<italic>Alu</italic>2 putative pairing alignment. This region has the potential to form an internal stem-loop structure (##FIG##4##Figure 5C##, upper alignments). Deletion of this region, replacement of this sequence with a similar stem-loop structure of a different sequence, creation of fully paired stem structure, or disruption of the stem structure all caused full exon 3 skipping. These results imply that the sequence, rather than the potential secondary structure, of this region is important for the inclusion of exon 3. This sequence contains two putative SC35 binding sites. However, mutations that eliminated these potential binding sites without generating another known splicing regulatory sites had no effect on splicing of exon 3, indicating that this is not the sequence involved in the regulation (##FIG##4##Figure 5C##, lane 2). There are no other potential binding sites for known splicing regulatory factors in this sequence (based on ##REF##16793546##[62]##–##REF##15607979##[64]##). Finally, addition of a potential complementary sequence to <italic>Alu</italic>1 did not effect splicing of exon 3 (not shown). Although <italic>Alu</italic>2 functions primarily to inhibit exon 3 selection, the above sequence within <italic>Alu</italic>2 enhances the inclusion of exon 3. Formation of a duplex between <italic>Alu</italic>1 and <italic>Alu</italic>2 is needed in order to present this intronic enhancer sequence properly for its effect on splicing of exon 3.</p>
<title>Discussion</title> <p>There are over 0.5 million copies of <italic>Alu</italic> elements in introns of human protein coding genes ##REF##17986453##[65]##, yet their function in regulation of gene expression is largely unknown. Here we show that intronic <italic>Alu</italic>s are not ‘neutral’ elements; they affect splicing of flanking exons. Some of these effects can be directly linked to the shift from constitutive to alternative splicing during primate evolution. The regulation demonstrated here involves both positive and negative effects of <italic>Alu</italic> element in antisense orientation, in close proximity, and upstream to the regulated exon. This complex regulation causes the downstream exon to shift from constitutive to alternative splicing. There are several examples ##REF##12884004##[53]##–##REF##10737976##[57]## of <italic>de novo</italic> insertions of <italic>Alu</italic> elements within introns that result in skipping of the adjacent exons. In three of the reported cases, the insertion of the <italic>Alu</italic> in the antisense orientation caused a total skipping of the adjacent exon.</p> <p>Exon 3 of <italic>RABL5</italic> gene, analyzed in this study, is alternatively spliced in human and constitutively spliced in mouse, rat, dog, chicken, and zebrafish. Six <italic>Alu</italic>s have been inserted into the flanking introns of exon 3 since the last common ancestor of human and mouse. <italic>Alu</italic>2 was inserted in the antisense orientation just upstream of exon 3 and functions as a negative element that suppresses exon 3 selection. This negative effect is partially reversed by another <italic>Alu</italic> present in the same intron in the sense orientation. Although we were not able to fully resolve the mechanism by which the two <italic>Alu</italic>s regulate alternative splicing of the downstream exon, we provide evidence that regulation requires the formation of a double-stranded region between the two <italic>Alu</italic>s and a combination of negative and positive sequences located in <italic>Alu</italic>2. The end result of this complex regulation is a shift from constitutive to alternative splicing of the downstream exon. This results in a new primate-specific mRNA isoform that could acquire novel functions, as well as maintaining the original mRNA. Moreover, such a shift could introduce a premature termination codon resulting in truncated proteins that might have regulatory roles at certain times or in specific cell types as could be delicately determined by the levels of splicing regulatory proteins ##REF##17361132##[59]##,##UREF##2##[60]##.</p> <p>Introns in humans are considerably longer than their mouse counterparts, mostly due to the presence of <italic>Alu</italic> elements ##REF##16005101##[66]##. Introns that flank alternatively spliced exons are longer than introns that flank constitutively spliced ones ##REF##17158149##[4]##,##REF##16260721##[67]##. We found a correlation between the splicing pattern of exons and the presence of <italic>Alu</italic>s in the flanking introns. First, alternatively spliced exons are flanked by introns containing more <italic>Alu</italic>s compared with introns flanking constitutively spliced ones, even when controlling for the difference in intron lengths. Second, more <italic>Alu</italic>s are present in human introns than are corresponding mouse B1 elements in the orthologous mouse introns. This second observation correlates with the finding that there are more species-specifically, alternatively spliced exons in human than in mouse (354/612 alternatively spliced events in human and 258/612 alternatively spliced events in mouse; χ<sup>2</sup>, p –value &lt;0.01).</p> <p>It was suggested that intron complementarities formed by multiple copies of <italic>Alu</italic> could help define and increase the splicing efficiency of very large metazoan introns ##REF##9356473##[68]##. However, it may be also possible that formation of a long and stable double-stranded structure in the upstream intron, especially near the splice site as in the case studied in this manuscript, reduces the ability of the splicing machinery to properly recognize the downstream exon, leading to slower splicing kinetics or suboptimal exon selection and, thus, to intron retention or exon skipping. Supporting the hypothesis that formation of dsRNA in introns might slow splicing is a recent publication showing that formation of dsRNA during pre-microRNA formation can slow splicing of the intron where the microRNA resides ##REF##17255951##[69]##.</p> <p>Although the effect of intronic retroelements on the splicing of flanking exons is presumably not a general trend that applies to all exons, it is relevant to a certain fraction of alternatively spliced exons (1.5% to 4% of the alternatively skipped exons in human).</p> <p>Our analysis indicated that the presence of <italic>Alu</italic> elements is correlated with the mode of splicing of adjacent exons. There is an ‘exclusion zone’ in intron sequences flanking exons, where insertion of <italic>Alu</italic> elements is presumably under purifying selection. The length of this ‘exclusion zone’ is similar to that of the human-mouse conserved sequences flanking alternatively spliced exons (∼80–150 nucleotides). This is presumably indicative of regions where the presence of intronic splicing regulatory sequences can affect alternative splicing of the adjacent exon ##REF##14992493##[5]##,##REF##17530917##[18]##,##REF##12840041##[70]##,##REF##16424921##[71]##. <italic>Alu</italic>s might be excluded from the proximal intronic sequences flanking constitutively spliced exons because <italic>Alu</italic>s were never inserted into these regions or because <italic>Alu</italic>s were inserted in an equal proportion in all gene regions (intronic and exonic) but we currently observe only those <italic>Alu</italic>s that have escaped purifying selection. The major burst of <italic>Alu</italic> retroposition took place 50–60 million years ago and has since dropped to a frequency of one new retroposition for every 20–125 new births ##UREF##3##[72]##,##REF##16522357##[73]##. As some of these insertions were deleterious and thus selected against, we probably detect intronic <italic>Alu</italic>s that are neutral, mildly deleterious, or beneficial to human fitness. Some of these beneficial intronic <italic>Alu</italic>s presumably altered splicing of the flanking exons and resulted in the generation of new isoforms that presented an advantage during primate evolution and were thus fixated in our genome. The research described here sheds light on how <italic>Alu</italic> elements have shaped the human genome.</p>
<p>Conceived and designed the experiments: GLM OR GA. Performed the experiments: GLM OR. Analyzed the data: GLM OR EK NS AG GA. Contributed reagents/materials/analysis tools: GLM GA. Wrote the paper: GLM EK NS EYL GA.</p> <p>Examination of the human transcriptome reveals higher levels of RNA editing than in any other organism tested to date. This is indicative of extensive double-stranded RNA (dsRNA) formation within the human transcriptome. Most of the editing sites are located in the primate-specific retrotransposed element called <italic>Alu</italic>. A large fraction of <italic>Alu</italic>s are found in intronic sequences, implying extensive <italic>Alu</italic>-<italic>Alu</italic> dsRNA formation in mRNA precursors. Yet, the effect of these intronic <italic>Alu</italic>s on splicing of the flanking exons is largely unknown. Here, we show that more <italic>Alu</italic>s flank alternatively spliced exons than constitutively spliced ones; this is especially notable for those exons that have changed their mode of splicing from constitutive to alternative during human evolution. This implies that <italic>Alu</italic> insertions may change the mode of splicing of the flanking exons. Indeed, we demonstrate experimentally that two <italic>Alu</italic> elements that were inserted into an intron in opposite orientation undergo base-pairing, as evident by RNA editing, and affect the splicing patterns of a downstream exon, shifting it from constitutive to alternative. Our results indicate the importance of intronic <italic>Alu</italic>s in influencing the splicing of flanking exons, further emphasizing the role of <italic>Alu</italic>s in shaping of the human transcriptome.</p> <title>Author Summary</title> <p>The human genome is crowded with over one million copies of primate-specific retrotransposed elements, termed <italic>Alu</italic>. A large fraction of <italic>Alu</italic> elements are located within intronic sequences. The human transcriptome undergoes extensive RNA editing (A-to-I), to higher levels than any other tested organism. RNA editing requires the formation of a double-stranded RNA structure in order to occur. Over 90% of the editing sites in the human transcriptome are found within <italic>Alu</italic> sequences. Thus, the high level of RNA editing is indicative of extensive secondary structure formation in mRNA precursors driven by intronic <italic>Alu</italic>-<italic>Alu</italic> base pairing. Splicing is a molecular mechanism in which introns are removed from an mRNA precursor and exons are ligated to form a mature mRNA. Here, we show that <italic>Alu</italic> insertions into introns can affect the splicing of the flanking exons. We experimentally demonstrate that two <italic>Alu</italic> elements that were inserted into the same intron in opposite orientation undergo base-pairing, and consequently shift the splicing pattern of the downstream exon from constitutive inclusion in all mature mRNA molecules to alternative skipping. This emphasizes the impact of <italic>Alu</italic> elements on the primate-specific transcriptome evolution, as such events can generate new isoforms that might acquire novel functions.</p>
<title>Supporting Information</title>
<fig id="pgen-1000204-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000204.g001</object-id><label>Figure 1</label><caption><title>Bioinformatic analysis of <italic>Alu</italic> insertion within the flanking introns.</title><p>Conserved constitutively spliced exons were analyzed for the differences in the location of antisense and sense <italic>Alu</italic>s within the upstream and downstream introns (left and right panels, respectively). The x-axis is the distance in base pairs from the exon; the y-axis is the number of <italic>Alu</italic>s found within this distance. Antisense <italic>Alu</italic>s are marked in blue and sense <italic>Alu</italic>s are marked in red.</p></caption></fig> <fig id="pgen-1000204-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000204.g002</object-id><label>Figure 2</label><caption><title>The effect of intronic <italic>Alu</italic>s on the splicing of a flanking exon.</title><p>(A) Schematic illustration of the <italic>RABL5</italic> minigene containing three exons and two introns. The intronic <italic>Alu</italic>s (1 through 5) are marked by boxes with a point indicating the orientation of the <italic>Alu</italic> relative to the pre-mRNA. <italic>Alu</italic>4 is an <italic>Alu</italic>Sx inserted between the two arms of <italic>Alu</italic>Jo. (B) The indicated wild-type (wt) and mutant plasmids were transfected into 293T cells, total RNA was extracted, and splicing products were separated on a 2% agarose gel after RT-PCR analysis. Lane 1, splicing products of wt <italic>RABL5</italic> minigene. Lanes 2–26, splicing products of the indicated mutants. The following abbreviations were used: Δ indicates deletion of the specified <italic>Alu</italic> element, X 1w3 indicates replacement of <italic>Alu</italic>1 with <italic>Alu</italic>3 (i.e., the sequence of <italic>Alu</italic>3 was inserted instead of that of <italic>Alu</italic>1 and in the same orientation as <italic>Alu</italic>3), and 1int specifies replacement of the <italic>Alu</italic>1 sequence with a non-<italic>Alu</italic> intronic fragment. The two mRNA isoforms are shown on the right. Numbers on top of the gel indicate percentage of exon inclusion as determined using ImageJ software. PCR products were sequenced.</p></caption></fig> <fig id="pgen-1000204-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000204.g003</object-id><label>Figure 3</label><caption><title>Editing sites within the intronic <italic>Alu</italic>s.</title><p>(A) Schematic illustration of exons 2 to 3 of the <italic>RABL5</italic> gene. Exons are depicted as black boxes; the intronic <italic>Alus</italic>, derived from <italic>Alu</italic>Jo and <italic>Alu</italic>Sx, in sense and antisense orientations, respectively, are shown in the middle gray-shaped boxes. The intronic, antisense <italic>Alu</italic> sequence (<italic>Alu</italic>Sx) is 102 nucleotides downstream of the sense <italic>Alu</italic>Jo and <italic>Alu</italic>Sx is 24 nucleotides upstream of the junction of exon 3. Sense and antisense <italic>Alu</italic>s are expected to form a double-stranded secondary structure, thus allowing RNA editing. (B) Editing sites were inferred from alignment of five cDNAs (accession numbers BC050531, BC038668, BI547904, BI548328, and DB495755) to the human genomic DNA. RNA editing occurs at eight positions within the antisense <italic>Alu</italic>Sx and at eleven positions within the sense <italic>Alu</italic>Jo. Based on these editing sites, the pairing between the sense and antisense <italic>Alu</italic> sequences was inferred (upper and lower lines, respectively). The region in which editing occurs starts at the middle of the right arm (position 232 in the <italic>AluSx</italic> consensus) and ends at the beginning of the left arm of <italic>Alu</italic>2 (position 101 in the <italic>AluSx</italic> consensus). Panel B shows only this corresponding region, while the entire <italic>Alu</italic>-<italic>Alu</italic> potential dsRNA is shown in ##SUPPL##2##Figure S3##. (C) To further confirm the editing activity, total RNA was extracted from a neuroblastoma (SH-SY5Y) cell-line and treated with DNaseI, followed by RT-PCR analysis using primers to exon 2 and exon 3, and to intron 2 and exon 3 (lanes 1 and 2, respectively; see also <xref ref-type="sec" rid="s4">Materials and Methods</xref>). The PCR products were sequenced. (D) The upper PCR product shown in panel C lane 2 was cloned and sequenced. The Chromas sequence is shown with the editing sites, found in <italic>Alu</italic>Sx, marked by boxes. (E) Editing in <italic>Alu</italic>2 requires the presence of <italic>Alu</italic>1. Wild type RABL5 minigene (WT) and a mutant in which <italic>Alu</italic>1 was deleted (Δ<italic>Alu</italic>Jo) were transfected into 293T cells. RNA was extracted, treated with DNase I, and amplified using set of primers flanking <italic>Alu</italic>2 and designed to amplify only exogenic transcripts. Sequencing chromatograms of four nucleotides in AluSx are shown (this editing site is also marked in green in panel B).</p></caption></fig> <fig id="pgen-1000204-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000204.g004</object-id><label>Figure 4</label><caption><title>Distance effect of <italic>Alu</italic> elements on the alternative splicing pattern.</title><p>(A) A schematic illustration of the genomic region between exons 2 and 3 of <italic>RABL5</italic> gene. Arrows marked A and B indicate two positions where an intronic sequence was inserted. (B) An 800-nucleotide intronic sequence was inserted in site B. The 800-nucleotide sequence was gradually shortened to the size shown above each lane. The indicated wt and chimeric plasmids were transfected into human 293T cells, total RNA was collected and examined by RT-PCR analysis (lanes 1–9). Lanes 10–12 show insertions of a different sequence, containing 25 nucleotides without any known splicing regulatory sequences, into the same site. This sequence was duplicated and triplicated to generate 50- and 75-nucleotide inserts. These mutant RABL5 minigenes were examined as above. (C) Similar analysis as in panel B, except that the 800-nucleotide sequence and the shorter sequences were inserted into site A. Spliced products are shown on the right and each PCR product was confirmed by sequencing. Splicing products were separated on a 2% agarose gel. Numbers on top of the gel indicate percentage of exon inclusion as determined using ImageJ software.</p></caption></fig> <fig id="pgen-1000204-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000204.g005</object-id><label>Figure 5</label><caption><title>The effect of <italic>Alu</italic>2 on the alternative exon.</title><p>(A) The sequence of the antisense <italic>Alu</italic>2. Mutated putative 5′ss is shown in red. A sequence of 24 nucleotides that was deleted is underlined. Three putative 3′ss that were mutated are in bold and underlined. In yellow are the right and left PPT regions with the downstream AG that were deleted. The green sequence is a stem and loop region of 18 nucleotides that was examined as shown in panel C (referred as ‘A’ region). Underlined in that region are two overlapping SC35 potential binding sites (the gray in the middle indicates the overlap region). (B) The indicated wt and mutant plasmids were transfected into 293T cells, total RNA was extracted, and splicing products were separated on a 2% agarose gel after RT-PCR analysis. Lane 1, splicing products of wt <italic>RABL5</italic> minigene. Lanes 2–13, splicing products of the indicated mutants. The two mRNA isoforms are shown on the right. Numbers on top of the gel indicate percentage of exon inclusion, as determined using ImageJ software. (C) The upper part illustrates the putative secondary structure formed by <italic>Alu</italic>Jo and <italic>Alu</italic>Sx, as predicted using the Vienna secondary structure web site (<ext-link ext-link-type="uri" xlink:href="http://www.tbi.univie.ac.at/RNA">http://www.tbi.univie.ac.at/RNA</ext-link>). The green arrows in the right panel indicate the start and end positions of the stem and loop structure, marked as ‘A’. The lower part shows the effect of the wt and mutant plasmids that were analyzed, as in panel B. Lane 1, wilt type. Lane 2, elimination of the two SC35 putative binding sites by mutating their overlapping sequence (AA were mutated to CC). Lane 3, deletion of the entire stem-loop ‘A’ sequence. Lane 4, replacement of the ‘A’ sequence by a random sequence lacking putative splicing binding sites. Lane 5, replacement of the loop sequence by a perfect complementary sequence, so that the ‘A’ element is entirely in a stem structure. Lane 6, disruption of the stem part of ‘A’, so that ‘A’ is entirely in a loop structure.</p></caption></fig>
<table-wrap id="pgen-1000204-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000204.t001</object-id><label>Table 1</label><caption><title>Diseases resulting from <italic>Alu</italic> insertion within an intron.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene</td><td align="left" rowspan="1" colspan="1">Disease</td><td align="left" rowspan="1" colspan="1">
<italic>Alu</italic>
</td><td align="left" rowspan="1" colspan="1">Intron insertion</td><td align="left" rowspan="1" colspan="1">Orientation</td><td align="left" rowspan="1" colspan="1">Distance from SS</td><td align="left" rowspan="1" colspan="1">Effect</td><td align="left" rowspan="1" colspan="1">Reference</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Fas</td><td align="left" rowspan="1" colspan="1">ALPS syndrome</td><td align="left" rowspan="1" colspan="1">Sb1</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">antisense</td><td align="left" rowspan="1" colspan="1">50 bp upstream of exon 8</td><td align="left" rowspan="1" colspan="1">Skipping of exon 8</td><td align="left" rowspan="1" colspan="1">Tighe et al., 2002 ##REF##12215906##[55]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">GK</td><td align="left" rowspan="1" colspan="1">Glycerol kinase deficiency</td><td align="left" rowspan="1" colspan="1">Y</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">antisense</td><td align="left" rowspan="1" colspan="1">52 bp upstream of exon 5</td><td align="left" rowspan="1" colspan="1">Alternative Splicing</td><td align="left" rowspan="1" colspan="1">Zhang et al., 2000 ##REF##10737976##[57]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">FGFR2</td><td align="left" rowspan="1" colspan="1">Apert syndrome</td><td align="left" rowspan="1" colspan="1">Ya5</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">antisense</td><td align="left" rowspan="1" colspan="1">19 bp upstream of exon 9</td><td align="left" rowspan="1" colspan="1">Alternative Splicing</td><td align="left" rowspan="1" colspan="1">Oldridge et al., 1999 ##REF##9973282##[54]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">NF1</td><td align="left" rowspan="1" colspan="1">Neurofibromatosis type1</td><td align="left" rowspan="1" colspan="1">Ya5</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">antisense</td><td align="left" rowspan="1" colspan="1">44 bp upstream of exon 6</td><td align="left" rowspan="1" colspan="1">Skipping of exon 6</td><td align="left" rowspan="1" colspan="1">Wallace et al., 1991 ##REF##1719426##[56]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">Factor VIII</td><td align="left" rowspan="1" colspan="1">Hemophilia A X-linked severe bleeding disorder</td><td align="left" rowspan="1" colspan="1">Yb9</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">antisense</td><td align="left" rowspan="1" colspan="1">19 bp upstream of exon 19</td><td align="left" rowspan="1" colspan="1">Skipping of exon 19</td><td align="left" rowspan="1" colspan="1">Ganguly et al., 2003 ##REF##12884004##[53]##
</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula></disp-formula>
<supplementary-material content-type="local-data" id="pgen.1000204.s001"><label>Figure S1</label><caption><p>Intron length distribution of human introns.</p><p>(0.17 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000204.s002"><label>Figure S2</label><caption><p>A screen shot created by the UCSC genome browser.</p><p>(1.44 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000204.s003"><label>Figure S3</label><caption><p>Potential dsRNA of AluJo+ and AluSx-.</p><p>(0.51 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000204.s004"><label>Text S1</label><caption><p>Minigenes' sequence.</p><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000204.s005"><label>Text S2</label><caption><p>270 nucleotides non-<italic>Alu</italic> intronic sequence from intron 20 of IKBKAP gene.</p><p>(0.02 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000204.s006"><label>Text S3</label><caption><p>800 nucleotides intronic sequence from intron 11 of IMP gene.</p><p>(0.02 MB DOC)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work was supported by the Israel Science Foundation (1449/04 and 40/05), GIF, ICA (through the Ber-Lehmsdorf Memorial Fund), DIP, and EURASNET. EK is a fellow of the Clore Scholars Programme. AG is supported by the Adams Fellowship Program of the Israel Academy of Sciences and Humanities.</p></fn></fn-group>
<graphic xlink:href="pgen.1000204.g001"/> <graphic id="pgen-1000204-t001-1" xlink:href="pgen.1000204.t001"/> <graphic xlink:href="pgen.1000204.g002"/> <graphic xlink:href="pgen.1000204.g003"/> <graphic xlink:href="pgen.1000204.g004"/> <graphic xlink:href="pgen.1000204.g005"/> <graphic xlink:href="pgen.1000204.e001.jpg" mimetype="image" position="float"/>
<media xlink:href="pgen.1000204.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000204.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000204.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000204.s004.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000204.s005.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000204.s006.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["29"], "element-citation": ["\n"], "surname": ["Jepson", "Reenan"], "given-names": ["JE", "RA"], "year": ["2007"], "article-title": ["RNA editing in regulating gene expression in the brain."], "source": ["Biochim Biophys Acta"]}, {"label": ["48"], "element-citation": ["\n"], "surname": ["Karolchik", "Kuhn", "Baertsch", "Barber", "Clawson"], "given-names": ["D", "RM", "R", "GP", "H"], "year": ["2007"], "article-title": ["The UCSC Genome Browser Database: 2008 update."], "source": ["Nucleic Acids Res"]}, {"label": ["60"], "element-citation": ["\n"], "surname": ["McGlincy", "Smith"], "given-names": ["NJ", "CW"], "year": ["2008"], "article-title": ["Alternative splicing resulting in nonsense-mediated mRNA decay: what is the meaning of nonsense?"], "source": ["Trends Biochem Sci"]}, {"label": ["72"], "element-citation": ["\n"], "surname": ["Deininger", "Batzer"], "given-names": ["P", "MA"], "year": ["1993"], "article-title": ["Evolution of retroposons."], "source": ["Evol Biol"], "volume": ["27"], "fpage": ["157"], "lpage": ["196"]}, {"label": ["76"], "element-citation": ["\n"], "surname": ["Smit", "Hubley", "Green"], "given-names": ["AF", "R", "P"], "year": ["1996"], "article-title": ["RepeatMasker Open-3.0"]}]
PLoS Genet. 2008 Sep 26; 4(9):e1000204