Footprints of primordial introns on the eukaryotic genome

Another expected legacy of ancient phase-zero introns is that putatively ancient introns (those widely phylogenetically distributed) should have a stronger phase-zero bias. We tested this notion by examining widely distributed C. elegans intron positions: those present in two or more kingdoms. We used a database of ancient genes with known corresponding protein structures (culled from the Protein Data Bank) and the introns from their homologues, compiled by Fedorov and collaborators (unpublished). Of 6568 total intron positions, C. elegans introns are found at 1539 positions within 280 gene families. At 233 of these C. elegans positions there are also introns in the copies of the gene from non-animal organisms. Such coincident positions are candidates for the positions of ancient introns, based on this very wide phylogenetic distribution. (These coincident intron positions do, in fact, show a stronger correlation with module boundaries 6xCentripetal modules and ancient introns. Roy, S.W. et al. Gene. 1999; 238: 85–91Crossref | PubMed | Scopus (40)See all References6.) If these positions are indeed ancient, and if the phase-zero bias of introns is due to the presence of ancient introns, this ancient subset of C. elegans introns should have a higher phase-zero bias than should the set of animal-unique C. elegans introns.Table 2Table 2 shows that this ancient subset does have a significantly stronger phase-zero bias (62.2% versus 50.8%), with a χ2 of 10.4 (P = 0.0013). This is a striking verification of the expectation that part of the phase-zero bias is due to the presence of truly ancient introns. Table 3Table 3 shows further that this is a general phenomenon. There is a highly-significant difference between the phase bias of animal intron positions at which introns are found in two or more kingdoms and those positions that are exclusively animal. The same is true for invertebrates and vertebrates separately, for plants, and for fungi, and there is a significant difference for protists. This trend is seen even further in positions where introns are found in three or more kingdoms. This result suggests that there were in fact ancient intron positions and that the observed excess of phase-zero positions is, at least in part, due to their presence in modern eukaryotes.Table 2Phase patterns of introns in Caenorhabditis elegans genesPhase zeroχ2Phase onePhase twoTotal intronsExclusively animal C. elegans positions663 50.8%306 23.4%337 25.8%1306Widely distributed C. elegans positions145 62.2%10.446 19.7%42 18.0%233Widely distributed signifies that an intron is also found at that position in a non-animal copy of the gene. The introns are from the database of 6568 intron positions in the set of ancient genes that have three-dimensional structures for the corresponding proteins compiled by Fedorov and collaborators (unpublished). The χ2 is calculated comparing phase zero with (phase one + phase two).Table 3Phase patterns of introns in different kingdomsaPhase zeroχ2Phase onePhase twoTotal intronsAnimalExclusively animal1607 50.2%860 26.9%735 23.0%3202Animal and other kingdom(s)329 59.8%17.4127 23.1%94 17.1%550Exclusively vertebrate504 48.2%324 31.0%217 20.8%1045Vertebrate and other kingdom(s)194 61.8%17.869 22.0%51 16.2%314Exclusively invertebrate906 51.6%423 24.1%428 24.4%1757Invertebrate and other kingdom(s)246 58.9%7.296 23.0%76 18.2%418 PlantExclusively plant1166 59.9%398 20.5%381 19.6%1945Plant and other kingdom(s)323 66.1%6.199 20.2%67 13.7%489 FungiExclusively fungi303 42.7%236 33.3%170 24.0%709Fungi and other kingdom(s)118 54.1%8.756 25.7%44 20.2%218 ProtistExclusively protist34 38.6%28 31.8%26 29.5%88Protists and other kingdom(s)41 56.9%5.316 22.2%15 20.8%72 AllOnly one kingdom3110 52.3%1522 25.6%1312 22.1%5944Two or more kingdoms371 60.4%14.7142 23.1%101 16.4%614Three or more kingdoms51 64.6%13 16.5%15 19.0%79The introns are from the database of 6568 intron positions in the set of ancient genes that have three-dimensional structures for the corresponding proteins compiled by Fedorov and collaborators (unpublished). The χ2 is calculated comparing phase zero with (phase one + phase two).These twin findings – first, that the intron phases in ancient genes, defined by BLAST scores, differ from those of recently transferred genes by having a stronger phase-zero bias, and second, that phylogenetically ancient intron positions also show a stronger phase-zero bias – support the theory that introns, mostly in phase zero, were present in the LUCA.Furthermore, as phase-zero, widely phylogenetically-distributed intron positions are highly correlated with the boundaries of modules 6xCentripetal modules and ancient introns. Roy, S.W. et al. Gene. 1999; 238: 85–91Crossref | PubMed | Scopus (40)See all References6, the subset of introns possessing one characteristic expected of ancient intron positions – wide phylogenetic distribution – also possesses two other traits expected of ancient positions. This is a result easily explained on a mixed intron-origin model, but not easily explained on strictly insertional ones.