# Detection of Introns in Eukaryotic Small Subunit Ribosomal RNA Gene Sequences

Authors: Dipankar Bachar; Laure Guillou; Richard Christen
Journal: Dataset Papers in Biology
 (2013)

---

## Abstract

The gene encoding SSU-rRNA sequences is the tool of choice for phylogenetic analyses and environmental biodiversity analyses of bacteria, Archaea but also unicellular Eukaryota. In Eukaryota, gene sequences may often be interrupted by long or several introns. Searching in GenBank release 188, we found descriptions of 3638 such sequences. Using a database of 180 000 SSU-rRNA sequences well annotated for taxonomy and a C++ program written for that purpose, we computed the presence of 18 691 introns (among which the 3638 described introns). Filtering on length and sequence quality, 3646 sequences were retained. These introns were clustered; clusters were analyzed for the presence of single or multiple clades at various levels of taxonomic depth, allowing future analyses of horizontal transfers. Various analyses of the results are provided as tabulated files as well as FASTA files of described or computed introns. Each sequence is annotated for cellular location (nuclear, chloroplast, and mitochondria), positions at which they were found in the SSU-rRNA sequences and taxonomy as provided by GenBank.


---

## Body

## 1. Introduction

The gene that microbiologists use to determine the taxonomic affiliations of microbes using molecular methods needs to meet a number of requirements. It has to be conservative in its function and present in every organism analyzed. Often the presence of conserved domains is required, allowing the design and use of universal PCR primers. Finally, sequences of most known organisms must be available in the public databases (i.e., International Nucleotide Sequence Database Collaboration (INSDC) between Japanese, European, and American nucleotide databases, resp., DDBJ, ENA, and GenBank,http://www.insdc.org/). A gene meets these requirements to a high degree: the gene for one of the RNA subunits that together form the ribosome, also known as the small subunit ribosomal RNA (SSU-rRNA) gene. For that reason, the gene encoding the SSU-rRNA serves as a prominent tool for phylogenetic and environmental biodiversity analyses of bacteria, Archaea but also unicellular Eukaryota [1–3].SSU-rRNA gene sequences may contain numerous self-splicing introns of variable lengths [4–23]. The SSU-rRNA genes can thus be enlarged to up to 3.5 kb. Introns have rarely been identified in bacterial SSU-rRNA gene sequences (see one example in Thiomargarita namibiensis [5]), but they are often present in SSU-rRNA gene sequences of Eukaryota (see the aforementioned part). Such length heterogeneity of SSU-rRNA gene sequences has so far seldomly been considered when constructing SSU-rRNA-based clone libraries or more recently for deep sequencing of SSU-rRNA amplicons. The detection of elongated SSU-rRNA gene sequences may have implications for analyses in molecular ecology and may cause systematic biases. When dealing with clone libraries, the usually low number of sequences allows for a manual analysis, and introns in new sequences can be easily detected and not taken into account. However, with the advent of Next Generation Sequencing (NGS) technologies, hundreds of thousands of sequences are generated during each experiment, and these data need to be analyzed via automated pipelines. If a sequence from a given organism differs from the sequence of the same organism present in the database via intron insertions/deletions, proper identification will frequently fail. Also, NGS pipelines apply filters on amplicon lengths, and a bias due to the PCR amplification step will discriminate against longer sequences. As a result, some groups of organisms may be simply discarded because the domains amplified become too long due to the insertion of introns.In this paper, we assessed the frequency of such introns using the release 188 of GenBank, for every sequence that contained at least one domain identified as a partial or complete SSU-rRNA sequence of eukaryotic origin. We also identified the origin of these sequences as nuclear, mitochondria, or chloroplast encoded. In a first step, we selected and extracted subsequences that had been annotated to contain introns. We then used these sequences to identify and extract putative introns from the entire set of SSU-rRNA gene sequences of eukaryotic origin. Finally, for each category, we provide the longest described intronic sequence. This dataset can be used to feed automated pipelines dealing with NGS datasets.


## 2. Methodology

A database of sequences was built as follows. Reference sequences were retrieved from GenBank release 188 (152 177 733 sequences) using a series of 175 keywords. 674 146–688 113 subsequences were retrieved using ACNUC [24]. Filtering out subsequences shorter than 800 nt and longer than 8000 nt provided 237 967 sequences. Filtering for Eukaryotic origin provided 145 968 sequences. These sequences were contained in 145 163 GenBank files which were used for the analyses of SSU-RNA gene sequences. Introns were retrieved from the same query but without filtering on sequence length, that is, using 471 263 GenBank files (entries) (598 851 rRNA subsequences).A given entry may contain a single SSU rRNA gene sequences, sequences of several genes, or a complete genome. When extracting SSU rRNA gene sequences and looking, we had for a long time considered two different cases.(i)
The rRNA is solely described as to be joined, with no description of introns as in accession number EF689886  (939 occurrences found):rRNA
join(<1..213,559..660,712..146,1790..2004,
2394..2407,2721..2904,3130..>3246)
/product =“18S ribosomal RNA”(ii)
Both join and introns are described as in accession number U54637 (2187 occurrences found):rRNA
join(650..2046,3427..3926)
/product = “small subunit ribosomal RNA”
intron
2047..3426
/note = “small subunit rRNA; group-IC2”This led us to build a database of SSU-rRNA annotated in the following way:(i)
_U|: sequence simply annotated as SSU-rRNA;
(ii)
_R|: sequence described by a join (see the following) and the corresponding sequence of that of rRNA (a consequence of the older default case for protein coding sequences);
(iii)
_G|: the genomic sequences (rDNA gene) corresponding to the _R| sequences including putative introns.During the time course of this study, we discovered several other possibilities:(1)
a simple description of the presence of intron(s) such as that of accession number AM231292, 512 occurrences, but the join to describe the rRNA sequence is missing:gene
<1..>3533
/gene = “18S rRNA”
rRNA
<1..>3533
/gene = “18S rRNA”
/product = “18S ribosomal RNA”
intron
569..1716
/gene = “18S rRNA”
/note = “S529 group I intron”
intron
2837..3365
/gene = “18S rRNA”
/note = “S1389 group I intron”(2)
some strange annotations such as those of  accession number X80710, for which exons are described, but within introns of length 0:rRNA
join(1..528,529..1134,1135..1743)
/product = “18S ribosomal RNA”
exon
1..528
/number = 1
exon
529..1134
/number = 2
exon
1135..1743
/number = 3(3)
still another strange case such as that of accession number AB016009, which creates two separate sequences for the SSU rRNA gene:rRNA
<1..1618
/product = “18S ribosomal RNA”
intron
1619..2115
/note = “Group IC1 intron (site 1389);
Insertion site is site 1389 corresponding toE. coli 16S
rDNA"
rRNA
2116..>2254
/product = “18S ribosomal RNA”(4)
note finally peculiar cases such as that of accession number AF274109, in which introns are described as spliced sequences:source
1..2126
/organism = “Ochrolechia parella”
/mol_type = “genomic DNA”
/db_xref = “taxon:129506”
rRNA
join(<1..384,443..1469, 1825..>2126)
/product = “18S ribosomal RNA”
intron
join(385..442,1470..1824)(5)
the creation of chimeric sequence accession number AB259428, because the sequences extracted are a concatenation of 5′ and 3′ ends of the molecule, which could suggest that an intron has been found in between, when this is only a series ofN inserted by the authors:source
1..525
/organism = “Didymium panniforme”
/mol_type = “genomic DNA”
/specimen_voucher = “TNS-M-Y-16880 (TNS)”
/db_xref=“taxon:385452”
rRNA
join(<1..233,334..>525)
/product = “small subunit ribosomal RNA”
gap
234..333
/estimated_length = unknown(6)
finally some sequences such as AJ506970 which are quite difficult to treat:source
1..1906
/organism = “Physcia dimidiata”
/mol_type = “genomic DNA”
/specimen_voucher = “Mayrhofer & Litterski 13.932 (GZU)”
/db_xref = “taxon:116814”
gene
1..1906
/gene = “nrSSU”
gene
order(<1..507,825..1262,1476..1906,AJ507615.1:1..12, AJ507615.1:206..565)
/gene = “nrSSU”
rRNA
join(<1..507,825..1262,1476..1906,AJ507615.1:1..12,AJ507615.1:206..565)
/gene = “nrSSU”
/product = “18S ribosomal RNA”
exon
1..507
/gene = “nrSSU”
/number = 1
intron
508..824
/gene = “nrSSU”
/number = 1
exon
825..1262
/gene = “nrSSU”
/number = 2
intron
1263..1475
/gene = “nrSSU”
/number = 2
exon
1476..1906
/gene=“nrSSU”
/number = 3Annotation extends the sequence over two or more accession numbers. In our present annotation and extraction pipeline, the _U| sequence will cover the extension over the various accession numbers, while the _G| sequence will be extracted solely from the sequence contained in the primary accession number (the one contained in the entry).Each of these different annotations was used to build a database on described intronic sequences.Reference taxonomy was built as follows. rRNA sequences of 800 nt or more were used to build a taxonomic database, using manual analyses and with homogenization of the taxonomy to 8 fields from domain (nuclear or organelle), kingdom, phylum, and so forth to genus and finally species. This database was then used to estimate the origins of introns retrieved. Note however that in most files provided in the dataset, the original GenBank taxonomy was kept.Introns were computed as follows. We developed a special C++ program that allows to search a file of sequences within another file of sequences, similarly to Blast, except that a Needleman-Wunsch global alignment is used. A search is done with the following parameters:(i)
n which is the maximal number of best hits to return,
(ii)
p which is the minimum percentage of similarity to return a hit,
(iii)
similarity which is computed from the alignment, not taking into account terminal gaps.After different trials, parametersn
=
30 and p
=
90 were retained. Every SSU-rRNA sequence was searched for the presence of introns; for each sequence introns found were sorted by decreasing length and successively removed from the sequence when present. When overlapping introns are found, only the longest one is thus taken as a true intron. We are aware that this approach is still crude, but in absence of inspection of the present results by experts, it is difficult to improve the algorithm without sound proofs of false positive and false negative results. Checking the algorithm on protein-coding sequences could be an alternate solution. However checking that was done on sequences that had been described as containing introns showed the results to be good (see the results). These results were then analyzed using a series of dedicated Python scripts that produced the dataset provided.598 851 entries were retrieved from GenBank release 188 (April 2012), which contained annotations identifying at least one subsequence as being an SSU-rRNA gene sequence. In this paper, the termentry will be used to refer to every information available for a given accession number. A given entry may contain a single SSU-rRNA gene sequence, sequences of several genes, or a complete genome. When extracting SSU-rRNA gene sequences, we considered three different cases depending on how SSU-rRNA sequences are annotated for introns, and each sequence was labeled as follows:(i)
_U|: for sequences simply annotated as SSU-rRNA,
(ii)
_R|: for sequences described by a join (implicit description), which implies that the gene contains intron(s), or with introns described in the features (explicit description),
(iii)
_G|: for the genomic sequences (rDNA gene sequence) corresponding to the _R| sequence including putative introns.We thus built a database of 184 827 sequences of 800 nt or more (_U|: 182 107 sequences; _G| and _R|: 1360 sequences, some of which may contain several introns). Extraction of introns—either explicitly or implicitly—described led us to build a database of 3638 intronic sequences for intronic sequences of 40 nt or more and not having at least 10 consecutiveN (to which we then added some more sequences as described previously; total sequences: 3646); see Table 1.Table 1
(a)
Annotations used to extract intronic sequences.
Annotation
Number
Sequences with both join and introns
2187
Sequences with only introns
512
Sequences with only join
939
(b)
Genomic localization of described introns.
Annotation
Number
Mitochondria encoded
50
Nuclear encoded
3560
Chloroplast encoded
36We also extracted the annotation that described the type of each intronic sequence (see Table2).Table 2
Category of well-annotated intronic sequences.
Intron Type
Number of Sequences
Group I intron
1461
Group IA2 intron
1
Group IA3 intron
1
Group IC intron
1
Group IC1 intron
42
Group IC2 intron
3
Group IE intron
5
Group II intron
10
Group IIA1 intron
8
Group IIB1 intron
10As shown in Table2, the vast majority of introns are described as group I introns. We used the Vienna RNA Web-Servers (http://rna.tbi.univie.ac.at/) to visualize predictions of secondary structures of these sequences, but we were not able to find any obvious conserved secondary structure among a given category (group I or group II), but we did not use any extensive search as described, for example, in [25].We clustered these sequences as described previously and looked at the clusters obtained. Figure1 shows an example of the result.Figure 1
Comparisons of 5′ and 3′ boundaries of sequences submitted by different authors. AY766358 was submitted in 2007 by Milstein and Oliveira as “group I intron.” AY573851 was submitted in 2004 by Teasdale, Klein, West, and Mathieson as “group ICI intron.” Remaining sequences were submitted in 2002 by Broom et al. as “putative group IC1 intron.”The results shown in Figure1 are a typical example of what could be seen: in a given cluster, some sequences differ for the boundaries of the intronic sequences. It is remarkable that for sequences submitted by the same author(s), boundaries are identical but differ from the boundaries found by different authors. We can also observe that the group to which an intron is assigned can be more or less precise.Described introns are predominantly found within the 50–600 nt range. Only 5 introns had a length of more than 2000 nt (Table3).Table 3
Longest introns described (>2000 nt). CG: complete genome; rRNA: only rRNA(s); fragment: long genomic fragment; 1: join only; 2: intron only; 3: both.
Accession Number and Genomic Locations
Cell Location
Clade
Species
Genome Fragment
Description of Intron
M68929.153546.157925
mito-SSU
Viridiplantae
Marchantia polymorpha
CG
3
HQ292070.249.3107
mito-SSU
Fungi
Ophiostoma minus
rRNA
3
HQ292073.249.3106
mito-SSU
Fungi
Ophiostoma minus
rRNA
3
HQ292072.249.3106
mito-SSU
Fungi
Ophiostoma minus
rRNA
3
HQ343317.195.3047
Nuclear
Fungi
Ophiostoma torulosum
rRNA
3
X91267.36.2342
Nuclear
Viridiplantae
Scenedesmus pupukensis
rRNA
3
AF029891.2413.4594
mito-SSU
Fungi
Cryphonectria parasitica
fragment
3
EF689875.1460.3551
Nuclear
Fungi
Verrucaria marmorea
rRNA
1The following describes the computation for retrieval of undescribed introns. We next used these 3638 sequences to analyze our database of 184 827 SSU-rRNA sequences with a length of 800 nt or more for the presence of such intronic sequences. Since some of the _G| (annotated in GenBank for the presence of introns) could not be retrieved, we manually added back some sequences that had been filtered out because they either were too short or contained too manyN; this provided a new database of 3638 intronic sequences that we again used to find introns in rRNA-SSU sequences. Numbers found are described in Table 4.Table 4
Type of sequence and number of introns found. _UC: introns found in sequences that had no description of introns; _GC: introns found in genomic sequences, described for presence of introns; _RC: introns found in sequences after removal of introns as described in annotations.
Sequence Type Found
Number of Introns
_UC
15 822
_GC
2162
_RC
607Introns found in _G| sequences are interesting to test the efficiency of our algorithm, since we would expect for each _G| sequence to retrieve introns as described in the GenBank annotations (usually and for finding introns, the authors compared the genomic sequence to the rRNA sequence). Results of this analysis were as follows: in 1555 cases, we exactly retrieved the boundaries described in the entries. However, in 357 cases, there was a discrepancy. Dataset Item 29 (Table) shows what occurred. In most cases, the difference lied in one or a few positions at either or both ends. This problem results from the extreme difficulty to exactly determine intron boundaries in non-protein-coding sequences. We looked for a confirmation of this by clustering annotated introns allowing from 5 to 20 differences between sequences (taking or not external gaps as differences). Looking at aligned sequences confirmed the results of Table1 (see Figure 2).Figure 2
Distribution length of described SSU-rRNA introns.The _RC sequences are puzzling, since we did not expect to retrieve many “new” introns from rRNA sequences which had been carefully studied for the presence of introns. We randomly and manually analyzed a small number of cases, but we were not able to decide whether or not these “new” introns were real cases or artifacts. The main problem lies in the methods section of the publications associated with the sequences. In no case, we were able to find a proper description of how introns had been identified, that is, either by comparisons with public sequences of closely related organisms or by comparison between the genomic sequence and the direct sequencing of the rRNA sequence. Since rRNA sequencing is tedious, we expect the first case to be true in most cases.Since new intronic sequences are differently annotated according to the kind of sequences they extracted, a conservative approach would be not to take into account the _RC sequences as real introns. On the other hand, keeping every intronic sequence would allow identification of yet unknown introns in the genomes of a given clade (Table5).Table 5
Taxonomy and location of described introns.
Number of Introns
Cell Location
Clade
33
chloro-SSU
Viridiplantae
2
mito-LSU
Fungi
7
mito-LSU
Viridiplantae
1
mito-LSU
Stramenopiles
40
mito-SSU
Fungi
1
mito-SSU
Haptophyceae
1
mito-SSU
Ichthyosporea
7
mito-SSU
Viridiplantae
9
Nuclear
Alveolata
253
Nuclear
Amoebozoa
1
Nuclear
Euglenozoa
1928
Nuclear
Fungi
24
Nuclear
Heterolobosea
1
Nuclear
Metazoa
10
Nuclear
Rhizaria
338
Nuclear
Rhodophyta
968
Nuclear
Viridiplantae
13
Nuclear
StramenopilesWe then clustered the introns found and analyzed clusters of 10 or more sequences which we extracted as FASTA sequences. Differences observed at boundaries confirmed our previous observations and that finding boundaries by computation alone is even more difficult both because known introns may differ slightly at boundaries and because computation is difficult, even using a global alignment algorithm.The problem appears more acute on the 3′ end, most likely due to a not accurate enough detection of the intron endings (Figure3).Comparisons of 5′ (a) and 3′ boundaries (b) of intron sequences extracted by our pipeline and grouped in the same cluster.
(a)
(b)Our final analysis was to examine each cluster for taxonomic discrepancies. Because horizontal transfers of rRNA introns seem to be rare events [19, 26, 27], we expected that the taxonomy would be homogeneous within each cluster, at least down to the family level. On the other hand, the presence of a given intron in distinct clades would suggest an old origin of this intron, followed by vertical transfer. Because the taxonomy associated with unicellular organism and provided by GenBank is often very poor, for sequences retrieved in various clusters we developed our own taxonomy, based on manual analyses (see previous explanation). These sequences were all annotated using 8 successive fields ranging from the domain (Eukaryota or Organelle) to the species level.We provide a series of FASTA files in which only sequences extracted from abnormal clusters are shown; each sequence is annotated with our taxonomy. A series of files is thus provided, each with sequences clustered at a given level of dissimilarities and for discrepancies found at successive levels of the taxonomy. Note however that some clusters can be found abnormal when they are probably not. This is true in particular when a discrepancy is found as, for example, between “Eukaryota Opisthokonta” and “Eukaryota Eukaryota_X”, the latter annotation resulting from our failing to properly assign an SSU-rRNA sequence. When such case is encountered, it can probably be manually solved.


## 3. Dataset Description

The dataset associated with this Dataset Paper consists of 69 items which are described as follows.1
The longest sequence in each group of introns (groups I, IA2, IA3, IC, IC1, IC2, IE, II, IIA1, and IIB1).2
Introns described in GenBank files (longer than 40 nt). Most introns containing long stretches ofN have been removed (some were kept because they had unique sequences). Each sequence is described as follows: >accession.p1.p2|type|taxonomy, where accession is the accession number, p1 is the 5′ position of the introns in the entries of sequences, p2 is the 5′ position of the introns in the entries of sequences, type describes the intron origin, and taxonomy is the taxonomy provided by GenBank (| separated fields).3
Introns described in GenBank files (longer than 40 nt and sorted by increasing length).4
Clusters of known introns at several dissimilarities (5 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.5
Clusters of known introns at several dissimilarities (10 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.6
Clusters of known introns at several dissimilarities (15 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.7
Clusters of known introns at several dissimilarities (20 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.8
Clusters of known introns at several dissimilarities (25 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.9
Clusters of known introns at several dissimilarities (30 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.10
Clusters of known introns at several dissimilarities (5 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences.11
Clusters of known introns at several dissimilarities (10 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences.12
Clusters of known introns at several dissimilarities (15 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences.13
Clusters of known introns at several dissimilarities (20 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences.14
Clusters of known introns at several dissimilarities (25 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences.15
Clusters of known introns at several dissimilarities (30 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences.16
Sequences for clusters of known introns containing 10 or more sequences of clusters of 5 differencesnottaking into account terminal gaps as differences. Each cluster is separated by an empty sequence.17
Sequences for clusters of known introns containing 10 or more sequences of clusters of 10 differencesnottaking into account terminal gaps as differences. Each cluster is separated by an empty sequence.18
Sequences for clusters of known introns containing 10 or more sequences of clusters of 15 differencesnottaking into account terminal gaps as differences. Each cluster is separated by an empty sequence.19
Sequences for clusters of known introns containing 10 or more sequences of clusters of 20 differencesnottaking into account terminal gaps as differences. Each cluster is separated by an empty sequence.20
Sequences for clusters of known introns containing 10 or more sequences of clusters of 25 differencesnottaking into account terminal gaps as differences. Each cluster is separated by an empty sequence.21
Sequences for clusters of known introns containing 10 or more sequences of clusters of 30 differencesnottaking into account terminal gaps as differences. Each cluster is separated by an empty sequence.22
Sequences for clusters of known introns containing at least 3 sequences at several dissimilarities (5 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 3 sequences are shown.23
Sequences for clusters of known introns containing at least 3 sequences at several dissimilarities (10 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 3 sequences are shown.24
Sequences for clusters of known introns containing at least 3 sequences at several dissimilarities (15 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 3 sequences are shown.25
Sequences for clusters of known introns containing at least 3 sequences at several dissimilarities (20 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 3 sequences are shown.26
Sequences for clusters of known introns containing at least 3 sequences at several dissimilarities (25 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 3 sequences are shown.27
Sequences for clusters of known introns containing at least 3 sequences at several dissimilarities (30 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 3 sequences are shown.28
Variations in sequence length of described introns.
0
3
0
Minimum Length of Sequence
Maximum Length of Sequence
Number of Sequences29
Differences between boundaries described for introns (in GenBank entries) and boundaries found by computation.
0
3
0
Described Introns
Found 5′
Found 3′30
Tabulated file for computed introns with description and sequence in which an intron was found with the following fields: ID of sequence in which an intron was found by computation, ID of sequence in which this intron had been described (accession.p1.p2: p1, p2 positions of introns in accession sequence), and sequences of introns found.31
Clusters of computed introns at several dissimilarities (5 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.32
Clusters of computed introns at several dissimilarities (10 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.33
Clusters of computed introns at several dissimilarities (15 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.34
Clusters of computed introns at several dissimilarities (20 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.35
Clusters of computed introns at several dissimilarities (25 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.36
Clusters of computed introns at several dissimilarities (30 differences between the longest sequences and other sequences), taking into account terminal gaps as differences.37
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (5 differences between the longest sequences and other sequences), taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.38
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (10 differences between the longest sequences and other sequences), taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.39
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (15 differences between the longest sequences and other sequences), taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.40
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (20 differences between the longest sequences and other sequences), taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.41
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (25 differences between the longest sequences and other sequences), taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.42
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (30 differences between the longest sequences and other sequences), taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.43
Clusters of computed introns at several dissimilarities (5 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences, but clustering not counting end gaps as differences.44
Clusters of computed introns at several dissimilarities (10 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences, but clustering not counting end gaps as differences.45
Clusters of computed introns at several dissimilarities (15 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences, but clustering not counting end gaps as differences.46
Clusters of computed introns at several dissimilarities (20 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences, but clustering not counting end gaps as differences.47
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (5 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.48
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (10 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.49
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (15 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.50
Sequences for clusters of computed introns containing at least 5 sequences at several dissimilarities (20 differences between the longest sequences and other sequences),not taking into account terminal gaps as differences. Each cluster is separated by an empty sequence. Only clusters with at least 5 sequences are shown.51
Introns found in genomes of chloroplasts.52
Introns found in genomes of mitochondria.53
Introns found in nuclear genomes.54
Introns found in nuclear genomes. Sequences are sorted by decreasing length.55
Sequences extracted from clusters that showed 5 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 2. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.56
Sequences extracted from clusters that showed 5 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 4. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.57
Sequences extracted from clusters that showed 5 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 5. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.58
Sequences extracted from clusters that showed 5 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 6. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum; …; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.59
Sequences extracted from clusters that showed 10 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 2. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.60
Sequences extracted from clusters that showed 10 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 3. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.61
Sequences extracted from clusters that showed 10 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 4. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.62
Sequences extracted from clusters that showed 10 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 5. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.63
Sequences extracted from clusters that showed 10 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 6. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.64
Sequences extracted from clusters that showed 15 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 2. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.65
Sequences extracted from clusters that showed 15 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 3. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.66
Sequences extracted from clusters that showed 15 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 4. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.67
Sequences extracted from clusters that showed 15 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 5. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.68
Sequences extracted from clusters that showed 15 differences or less between the longest sequences in cluster and other sequences in taxonomy at level 6. The taxonomic levels used to asses a discrepancy are as follows: level 1, nuclear or organelle; level 2, kingdom; level 3, phylum;…; level 7, genus; level 8, species. This is based on our own manually derived and homogeneous taxonomy for Eukaryota. There is no space in definition line.69
Sequence alignments that allowed to detect (compute) introns. Gene sequence is followed by intron sequence.


## 4. Concluding Remarks

Most of the introns contained in SSU-rRNA gene sequences have been described as group I introns, only 28 sequences being described as group II. Group I introns are a distinct class of RNA enzymes (ribozymes) characterized by a conserved RNA primary and secondary structures (http://www.rna.icmb.utexas.edu/) essential for splicing and are often capable of self-splicing. The sporadic and wide distribution of group I introns in the nuclear-encoded or organelle-encoded rRNA gene sequences of green and red algae, fungi, ciliates, and different amoebae suggests that these sequences have been highly successful in invading and maintaining in eukaryotic genomes [19, 28].Phylogenetic distribution of group I introns is widespread but often sporadic, which strongly suggests that these sequences are mobile genetic elements capable of some although rare horizontal transfers between evolutionarily distinct lineages or alternatively that they are old elements that can be often deleted.Despite extensive literature and documentation, it was not our purpose to enter a detailed analysis of thenew introns we found in clades in which they had not been yet described. This is why we separated our dataset in three categories.(i)
_GC. These sequences can be used in the future to improve our algorithm, provided that introns described in these entries have been properly described. Without any mention in the exact method used to identify introns that have been described, we were not able to pursue this issue.
(ii)
_UC. These sequences are the most likely sequences containing genuine introns. These sequences could be added to any database of SSU-rRNA sequences when biodiversity analyses are done from amplicons derived from rRNA sequences rather than rRNA gene sequences.
(iii)
_RC. These sequences were unexpected (at least in such amount) but could not be presently explored for the same reasons as the _GC sequences.We present a large series of files that could be helpful for future manual analyses of introns in SSU-rRNA sequences. We hope our analyses will be helpful for several studies.(i)
They will help in the analyses of intronic sequences themselves, in terms of structure, insertion sites, and so forth.
(ii)
They will help in the analyses of intronic sequences in terms of history, that is, deciphering in terms of phylogeny, the time at which a given type of intron invaded the SSU-rRNA gene from another location and was then either kept of lost (using the sets of clusters having a heterogeneous taxonomy).
(iii)
Comparative biodiversity analyses: such analyses can be based either on amplicons derived from gene sequences (DNA) or on amplicons derived from rRNA sequences following a reverse-transcription step. Comparisons of both datasets are, for example, used to compare total diversity (DNA based) and diversity of most active organisms (rRNA based). Searching introns in DNA-based sequences should help such comparisons, because we expect sequences with long introns to be easily retrieved as rRNA but no so as rDNA, which would lead to fallacious estimates of which species are active or not.
(iv)
They will be helpful also in the building of well-annotated reference sequence databases of rDNA and rRNA sequences to be used for biodiversity analyses of NGS datasets for which a manual approach is clearly not feasible.The next step envisioned in our analysis will be to find out PCR primers which are most often used to obtain amplicons of SSU-rRNA gene sequences, build a database of such amplicons, and analyze domains often used for diversity analyses and for which long introns are present in some clades. Indeed with the new NGS approaches of biodiversity, a filter on lengths is de facto applied and absence of a given clade in a sample could only result from long introns or multiple introns in the gene-amplified domain.



---

*Source: 854869-2013-02-06.xml*