Secondary loss of a cis- spliced intron during the divergence of Giardia intestinalis assemblages
© Kamikawa et al.; licensee BioMed Central Ltd. 2014
Received: 19 November 2013
Accepted: 20 June 2014
Published: 30 June 2014
Giardia intestinalis is a parasitic unicellular eukaryote with a highly reduced genome, in which only six cis- spliced and four trans- spliced introns have been discovered. However, we anticipate that more cis- and trans- spliced introns likely remain unidentified in genes encoding hypothetical proteins that occupy ca. 2/3 of all of the open reading frames (ORFs) in the Giardia genome. Consequently, comprehensive surveys of introns in ORFs for hypothetical proteins are critical for better understanding of the intron evolution in this organism.
In this study, we identified two novel cis-spliced introns in the draft genome data of G. intestinalis strain WB, by surveying the conserved sequence motifs shared amongst the previously known introns. G. intestinalis strains can be divided into phylogenetically distinct assemblages A–H, and all the introns identified in past studies are shared among the published genome data from strains WB, DH, GS, and P15 representing assemblages A1, A2, B, and E, respectively. Nevertheless one of the two novel introns identified in this study was found to be absent in strain P15.
By considering the organismal relationship among G. intestinalis assemblages A1, A2, B, and E, one of the two introns identified in this study has highly likely been lost after the divergence of the assemblages. On the basis of a sequence comparison between the intron-bearing loci in WB, DH, and GS genomes and the homologous but intron-free locus in P15 genome, we propose that the loss of this particular intron was mediated by integration of the DNA fragment reverse-transcribed from mature mRNAs.
Spliceosomal introns, which are excised from pre-mature mRNAs by RNA-protein complexes called spliceosomes , are one of the features exclusively found in eukaryotic genomes. However, a large variety in intron density has been found across eukaryotic genomes sequenced to date . In the human genome, for example, 8.4 introns on average are annotated per gene , and the mean intron size is ca. 3,000 bp in length . In contrast, Giardia intestinalis, a unicellular eukaryotic parasite belonging to the Diplomonadida (Excavata) is known to possess a highly reduced genome of only 12 Mbp in length . One of the prominent natures of the Giardia genome is its low intron density—only 6 cis-spliced introns and 4 trans-spliced introns (split introns) have been identified prior to this study [4–12]. Henceforth here, we simply designate cis-spliced introns as ‘introns,’ and trans-spliced introns as ‘splintrons’ .
Most of introns/splintrons in the Giardia genome were identified principally as non-coding stretches intervening in open reading frames (ORFs) encoding proteins shared amongst phylogenetically diverse eukaryotes. However, the simple procedure described above may be problematic for distinguishing the coding and non-coding regions (i.e. exons and introns/splintrons) in functionally unidentified ORFs encoding Giardia-specific proteins. Since unidentified ORFs occupy approximately 2/3 of the ca. 9,000 ORFs encoded in the Giardia genome , a large fraction of introns/splintrons in the Giardia genome may have been overlooked by pioneering surveys principally based on sequence similarity.
To shed light on introns/splintrons veiled in unidentified ORFs in the genome of G. intestinalis strain WB, we conducted an intron survey based on the conserved sequence motifs in introns/splintrons, and successfully detected two novel introns in unidentified ORFs (Note that our approach is not technically applicable to survey splintrons). The two ORFs, which harbor introns in the WB genome, were identified in the genomes of G. intestinalis strains DH , GS  and P15  as well, but one of these in the P15 genome were found to be intron-free. We propose a scenario to explain the presence/absence of the particular intron in the four G. intestinalis strains.
In silico detection of conserved intron sequences
Giardia introns/splintrons known to date bear conserved sequence motifs at the 5′ and 3′ termini, 5′-STATG-3′ and 5′-HCTRACMCVCAG-3′ (R = A or G; H = A, T, or C; M = A or C; V = A, C, or G; S = G or C), respectively. Furthermore, the two motifs may be flanked with each other within 300 bp, since all of the known introns in the Giardia genome range from 29 to 220 bp in length. We searched for genome segments that satisfied the above criteria in the draft genome data of G. intestinalis strain WB (GiardiaDB, http://www.giardiadb.org/giardiadb/).
Cells, DNA, RNA, and reverse transcription
G. intestinalis strain WB (ATCC50803) was cultivated as described previously . Genomic DNA (gDNA) was extracted by cetyl trimethylammonium bromide buffer  from the harvested cells. Total RNA was isolated from the cells with the RNeasy Plant Mini kit (QIAGEN) following the manufacturer’s instruction. To synthesize cDNA from total RNA, reverse transcription was performed by the 3′ rapid amplification of cDNA ends kit (Invitrogen) following the manufacturer’s instruction.
Detection of intron splicing
We designed exact-match primers at the 5′ and 3′ flanking regions of intron-like sequences nominated by the in silico survey (see above), and performed two separate PCRs, one with total gDNA as the template (gDNA-based PCR) and the other with cDNA as the template (cDNA-based PCR). If the particular candidate is an intron, the amplicons from cDNA-based PCR should be shorter than those from gDNA-based PCR by the intron length. We examined all intron-like sequences by comparing size difference between gDNA-based and cDNA-based PCR amplicons. In addition, we sequenced the cDNA-based PCR amplicons to assess whether intron-like sequences were excised. The experimental procedures above identified that two intron-like sequences, one in ORF no. AACB02000068-1-10039-10248 and the other in ORF no. AACB02000001-6-305427-304747, were excised in vivo (see below). Sets of primers 5′-GAAAAAAAATCCAGAGATGGC-3′ and 5′-TTGCAAAGTGCAATGAAAGC-3′, and 5′-AAACAGGTTCGTCAATATCAC-3′ and 5′-AGGATACGAAGCGTTGCGAA-3′ were used for examining the former and latter introns, respectively. Amplified PCR products were cloned into the pGEM T-easy vector (Promega) and sequenced completely.
We experimentally determined the 5′ ends of the mRNAs from ORFs AACB02000068-1-10039-10248 and AACB02000001-6-305427-304747 by using the 5′ rapid amplification of cDNA ends kit (Invitrogen) following the manufacturer’s instruction. Cloning and sequencing were performed as described above.
Results and discussion
The two introns were found in ORFs no. AACB02000068-1-10039-10248 and AACB02000001-6-305427-304747, shown in Figure 1C and D, respectively. Hereafter, we designate the ORFs AACB02000068-1-10039-10248 and AACB02000001-6-305427-304747 as orfA and orfB, respectively. Each of the novel introns locates at the 5′ terminal region of the corresponding ORF, as seen in the previously identified Giardia introns, except for that found at the 3′ terminal region of rpl7A gene. In terms of intron length, both the intron in orfA (41 bp) and that in orfB (33 bp) are comparable with other Giardia introns (29–36 bp), except for the intron in the ORF encoding hypothetical protein GL50803_35332 (220 bp) and that in rpl7A gene (109 bp). Based on the experimentally confirmed 5′ ends of the mRNAs, the gene models for the two ORFs, as well as their intron-exon boundaries, were refined (Figure 1C and D).
There are three major models to explain how eukaryotic genomes lost spliceosomal introns: (i) ‘de-intronization’ by mutations, (ii) non-homologous end joining (NHEJ) repair of double strand break (DSB) in an intron sequence, and (iii) homologous recombination of the cDNA e.g., [21–23]; see also Additional file 2: Figure S2A-C]. The first model assumes the conversion of an intron sequence to an exon sequence by nucleotide substitutions, which results in extension of the corresponding ORF (Additional file 2: Figure S2A). Nevertheless, the length of orfB was found to be uniform among WB, DH, GS, and P15 genomes (Figure 2; see also Additional file 1: Figure S1B), suggesting that the loss of orfB intron cannot be rationalized by deintronization. The second model demands ‘microhomology’ pairing between 5′ and 3′ splice sites to anchor the upstream and downstream exons, which are split by DSB in the intron, during NHEJ repair ; see also Additional file 2: Figure S2B]. Importantly, both 5′ and 3′ splice sites need to be 5′-AG/GT-3′ (the slash indicates the boundary between intron and exon) in the second intron loss model (Additional file 2: Figure S2B). As the key assumption does not fit to the splice sites of the orfB introns (Figure 2), suggesting that the intron has not been eliminated from the P15 genome by NHEJ repair. The last model invokes integration of a reverse-transcribed mRNA (i.e. cDNA; intron-free) into the original, intron-containing locus through homologous recombination (Additional file 2: Figure S2C). We regard that the homologous recombination of the cDNA, which can eliminate the entire intron sequence but does not require sequence conservation at the 5′ and 3′ splice sites, is more appropriate to explain the loss of orfB intron in the P15 genome than the two models described above, (Figure 2). It is intriguing to point out the presence of putative reverse transcriptase genes in the genome data of the four G. intestinalis assemblages [e.g., Genbank/EMBL/DDBJ accession nos. AF434198 (WB), AHGT01000152 (DH), EES99684 (GS), and EFO60876 (P15)], although reverse transcription activity of these encoded proteins has yet to be experimentally confirmed in G. intestinalis cells .
In this study, we found two novel cis- spliced introns and their punctate distribution in the genomes of G. intestinalis assemblies. Together with the recently found trans-spliced introns, the data presented here suggest that the intron evolution in this organism is more complex than we previously thought.
We thank Dr. Aaron A. Heiss (Univ. Tsukuba) for critical comments and English-language corrections. This work was supported in part by grants from JSPS awarded to RK (no. 24870004), YI (no. 21370031, 23117006) and TH (no. 23117005, 23405013, 23247038).
- Wahl MC, Will CL, Lührmann R: The spliceosome: design principles of a dynamic RNP machine. Cell. 2009, 136: 701-718.PubMedView ArticleGoogle Scholar
- Keeling PJ, Slamovits CH: Causes and effects of nuclear genome reduction. Curr Opn Genet Dev. 2005, 15: 601-608.View ArticleGoogle Scholar
- Rodriguez-Trelles F, Tarrio R, Ayala FJ: Origins and evolution of spliceosomal introns. Annu Rev Genet. 2006, 40: 47-76.PubMedView ArticleGoogle Scholar
- Morrison HG, McArthur AG, Gillin FD, Aley SB, Adam RD, Olsen GJ, Best AA, Cande WZ, Chen F, Cipriano MJ, Davids BJ, Dawson SC, Elmendorf HG, Hehl AB, Holder ME, Huse SM, Kim UU, Lasek-Nesselquist E, Manning G, Nigam A, Nixon JE, Palm D, Passamaneck NE, Prabhu A, Reich CI, Reiner DS, Samuelson J, Svard SG, Sogin ML: Genomic minimalism in the early diverging intestinal parasite Giardia lamblia. Science. 2007, 317: 1921-1926.PubMedView ArticleGoogle Scholar
- Nixon JE, Wang A, Morrison HG, McArthur AG, Sogin ML, Loftus BJ, Samuelson J: A spliceosomal intron in Giardia lamblia. Proc Natl Acad Sci U S A. 2002, 99: 3701-3705.PubMedPubMed CentralView ArticleGoogle Scholar
- Russell AG, Shutt TE, Watkins RF, Gray MW: An ancient spliceosomal intron in the ribosomal protein L7a gene (Rpl7a) of Giardia lamblia. BMC Evol Biol. 2005, 5: 45-PubMedPubMed CentralView ArticleGoogle Scholar
- Kamikawa R, Inagaki Y, Tokoro M, Roger AJ, Hashimoto T: Split introns in the genome of Giardia intestinalis are excised by spliceosome-mediated trans-splicing. Curr Biol. 2011, 21: 311-315.PubMedView ArticleGoogle Scholar
- Kamikawa R, Inagaki Y, Roger AJ, Hashimoto T: Splintrons in Giardia intestinalis: spliceosomal introns in a split form. Comm Integr Biol. 2011, 4: 454-456.View ArticleGoogle Scholar
- Nageshan RK, Roy N, Hehl AB, Tatu U: Post-transcriptional repair of a split heat shock protein 90 gene by mRNA trans-splicing. J Biol Chem. 2011, 286: 7116-7122.PubMedPubMed CentralView ArticleGoogle Scholar
- Kamikawa R, Inagaki Y, Hashimoto T: A novel spliceosome-mediated trans-splicing can change our view on genome complexity of the divergent eukaryote Giardia intestinalis. Biophys Rev. 2012, 3: 193-197.View ArticleGoogle Scholar
- Roy SW, Hudson AJ, Joseph J, Yee J, Russell AG: Numerous fragmented spliceosomal introns, AT-AC splicing, and an unusual dynein gene expression pathway in Giardia lamblia. Mol Biol Evol. 2012, 29: 43-49.PubMedView ArticleGoogle Scholar
- Franzén O, Jerlström-Hultqvist J, Einarsson E, Ankarklev J, Ferella M, Andersson B, Svärd SG: Transcriptome profiling of Giardia intestinalis using strand-specific RNA-Seq. PLoS Comput Biol. 2013, 9: e1003000-PubMedPubMed CentralView ArticleGoogle Scholar
- Adam RD, Dahlstrom EW, Martens CA, Bruno DP, Barbian KD, Ricklefs SM, Hernandez MM, Narla NP, Patel RB, Porcella SF, Nash TE: Genome sequencing of Giardia lamblia genotypes A2 and B isolates (DH and GS) and comparative analysis with the genomes of genotypes A1 and E (WB and Pig). Genome Biol Evol. 2013, 5: 2498-2511.PubMedPubMed CentralView ArticleGoogle Scholar
- Franzén O, Jerlström-Hultqvist J, Castro E, Sherwood E, Ankarklev J, Reiner DS, Palm D, Andersson JO, Andersson B, Svärd SG: Draft genome sequencing of Giardia intestinalis assemblage B isolate GS: Is human giardiasis caused by two different species?. PLoS Pathog. 2009, 5: e1000560-PubMedPubMed CentralView ArticleGoogle Scholar
- Jerlström-Hultqvist J, Franzén O, Ankarklev J, Xu F, Nohynkova E, Andersson JO, Svard SG, Andersson B: Genome analysis and comparative genomics of a Giardia intestinalis assemblage E isolate. BMC Genomics. 2010, 11: 543-PubMedPubMed CentralView ArticleGoogle Scholar
- Kamikawa R, Nishimura H, Sako Y: Analysis of the mitochondrial genome, transcripts, and electron transport activity in the dinoflagellate Alexandrim catenella (Gonyaulacales, Dinophyceae). Phycol Res. 2009, 57: 1-11.View ArticleGoogle Scholar
- Lasek-Nesselquista E, Welch DM, Sogin ML: The identification of a new Giardia duodenalis assemblage in marine vertebrates and a preliminary analysis of G. duodenalis population biology in marine systems. Int J Parasitol. 2010, 40: 1063-1074.View ArticleGoogle Scholar
- Nash TE, Keister DB: Differences in excretory-secretory products and surface antigens among 19 isolates of Giardia. J Infect Dis. 1985, 152: 1166-1171.PubMedView ArticleGoogle Scholar
- Monis PT, Andrews RH, Mayrhofer G, Ey PL: Molecular systematics of the parasitic protozoan Giardia intestinalis. Mol Biol Evol. 1999, 16: 1135-1144.PubMedView ArticleGoogle Scholar
- Monis PT, Caccio SM, Thompson RCA: Variation in Giardia: towards a taxonomic revision of the genus. Trends Parasiol. 2009, 25: 93-100.View ArticleGoogle Scholar
- Roy SW, Gilbert W: The evolution of spliceosomal introns: patterns, puzzles and progress. Nat Rev Genet. 2006, 7: 211-221.PubMedGoogle Scholar
- Farlow A, Meduri E, Schlotterer C: DNA double-strand break repair and the evolution of intron density. Trends Genet. 2011, 27: 1-6.PubMedPubMed CentralView ArticleGoogle Scholar
- Cohen NE, Shen R, Carmel L: The role of reverse transcriptase in intron gain and loss mechanisms. Mol Biol Evol. 2012, 29: 179-186.PubMedView ArticleGoogle Scholar
- Burke WD, Malik HS, Rich SM, Eickbush TH: Ancient lineages of non-LTR retrotransposons in the primitive eukaryote, Giardia lamblia. Mol Biol Evol. 2002, 19: 619-630.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.