The dynamic intein landscape of eukaryotes
Mobile DNA volume 9, Article number: 4 (2018)
Inteins are mobile, self-splicing sequences that interrupt proteins and occur across all three domains of life. Scrutiny of the intein landscape in prokaryotes led to the hypothesis that some inteins are functionally important. Our focus shifts to eukaryotic inteins to assess their diversity, distribution, and dissemination, with the aim to comprehensively evaluate the eukaryotic intein landscape, understand intein maintenance, and dissect evolutionary relationships.
This bioinformatics study reveals that eukaryotic inteins are scarce, but present in nuclear genomes of fungi, chloroplast genomes of algae, and within some eukaryotic viruses. There is a preponderance of inteins in several fungal pathogens of humans and plants. Inteins are pervasive in certain proteins, including the nuclear RNA splicing factor, Prp8, and the chloroplast DNA helicase, DnaB. We find that eukaryotic inteins frequently localize to unstructured loops of the host protein, often at highly conserved sites. More broadly, a sequence similarity network analysis of all eukaryotic inteins uncovered several routes of intein mobility. Some eukaryotic inteins appear to have been acquired through horizontal transfer with dsDNA viruses, yet other inteins are spread through intragenomic transfer. Remarkably, endosymbiosis can explain patterns of DnaB intein inheritance across several algal phyla, a novel mechanism for intein acquisition and distribution.
Overall, an intriguing picture emerges for how the eukaryotic intein landscape arose, with many evolutionary forces having contributed to its current state. Our collective results provide a framework for exploring inteins as novel regulatory elements and innovative drug targets.
Intervening protein elements, called inteins, are capable of self-excision from precursor polypeptides and ligation of the flanking sequences, termed exteins, to produce the functional host protein [1, 2]. As mobile elements, some inteins also contain a homing endonuclease domain (HEN), allowing movement at the DNA level [3, 4]. Inteins remove themselves from host polypeptides through a process known as protein splicing (Fig. 1a, Intein). This excision occurs in four steps using conserved splice junction residues, often highly reactive nucleophiles, like cysteine, serine, or threonine (Additional file 1: Figure S1).
In addition to inteins, eukaryotes encode other self-splicing proteins that utilize similar chemistry. Hedgehog proteins are structural and functional analogs of inteins, but are involved in cell signaling in higher eukaryotes exclusively  (Fig. 1a, Hedgehog; Additional file 1: Figure S2). These essential proteins are synthesized as inactive precursors and undergo nucleophilic attack similar to protein splicing in order to form a functional molecule with cholesterol . Eukaryotes also encode proteins with Hedgehog-intein (Hint) domains (Fig. 1a, Hint; Additional file 1: Figure S2). Hints employ similar autoproteolytic reactions, either mirroring activity at the N-terminus or the C-terminus of an intein, depending on where peptide cleavage occurs . All three of these eukaryotic self-splicing elements share a canonical Hint fold to promote peptide cleavage reactions .
Inteins were originally discovered while studying the vacuolar ATPase (VMA1) of Saccharomyces cerevisiae [9, 10]. A central region of the ATPase that had no homology to other known proton pumps was observed. This ‘spacer’ portion, present in the mRNA but not the final protein, turned out to be a 50-kDa self-splicing intein. Since this discovery, thousands of inteins have been documented mainly through sequence-based approaches. A striking outcome has been uncovering inteins across all three domains of life.
In addition to the intein in VMA1, inteins were discovered across other fungal phyla, interrupting diverse genes, such as the pre-mRNA processing factor 8 (Prp8) [11,12,13,14]. Experimental data have shown that these fungal inteins are splicing active, at least in foreign exteins, and some were shown to be mobile [15,16,17,18]. Additionally, the VMA1 intein was crystallized, providing the first insight into an intein HEN . Chloroplast genes from some algae also contain some of these self-splicing elements [20, 21], and various inteins in eukaryotic viruses were reported [22,23,24,25]. The existence of inteins in both prokaryotic and eukaryotic viruses lead to the hypothesis that some inteins may disseminate through this intermediary route [23, 24, 26, 27]. Nevertheless, inteins are found overwhelmingly in the genomes of bacteria and archaea, and a comprehensive picture of the eukaryotic intein landscape is currently lacking [2, 4, 28,29,30,31].
A large-scale survey was recently performed to analyze inteins across bacteria and archaea . The results revealed that about half of sequenced archaeal genomes contain at least one intein, whereas a quarter of bacteria are intein-positive among those sequenced. Within this broad distribution, over 60% of proteins containing inteins are involved in replication, recombination, and repair, with 70% localizing to ATP-binding proteins. A similar study investigated where inteins are distributed across bacteria and their phages, providing the first evidence that mycobacteriophages function as facilitators of intein dissemination across mycobacteria .
With the recent surge in available sequenced genomes, we next focused on inteins in eukaryotic genes to learn how these inteins might have evolved and spread. Eukaryotic inteins are not only present in nuclear genomes, but in genomes of the chloroplast, and in viral genomes. These inteins frequently interrupt proteins involved in replication, recombination, and repair, as well as RNA processing and energy metabolism, differing somewhat from the tendencies of bacterial and archaeal inteins. Curiously, inteins seem to be enriched in specific lineages of pathogenic fungi, such as Cryptococcus spp. More generally, the patterns of eukaryotic inteins suggest important roles in the regulation of the host protein to the potential benefit of the organism in which it resides. Eukaryotic inteins can be classified into three types based on size and presence of a HEN. Apparently, some of these HENs were or are active, based on connections between diverse inteins on a sequence similarity network that is indicative of movement. We discovered that a specific chloroplast intein is a relic of endosymbiosis several billion years ago, while other chloroplast inteins appear to be horizontally transferred by viral vectors. This work suggests a complex, evolving picture of inteins across the eukaryotic mobilome.
Eukaryotic inteins are scarce
To assess and characterize the diversity of eukaryotic inteins, we surveyed genomic sequences using BLASTp and previously developed pipelines [31, 32]. At the time of analysis (February 20, 2017), there were 11,001 eukaryotic genomes available through the National Center for Biotechnology Information (NCBI), including 2377 entries for nuclear genomic sequences and 8624 entries for organelle- and plasmid-only sequences. We also analyzed ~ 1500 viral reference proteomes available at ViralZone . The full list of intein-containing genomes and other relevant information is provided in Additional file 2: Table S1 and Additional file 3: Table S2. As an additional point of reference for understanding intein evolution in eukaryotes, we mined genomic sequences for other self-splicing protein elements, including Hedgehog, Hedgehog-like proteins, and Hedgehog-intein (Hint) domains (Fig. 1a; Additional file 1: Figure S2 and Table S3). In total, 257 inteins were identified in 231 eukaryotic species, either in their nuclear (nDNA) or chloroplast (cpDNA) genomes, and 136 inteins were found in proteomes from 98 viruses (vDNA) (Fig. 1b).
In general, eukaryotic inteins are scarce compared to bacterial and archaeal inteins. We find that they are widely and sporadically distributed across the eukaryotic tree with the highest number of intein-positive species observed in Fungi, mostly in Ascomycota (141 inteins, or 55% of total; Fig. 1b). Ascomycetes with inteins represent some notable pathogenic Fungi, such as Candida spp.  and Aspergillus spp. . There are 15 inteins (~ 6%) found in Basidiomycota, the close relatives of ascomycetes. Among others, intein-containing basidiomycetes include human pathogens, such as Cryptococcus neoformans and C. gattii [11, 12], and plant pathogens Tilletia indica and T. walkeri . A few intein-containing species are identified among other Fungi (Fig. 1b; Table 1 and Additional file 2: Table S1). A noteworthy intein-containing chytrid fungus is Batrachochytrium dendrobatidis [12,13,14], the cause of chytridiomycosis that is devastating amphibian populations .
A large set of inteins (71 or 28%) was identified in cpDNA from diverse algae and seaweeds (Fig. 1b). Inteins are present in 45 red algae (Rhodophyta), 14 green algae (Chlorophyta), three cryptophytes (Cryptophyta), and nine brown algae and seaweeds (Heterokonta) (Fig. 1b; Table 1; Additional file 2: Tables S1 and Additional file 1: Table S3). Interestingly, no inteins are found in mitochondrial genomes. One green algal species, Auxenochlorella protothecoides, has an intein in nDNA. Additionally, inteins are found in the nuclear genomes of some social amoebae (Amoebozoa), as well as in one protozoan species (Apusozoa) (Fig. 1b).
Remarkably, despite intein rarity among eukaryotes, there are species with more than one intein per genome, which might indicate intragenomic spread (discussed below). For example, among Fungi, two inteins are identified in nDNA of an arbuscular mycorrhizal fungus Rhizophagus irregularis (Mucoromycota), and a phytopathogenic fungus Fusarium fujikuroi (Ascomycota) carries three inteins in its genome.
Among known viruses of eukaryotes, there are 117 inteins across four families: Iridoviridae (14 inteins), Marseilleviridae (20 inteins), Phycodnaviridae (22 inteins), and Mimiviridae (61 inteins). These viral families belong to the nucleocytoplasmic large DNA viruses (NCLDV), also known as order Megavirales, representatives of which are characterized by extremely large genome sizes . An additional 19 inteins are identified in genomes of unclassified viruses (Fig. 1b).
Intein enrichment in pathogens is genus-specific
While compiling data on intein distribution, we noticed that inteins are seemingly prevalent in fungal pathogens (Additional file 1: Table S4). To investigate further, we analyzed a distribution of pathogen and non-pathogen representatives among sequenced Ascomycota and Basidiomycota genomes to reveal a potential bias introduced by sequencing pathogenic species (Fig. 2a; Additional file 1: Table S5). It is estimated that roughly 33% of all fungi are pathogens, and a majority of them are ascomycetes, with some basidiomycetes . Out of 689 sequenced genomes from Ascomytoca and Basidiomycota, only 30.2% (208 genomes) are from pathogenic fungi, suggesting that there is no bias in the dataset (Fig. 2a, left, gray versus black circles; Additional file 1: Table S5). There also appears to be no bias towards sequencing genomes from pathogenic species in the separated Ascomycota and Basidiomycota subsets (Fig. 2a, right, gray versus black circles; Additional file 1: Table S5). The overall percentage of intein-containing pathogens from a combined analysis of Ascomycota and Basidiomycota (37.0%) is higher than the percentage of intein-containing non-pathogens (16.4%) (Fig. 2a, left, red circles). The trend remains when looking at the two phyla separately. Ascomycota have 41.8% pathogenic species with inteins versus 27.4% non-pathogenic species with inteins (Fig. 2a, right, red circle). Dramatically, Basidiomycota alone have 6-fold more intein-containing pathogens than non-pathogens, with 18.6% compared to 3.2% (Fig. 2a, right, red circle).
To further investigate this phenomenon, we focused on two intein-containing fungal groups: Aspergillus/Neosartorya (Ascomycota) and tremellomycetous yeasts (Basidiomycota), including C. neoformans, C. gattii, and their close relatives (Fig. 2b) [39,40,41,42]. We found that within a subset of selected Aspergillus, there was no correlation between having an intein and a pathogenic lifestyle, despite several important pathogens having inteins (Fig. 2b, left; Additional file 1: Table S6). On the other hand, within the tremellomycetous yeasts, there was a positive correlation between pathogenic species and intein-containing species. Indeed, C. neoformans and C. gattii are the only two pathogenic species within the tremellomycetous yeast group, and are the only two intein-containing species in the analysis (Fig. 2b, right; Additional file 1: Table S7). It seems that pathogenic fungi do have a propensity for inteins, but the pattern is variable among specific genera.
Inteins are found in genes specific to genome type (nDNA, cpDNA or vDNA)
Next, we analyzed the distribution of inteins relative to the protein into which they are inserted (extein) (Fig. 3a–3c; Table 1). Interestingly, there is little or no overlap in exteins between different types of genomes that harbor inteins (Fig. 3a–3c; Table 1). While in nDNA, a majority of inteins are found in Prp8 and VMA1 (Fig. 3a; Additional file 1: Table S8), the replicative helicase DnaB and DNA-dependent RNA polymerase (DdRP) are the most common intein-containing proteins in cpDNA (Fig. 3b; Additional file 1: Table S9). Only red and brown algae contain DnaB inteins, as green algae lack a dnaB gene, and the reverse is observed for the ATP-dependent Clp protease, ClpP, and its inteins (Additional file 1: Figure S3). In viral genomes, DNA-dependent DNA polymerases (DdDPs) often harbor inteins (Fig. 3c; Additional file 1: Table S10). The only overlapping intein-containing proteins among all three groups of genomes are DdRPs. However, it should be noted that DdRPs are diverse and genome-specific. Also, the Prp8 intein exists in both nDNA and cpDNA genomes, although its presence in algae might be the result of fungal contamination.
We also sorted exteins into functional caterogies using the KOG (EuKaryotic Orthologous Groups) database for nDNA and vDNA, or the COG (Clusters of Orthologous Groups) database for cpDNA [43,44,45] (Fig. 3d; Table 1; Additional file 1: Tables S8-S11). Analysis of functional categories indicates a bias toward category A (RNA processing and modification) for intein-containing proteins from nDNA (108 proteins with intein insertion), and toward category L (replication, recombination, and repair) for exteins from cpDNA (56 proteins) and vDNA (69 proteins) (Fig. 3d; Additional file 1: Table S11). As expected, category K, represented by proteins involved with transcription, is shared by nDNA (21 proteins), cpDNA (12 proteins), and vDNA (10 proteins) for insertions in DdRP exteins. Additionally, category C (energy production and conversion) is prominent among nDNA exteins (42 proteins), and category F (nucleotide metabolism and transport) is overrepresented among vDNA exteins (35 proteins). The other functional categories are less represented (Fig. 3d; Table 1).
Inteins are located in conserved structure boundaries
Inteins are known to occupy sites with specific characteristics that allow splicing besides preceding a nucleophile . Ideally, intein insertion sites must allow proper folding and splicing, and permit the host protein to correctly fold post-splicing . Additionally, inteins are often found in protein active sites and in other conserved, functionally important domains within host proteins [31, 48, 49]. Thus, we examined native insertion sites of the eukaryotic inteins for the properties such as conservation and secondary structure (Fig. 4; Additional file 1: Figure S4). We note that it is a common practice to distinguish inteins based on their extein identity and insertion site . For example, inteins in the Prp8 protein are identified as Prp8i. There are six intein insertion points found in the Prp8 protein, which were designated as a-f based on the conventional intein insertion site classification scheme [50, 51]. Thus, Prp8i from insertion point a is indicated as Prp8i-a (Fig. 4).
To determine the local secondary structure of the intein insertion site, we used homology modeling (Additional file 1: Table S13). As demonstrated previously, intein insertion sites are more likely to occur in loop-structure boundary (54% of all sites) than in the middle of a β-sheet or α-helix (31% of sites) . Results of our secondary structure modeling are in agreement with these observations. As demonstrated in Fig. 4 above the 3D structures, the most common eukaryotic inteins Prp8i-a and VMA1i-a, as well as less widely distributed inteins in threonyl-tRNA synthetase (ThrRSi-a) and glutamate synthase (GLTi-a), are inserted either in flexible loops, or in close proximity to the loop/α-helix boundary (Fig. 4). Intein insertion in a flexible loop likely allows the intein to fold and splice properly, creating less strain on the host protein.
In order to assess insertion site conservation, the PDB structures of these exteins were uploaded to ConSurf , a bioinformatics tool for estimating the evolutionary conservation of amino acid positions in a protein based on the phylogenetic relations between homologous sequences. The ConSurf alignments of nDNA exteins Prp8, VMA1, GLT, and ThrRS reveal they are all in highly conserved sites, as previously noted by several groups [14, 34, 49] (Fig. 4, yellow highlighting; Additional file 1: Figure S4). The structure models also demonstrate the inteins insert at flexible boundaries, often between structured regions. VMA1i-a is a clear example, being inserted at a loop between a β-turn and an α-helix.
The presence of homing endonucleases (HEN) (see Fig. 5) likely facilitates mobility to conserved regions. Inteins likely home to such positions to decrease the chance of elimination and delay removal, as deletion would require absolute precision to maintain host protein function. This conservation assessment is in line with what has been reported for bacteria and archaea [27, 31, 49]. Two possibilities exist: limitations on mutating the homing site make it difficult for hosts to evolve immunity to the HEN, allowing inteins to propagate as selfish elements, or inteins are functionally important at these conserved sites and are retained to perform yet unknown regulatory roles.
Eukaryotic inteins vary in size
Many inteins carry a site-specific HEN in addition to the protein splicing domain  (Fig. 5a). A HEN renders an intein mobile by introducing a double-strand break into an inteinless allele and initiating a gene conversion event, which results in the copying of the intein coding sequence . This process, known as homing, is thought to be responsible for the horizontal transfer and spread of inteins . Thus, to elucidate the evolutionary dynamics of eukaryotic inteins, it is important to trace the presence/absence of intein-associated HENs and their key features (Fig. 5).
The length of eukaryotic inteins, a simple indicator of the presence of a HEN domain, varies greatly from 123 amino acid residues (DnaBi from red algae Dichotomaria marginata) up to 1065 amino acid residues (Prp8i-b from grass-infecting smut fungus Anthracocystis flocculosa) (Fig. 5a and b). When all identified inteins are plotted by residue length, the resulting plots show a step-wise pattern with varying intein lengths (Fig. 5b; Additional file 1: Figure S5). The lower step, up to ~ 200 amino acid residues, must correspond to inteins that lack HEN domains, referred to as mini-inteins or HEN(−). The middle step between ~200 to 700 amino acids, is likely represented by inteins carrying a HEN or inteins with splicing domains interrupted by a ‘linker’ sequence. These inteins were analyzed for endonuclease motifs and were indeed found to contain LAGLIDADG HENs. The most puzzling aspect of the plot is a group of inteins that are more than 800 amino acid residues in length (Fig. 5b and c). Although there are only a few extra-large inteins and all of them are in Prp8 (Fig. 5c, HEN(+)extra), they exist in diverse fungal species, including both ascomycetes (e.g. F. fujikuroi and A. fumigatus) and basidiomycetes (An. flocculosa). In addition to protein splicing domain and HEN, these extra-large inteins carry linker sequences of unknown origin and function, with some of these linkers having repeated motifs (Fig. 5a). Repetitive sequences are common in proteins for structural reasons or may participate in ligand binding [54,55,56]. The HEN(+)extra inteins were also analyzed for endonuclease type, and many were found to carry LAGLIDADG repetitive motifs (Additional file 2: Table S1 and Additional file 3: Table S2). In general, when inteins in specific exteins are plotted as a function of length, they tend to cluster together in the same HEN category (Fig. 5c).
Sequence similarity network and phylonetwork of eukaryotic inteins
Typically, phylogenetic reconstruction, which focuses on relationships resulting from vertical descent, is used for analysis of homologous genes (or proteins). However, inteins do not follow the rule of strict vertical descent [4, 14, 27, 30, 31], and classic phylogenetic analysis would pose significant limitations to the analysis of intein evolutionary dynamics. Thus, to further elucidate evolutionary dynamics of eukaryal inteins, we built a sequence similarity network (SSN), which does not assume input sequences are homologous (Fig. 6a) [57,58,59,60]. The resulting SSN consists of 410 nodes and 5782 edges (links) representing significant relationships between analyzed sequences (Fig. 6a). A majority of the SSN components consists of the sequences originating from the same genomic pool (nDNA, cpDNA or vDNA), and inteins cluster according to their exteins, with a few notable exceptions.
The largest SSN component (Fig. 6a, clusters 1a-1c) is formed by Prp8i (Fig. 6a, cluster 1a and 1b, red) and a group of viral inteins (Fig. 6a, cluster 1c, blue). A smaller cluster of Prp8i (Fig. 6a, cluster 1b), which is separated from the bulk, forms a bridge between a larger cluster of Prp8i (Fig. 6a, cluster 1a) and a cluster of viral inteins (Fig. 6a, cluster 1c). VMA1i (Fig. 6a, cluster 2, red) cluster into a mass of highly interconnected nodes, indicating a high degree of similarity and, together with yeast HO Hop endonuclease (Fig. 6a, cluster 2, yellow) , form a large SSN subnetwork. Other SSN components worth mentioning include a large group of viral DdDPi (Fig. 6a, cluster 7, blue), a subnetwork of viral ribonucleotide reductase inteins (RIRi) (Fig. 6a, cluster 10, blue), two smaller SSN components composed of DnaBi from cpDNAs of red algae (Rhodophyta, Bangiophyceae) (Fig. 6a, cluster 5a, green) and brown algae (Stramenopiles, Phaeophyceae) (Fig. 6a, cluster 5b, green). A small group of randomly sampled Hedgehog proteins (C-terminal domain only) (Fig. 6a, cluster 11, black) was included into the analysis and seems to form an outgroup.
There are several SSN components that are formed by inteins located in different proteins. For example, cpDNA inteins, ClpPi and DdRPi from green algae (Chlorophyta) fell within a single connected SSN component (Fig. 6a, cluster 6, green), and also, surprisingly, clustered together with viral DdRPi (Fig. 6a, cluster 6, blue), indicating potential intragenomic intein mobility as well as putative horizontal transfer between cpDNA and vDNA.
One of the most puzzling SSN components consists of fungal glutamate synthase inteins (GLTi) and chitin synthase inteins (CHSi) (Fig. 6a, cluster 3, red). Both GLTi and CHSi are found in genomes of F. fujikuroi and Podospora anserina. There are also overlaps in distribution among fungal genomes between Prp8i and CHSi, as well as between VMA1i and GLTi. While inteins within Prp8 and CHS are present in F. fujikuroi and Diaporthe helianthi, the yeast Debaryomyces hansenii carries VMA1i and GLTi. Thus, additional inter- and intraspecies pairwise comparative analysis of these inteins was warranted (Fig. 6b; Additional file 1: Table S12). The amino acid sequence similarity between GLTi and CHSi is unusually high, ranging between 48.3% and 56.9%, with an average of 52.3% (Fig. 6b). In contrast, the amino acid sequence similarities in Prp8i-versus-VMA1i, Prp8i-versus-CHSi, Prp8i-versus-GLTi, VMA1i-versus-CHSi, and VMA1i-versus-GLTi comparisons do not exceed 40%, even within the same genome (Fig. 6b). Unusually high similarity between GLTi and CHSi suggests relatively recent intragenomic intein mobility.
Finally, we were intrigued by the large number of DnaBi in cpDNA from taxonomically diverse and evolutionary distant species of eukaryotes (Fig. 1b, Table 1). In bacteria, DnaB is an intein hot-spot [30, 31]. Therefore, we hypothesized that eukaryal DnaBi was likely inherited from the cyanobacterial progenitor of chloroplasts by endosymbiosis [62,63,64]. Although cpDNA DnaBi show a relatively high degree of protein sequence similarity in pairwise comparisons (data not shown), an overwhelming majority is represented as singletons in the SSN (Fig. 6a, cluster 12, green), which is likely due to difficulties with obtaining meaningful local alignments during the all-against-all BLASTp stage of SSN reconstruction. Thus, to identify possible evolutionary relationships among DnaBi from cpDNA and bacteria, we utilized an alignment-based approach for phylogenetic network inference (Fig. 6c; Additional file 4) . As evident from the resulting phylonetwork (Fig. 6c; Additional file 1: Figure S6), some bacterial and cpDNA DnaBi are distantly related as they cluster together on the same branches, suggesting intein transfer by endosymbiosis.
Inteins can execute an autocatalytic protein splicing reaction and perform post-translational modification of a precursor (Fig. 1a). These features make intein-based tools indispensable in modern protein chemistry and bioengineering [66,67,68]. However, the nature of inteins, their origin, and possible biological roles are still poorly understood, especially in their native context. This work expands our understanding, by providing a comprehensive analysis of intein diversity, distribution, and dissemination in the eukaryotic world. Inteins are found in three genome types: nuclear (nDNA), chloroplast (cpDNA), and viral (vDNA) (Fig. 1b, Fig. 3a-c). Not surprisingly, eukaryotic inteins are inserted in several important, often functionally critical host proteins, matching a global intein trend [14, 27, 31, 49]. The occurrence of inteins in notable human and plant pathogens is a highlight, and will pave the way for further research into intein drug discovery (Fig. 2). Lastly, by applying sequence similarity network (SSN) and phylonetwork analyses, we find evidence for intragenomic mobility, endosymbiotic acquisition, and horizontal transfer of eukaryotic inteins, thus providing an updated picture for how the eukaryotic intein landscape might have evolved (Fig. 7).
Eukaryotic intein landscape and potential domestication
To gain greater insight into inteins, we previously utilized data mining to investigate the distribution of inteins in prokaryotes, and demonstrated that both bacteria and archaea are relatively intein-rich, with 25% and 50% intein-positive genomes, respectively . In contrast, this study found that fewer than 2% of eukaryotes harbor inteins (Fig. 1b, Fig. 3a-c). The paucity of inteins in eukaryotes suggests that their more complex genomes might have expunged these otherwise invasive elements.
Another possibility for the rarity is that inteins became ‘domesticated’ in eukaryotic genomes and lost some of their features. Indeed, Hedgehog-intein (Hint) domains are widely distributed among eukaryal phyla (Fig. 1a and Additional file 1: Figure S2) [6, 7]. The function of a majority of these Hint-containing proteins remains unknown, and the essentiality of protein splicing to the host proteins has yet to be established . An evolutionary link of inteins to Hedgehog proteins is also a tantalizing possibility [5, 8]. Although there is only limited sequence similarity between the two self-splicing sequences , Hedgehog proteins are structurally analogous to inteins, sharing a characteristic β-strand core known as the Hint module . This common architecture allows both Hedgehog and inteins to undergo post-translational cleavage and ligation by similar trans-esterification reactions to form active molecules. These obvious structural and functional similarities argue in favor of their ancestral relatedness.
Intein distribution suggests possible functional roles
It was observed that prokaryotic inteins localize to highly conserved sites of vital proteins [14, 31, 34]. In general, this trend holds true for inteins in eukaryotes and their viruses (Fig. 3a-c; Fig. 4). cpDNA and vDNA inteins parallel patterns of extein preference for prokaryotic inteins, with insertions overwhelmingly in proteins involved in replication, recombination, and repair. nDNA inteins are often located in proteins associated with RNA processing, polymerases, and energy production, again critical proteins (Fig. 3d; Table 1). The presence of inteins in essential genes, coupled with their occurrence at different sites within the same protein, argues that inteins may be important to the host protein. On the other hand, these data may also suggest that inteins are successful molecular parasites, by inserting at essential sites where they are less likely to be deleted.
Indeed, recent experimental work has shed light on inteins modulating the host protein through conditional splicing. Inteins in bacteria and archaea have been shown to respond to salt, redox state, temperature, and ssDNA [69,70,71,72,73,74]. Additionally, evidence from data mining led to the observation that functionally related, but evolutionarily distinct proteins carry inteins, indicating a bias that suggests maintenance due to a regulatory role . There are also often many independent intein insertions into different sites within the same protein or domain. For example, there are 104 varied inteins at different sites in the critical spliceosomal protein Prp8, suggesting several invasion events, and that the inteins are retained because they provide some function. A fascinating hypothesis is that protein splicing acts to regulate RNA splicing, to the selective advantage of the organism.
The finding that inteins are prevalent in human and plant pathogens is also noteworthy (Fig. 2). Given the absence of inteins from metazoan genomes, a pursuit for intein inhibitors as novel antimicrobials is underway [75, 76]. The discovery here of inteins in several more fungal human pathogens strengthens the rationale behind such drug screening efforts. Furthermore, the presence of inteins in agricultural pathogens makes them ripe for exploitation in drug discovery, at a time when fungal diseases are of increasing concern .
Multiple pathways for eukaryotic intein dissemination
To further elucidate eukaryotic intein dynamics, we performed comparative analyses of all known eukaryotic inteins. Phylogenetic analysis of large sets of inteins is often hindered by an inability to produce quality multiple sequence alignments due to a high level of intein sequence diversity. Moreover, there is mounting evidence that horizontal transfer plays a significant role in intein evolution [14, 27, 29,30,31, 58], but classic phylogenetic analyses focus primarily on vertical descent and traditional phylogenetic trees become inefficient. Here we chose to utilize a SSN-based approach (Fig. 6a), which is less constrained and represents simultaneous interrelationships of all sequences based on their pairwise alignments .
Based on the results of our SSN and additional phylonetwork analyses (Fig. 6), we propose a model for intein dissemination in eukaryotes involving vertical inheritance, intragenomic mobility, endosymbiotic and horizontal (extragenomic) acquisitions (Fig. 7). The wide distribution of nDNA inteins and Hint domains among eukaryotes suggests vertical descent from an ancestral sequence as far back as the last universal common ancestor (LUCA) (Fig. 7, pathway 1) [14, 48]. Once in the chromosome, the intein may become fixed in a population and vertically transmit from generation to generation, as is the case for VMA1i (Fig. 6a). Some nDNA inteins may also be mobilized intragenomically (Fig. 7, pathway 2), as is the case for CHSi and GLTi, an event that has been postulated before due to their high sequence similarity (Fig. 6a and b) . Surprisingly, a putative case of intragenomic transfer is also found within the chloroplast genome, suggested by strong relatedness of cpDNA DdRPi and ClpPi (Fig. 6a). To our knowledge, this is the first example of intragenomic intein transfer within an organellar genome.
The largest group of cpDNA inteins, DnaBi, was likely inherited via endosymbiosis from a cyanobacterial progenitor (Fig. 7, pathway 3). Previously, a cyanobacterial DnaBi was shown to have weak sequence similarity with cpDNA DnaBi from a red alga . In the present study, a phylonetwork was reconstructed based on a joint dataset of both cpDNA DnaBi and bacterial DnaBi, and indicates that some of these inteins are distantly related (Fig. 6c). Endosymbiotic acquisition of inteins from bacteria seems reasonable and represents an unprecedented example of transfer of inteins in eukaryotes.
Finally, extragenomic horizontal gene transfer seems to be pervasive and one of the key forces driving eukaryotic intein distribution and diversity (Fig. 7, pathway 4). This is not surprising, given that viruses have been shown to transfer inteins among prokaryotes [23, 24, 27], and research has shown that a family of giant viruses has an ongoing process of exchanging inteins among each other . In algal-like cells (Fig. 7), various instances of intein transmission seem to play a role in the evolution of the current landscape. Here, we again see vDNA inteins clustering with cpDNA inteins (Fig. 6a), suggesting that spread involved viral transfer at some point. Algae prey on bacteria and archaea, and are infected by viruses, providing a place for inteins to transfer directly from bacteria to eukaryotic viruses . The fact that vDNA inteins are found in a class of viruses that replicate in the cytoplasm may allow the opportunity for spread.
Eukaryotic inteins are scarce compared to their prokaryotic counterparts, although they mimic patterns of distribution that suggest functional importance. These eukaryotic inteins are present in three genome types, nuclear, chloroplast and viral, and these inteins appear to have distinct routes of acquisition. Vertical transmission was primarily observed for fungal nuclear inteins, indicating maintenance since the last common ancestor. HEN-based intragenomic transfer allows these inteins to move to new sites within nuclei, as well as within chloroplasts. Ancient acquisition is also suggested by intein inheritance through endosymbiosis in chloroplast genomes. Finally, horizontal gene transfer spreads inteins across all three genome types, likely mediated by unique eukaryotic viruses that replicate in the cytoplasm. This expansive analysis garnered insight into intein dissemination, and will guide future experiments to investigate intein function and their use as novel drug targets.
The first set of full-length precursors and intein protein sequences for eukaryotes was collected from the intein database, InBase , and the National Center for Biotechnology Information (NCBI) Protein Database (www.ncbi.nlm.nih.gov/protein/) using previously described pipeline [27, 31]. Next, NCBI BLASTp  (blast.ncbi.nlm.nih.gov/Blast.cgi) was used to source additional, mostly unannotated inteins: known intein sequences were used as queries in a series of BLASTp searches against a non-redundant protein database limited to Eukaryota (taxid:2759). The results of the intein search in eukaryotic genomes and other metadata are provided in Additional file 2: Table S1. A similar approach was used to search for viral inteins. However, in addition to mining from NCBI Protein Database, viral inteins were investigated among entries in viral reference proteomes available at ViralZone  (https://viralzone.expasy.org/). The results of the intein search in genomes of viruses of eukaryotes are available in Additional file 3: Table S2.
To assess if genomes of pathogenic fungal species were overrepresented in the analysis, fungal species listed on the MycoCosm web portal (http://jgi.doe.gov/data-and-tools/mycocosm/) were examined (February 20, 2017). For each species, a ‘pathogen’ was defined as any species causing plant or human disease. Anything listed as a ‘parasite’, whether of nematodes, plants, or other fungi, were also considered pathogenic for the purpose of this analysis. Any organisms listed as ‘rarely causing disease’ were classified as pathogens. Fungi with unknown medical relevance were not ranked as pathogens.
The primary sequence data obtained for precursor as a whole and intein(s) separately were analyzed further using the following resources. The presence of conserved protein domains and motifs was verified using the Conserved Domain Database (CDD; www.ncbi.nlm.nih.gov/cdd) and Conserved Domain Search Service (CD Search; www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) as well as InterPro protein analysis tool (www.ebi.ac.uk/interpro/). KOG annotation for nuclear (nDNA) and viral (vDNA) inteins was performed using protein function annotation by the KOG database tool available at WebMGA server [44, 45] (weizhong-lab.ucsd.edu/metagenomic-analysis/); COG annotation for inteins from chloroplast genomes (cpDNA) was performed using the COG database  (www.ncbi.nlm.nih.gov/COG/). The endonuclease size sorting and plotting was done using customized Python scripts, which are available upon request. The CD search tool was utilized to identify which type of endonuclease the HEN(+) and HEN(+)extra inteins encoded. Splicing blocks A, B, F, and G were annotated based on classical classification based on conserved residues [50, 79].
Pairwise and multiple sequence alignments of inteins were performed using Clustal Omega (www.ebi.ac.uk/Tools/msa/clustalo/) unless indicated otherwise. The sequence similarity network (SSN) was generated based on full-length protein sequences of all identified inteins using EFI-EST tool  (Additional file 5) (http://efi.igb.illinois.edu/efi-est/) and visualized using Cytoscape . The phylonetwork was built using an alignment of cpDNA DnaBi, ClpPi, and bacterial DnaBi and uploaded into SplitsTree using default settings.
Secondary structure modeling and insertion site conservation
Intein-containing exteins with known structures were accessed from the Protein Data Bank (www.rcsb.org/pdb/home/home.do). They had the following PDB numbers: Prp8 – 5GMK, VMA1 – 3J9T, ThrRS – 3UGQ, and GLT – 1EA0. The structures have a conserved intein insertion site, despite being from inteinless hosts. Each structure file was truncated to include the domain, or portion of the protein, with intein insertion. All shortened structure PDB files were uploaded to Pro-Origami software (http://munk.csse.unimelb.edu.au/pro-origami/) . For conservation, a single full-length extein sequence was chosen as a query and uploaded to ConSurf . An automatic multiple sequence alignment was performed using MAFFT and homologs were collected from UNIREF90 using a maximum percent identity of 95 and a minimum percent identity of 60. A PDF of the query sequence alone colored according to conservation was downloaded and used for further analysis.
ATP-dependent Clp protease, proteolytic subunit
Clusters of orthologous genes
DNA-directed DNA polymerase
DNA-directed RNA polymerase
- IF2 eIF5B:
Eukaryotic translation initiation factor 5B
Eukaryotic orthologous groups
Last universal common ancestor
Mismatch repair protein
Nucleocytoplasmic large DNA viruses
Protein data bank
pre-mRNA processing factor 8
Sequence similarity network
Threonyl tRNA synthetase
Vacuolar membrane ATPase
Saleh L, Perler FB. Protein splicing in cis and in trans. Chem Rec. 2006;6:183–93.
Novikova O, Topilina N, Belfort M. Enigmatic distribution, evolution, and function of inteins. J Biol Chem. 2014;289:14490–7.
Gimble FS, Thorner J. Homing of a DNA endonuclease gene by meiotic gene conversion in Saccharomyces cerevisiae. Nature. 1992;357:301–6.
Liu XQ. Protein-splicing intein: genetic mobility, origin, and evolution. Annu Rev Genet. 2000;34:61–76.
Hall TM, Porter JA, Young KE, Koonin EV, Beachy PA, Leahy DJ. Crystal structure of a hedgehog autoprocessing domain: homology between hedgehog and self-splicing proteins. Cell. 1997;91:85–97.
Burglin TR. The hedgehog protein family. Genome Biol. 2008;9:241.
Burglin TR. Evolution of hedgehog and hedgehog-related genes, their origin from hog proteins in ancestral eukaryotes and discovery of a novel hint motif. BMC Genomics. 2008;9:127.
Perler FB. Protein splicing of inteins and hedgehog autoproteolysis: structure, function, and evolution. Cell. 1998;92:1–4.
Hirata R, Ohsumk Y, Nakano A, Kawasaki H, Suzuki K, Anraku Y. Molecular structure of a gene, VMA1, encoding the catalytic subunit of H(+)-translocating adenosine triphosphatase from vacuolar membranes of Saccharomyces cerevisiae. J Biol Chem. 1990;265:6726–33.
Kane PM, Yamashiro CT, Wolczyk DF, Neff N, Goebl M, Stevens TH. Protein splicing converts the yeast TFP1 gene product to the 69-kD subunit of the vacuolar H(+)-adenosine triphosphatase. Science. 1990;250:651–7.
Butler MI, Goodwin TJ, Poulter RT. A nuclear-encoded intein in the fungal pathogen Cryptococcus neoformans. Yeast. 2001;18:1365–70.
Butler MI, Gray J, Goodwin TJ, Poulter RT. The distribution and evolutionary history of the PRP8 intein. BMC Evol Biol. 2006;6:42.
Goodwin TJ, Butler MI, Poulter RT. Multiple, non-allelic, intein-coding sequences in eukaryotic RNA polymerase genes. BMC Biol. 2006;4:38.
Poulter RT, Goodwin TJ, Butler MI. The nuclear-encoded inteins of fungi. Fungal Genet Biol. 2007;44:153–79.
Chong S, Xu MQ. Protein splicing of the Saccharomyces cerevisiae VMA intein without the endonuclease motifs. J Biol Chem. 1997;272:15587–90.
Elleuche S, Nolting N, Poggeler S. Protein splicing of PRP8 mini-inteins from species of the genus Penicillium. Appl Microbiol Biotechnol. 2006;72:959–67.
Pearl EJ, Bokor AA, Butler MI, Poulter RT, Wilbanks SM. Preceding hydrophobic and beta-branched amino acids attenuate splicing by the CnePRP8 intein. Biochim Biophys Acta. 2007;1774:995–1001.
Bokor AA, van Kan JA, Poulter RT. Sexual mating of Botrytis cinerea illustrates PRP8 intein HEG activity. Fungal Genet Biol. 2010;47:392–8.
Moure CM, Gimble FS, Quiocho FA. Crystal structure of the intein homing endonuclease PI-SceI bound to its recognition sequence. Nat Struct Biol. 2002;9:764–70.
Wang S, Liu XQ. Identification of an unusual intein in chloroplast ClpP protease of Chlamydomonas eugametos. J Biol Chem. 1997;272:11869–73.
Douglas SE, Penny SL. The plastid genome of the cryptophyte alga, Guillardia theta: complete sequence and conserved synteny groups confirm its common ancestry with red algae. J Mol Evol. 1999;48:236–44.
Pietrokovski S. Identification of a virus intein and a possible variation in the protein-splicing reaction. Curr Biol. 1998;8:R634–5.
Nagasaki K, Shirai Y, Tomaru Y, Nishida K, Pietrokovski S. Algal viruses with distinct intraspecies host specificities include identical intein elements. Appl Environ Microbiol. 2005;71:3599–607.
Ogata H, Raoult D, Claverie JM. A new example of viral intein in Mimivirus. Virol J. 2005;2:8.
Bigot Y, Piegu B, Casteret S, Gavory F, Bideshi DK, Federici BA. Characteristics of inteins in invertebrate iridoviruses and factors controlling insertion in their viral hosts. Mol Phylogenet Evol. 2013;67:246–54.
Culley AI, Asuncion BF, Steward GF. Detection of inteins among diverse DNA polymerase genes of uncultivated members of the Phycodnaviridae. ISME J. 2009;3:409–18.
Kelley DS, Lennon CW, Belfort M, Novikova O. Mycobacteriophages as incubators for intein dissemination and evolution. MBio. 2016;7:e01537–16.
Perler FB. InBase: the intein database. Nucleic Acids Res. 2002;30:383–4.
Swithers KS, Soucy SM, Lasek-Nesselquist E, Lapierre P, Gogarten JP. Distribution and evolution of the mobile vma-1b intein. Mol Biol Evol. 2013;30:2676–87.
Soucy SM, Fullmer MS, Papke RT, Gogarten JP. Inteins as indicators of gene flow in the halobacteria. Front Microbiol. 2014;5:299.
Novikova O, Jayachandran P, Kelley DS, Morton Z, Merwin S, Topilina NI, Belfort M. Intein clustering suggests functional importance in different domains of life. Mol Biol Evol. 2016;33:783–99.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Hulo C, de Castro E, Masson P, Bougueleret L, Bairoch A, Xenarios I, Le Mercier P. ViralZone: a knowledge resource to understand virus diversity. Nucleic Acids Res. 2011;39:D576–82.
Fernandes JA, Prandini TH, Castro MD, Arantes TD, Giacobino J, Bagagli E, Theodoro RC. Evolution and application of inteins in Candida species: a review. Front Microbiol. 2016;7:1585.
Frederick RD, Snyder KE, Tooley PW, Berthier-Schaad Y, Peterson GL, Bonde MR, Schaad NW, Knorr DA. Identification and differentiation of Tilletia indica and T. walkeri using the polymerase chain reaction. Phytopathology. 2000;90:951–60.
Olson DH, Aanensen DM, Ronnenberg KL, Powell CI, Walker SF, Bielby J, Garner TW, Weaver G, Fisher MC. Mapping the global emergence of Batrachochytrium dendrobatidis, the amphibian chytrid fungus. PLoS One. 2013;8:e56802.
Aherfi S, Colson P, La Scola B, Raoult D. Giant viruses of amoebas: an update. Front Microbiol. 2016;7:349.
Morrow CA, Fraser JA. Sexual reproduction and dimorphism in the pathogenic basidiomycetes. FEMS Yeast Res. 2009;9:161–77.
Findley K, Rodriguez-Carres M, Metin B, Kroiss J, Fonseca A, Vilgalys R, Heitman J. Phylogeny and phenotypic characterization of pathogenic Cryptococcus species and closely related saprobic taxa in the Tremellales. Eukaryot Cell. 2009;8:353–61.
Visagie CM, Hirooka Y, Tanney JB, Whitfield E, Mwange K, Meijer M, Amend AS, Seifert KA, Samson RA. Aspergillus, Penicillium and Talaromyces isolated from house dust samples collected around the world. Stud Mycol. 2014;78:63–139.
Liu XZ, Wang QM, Goker M, Groenewald M, Kachalkin AV, Lumbsch HT, Millanes AM, Wedin M, Yurkov AM, Boekhout T, Bai FY. Towards an integrated phylogenetic classification of the Tremellomycetes. Stud Mycol. 2015;81:85–147.
Kocsube S, Perrone G, Magista D, Houbraken J, Varga J, Szigeti G, Hubka V, Hong SB, Frisvad JC, Samson RA. Aspergillus is monophyletic: evidence from multiple gene phylogenies and extrolites profiles. Stud Mycol. 2016;85:199–213.
Koonin EV, Fedorova ND, Jackson JD, Jacobs AR, Krylov DM, Makarova KS, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, et al. A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes. Genome Biol. 2004;5:R7.
Nordberg H, Cantor M, Dusheyko S, Hua S, Poliakov A, Shabalov I, Smirnova T, Grigoriev IV, Dubchak I. The genome portal of the department of energy joint genome institute: 2014 updates. Nucleic Acids Res. 2014;42:D26–31.
Galperin MY, Makarova KS, Wolf YI, Koonin EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 2015;43:D261–9.
Apgar J, Ross M, Zuo X, Dohle S, Sturtevant D, Shen B, de la Vega H, Lessard P, Lazar G, Raab RM. A predictive model of intein insertion site for use in the engineering of molecular switches. PLoS One. 2012;7:e37355.
Poland BW, Xu MQ, Quiocho FA. Structural insights into the protein splicing mechanism of PI-SceI. J Biol Chem. 2000;275:16408–13.
Pietrokovski S. Intein spread and extinction in evolution. Trends Genet. 2001;17:465–72.
Swithers KS, Senejani AG, Fournier GP, Gogarten JP. Conservation of intron and intein insertion sites: implications for life histories of parasitic genetic elements. BMC Evol Biol. 2009;9:303.
Perler FB, Olsen GJ, Adam E. Compilation and analysis of intein sequences. Nucleic Acids Res. 1997;25:1087–93.
Monier A, Sudek S, Fast NM, Worden AZ. Gene invasion in distant eukaryotic lineages: discovery of mutually exclusive genetic elements reveals marine biodiversity. ISME J. 2013;7:1764–74.
Ashkenazy H, Abadi S, Martz E, Chay O, Mayrose I, Pupko T, Ben-Tal N. ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 2016;44:W344–50.
Koufopanou V, Goddard MR, Burt A. Adaptation for horizontal transfer in a homing endonuclease. Mol Biol Evol. 2002;19:239–46.
de La Fuente J, Garcia-Garcia JC, Blouin EF, Rodriguez SD, Garcia MA, Kocan KM. Evolution and function of tandem repeats in the major surface protein 1a of the ehrlichial pathogen Anaplasma marginale. Anim Health Res Rev. 2001;2:163–73.
Djinovic-Carugo K, Gautel M, Ylanne J, Young P. The spectrin repeat: a structural platform for cytoskeletal protein assemblies. FEBS Lett. 2002;513:119–23.
Staub E, Perez-Tur J, Siebert R, Nobile C, Moschonas NK, Deloukas P, Hinzmann B. The novel EPTP repeat defines a superfamily of proteins implicated in epileptic disorders. Trends Biochem Sci. 2002;27:441–4.
Atkinson HJ, Morris JH, Ferrin TE, Babbitt PC. Using sequence similarity networks for visualization of relationships across diverse protein superfamilies. PLoS One. 2009;4:e4345.
Bapteste E, Lopez P, Bouchard F, Baquero F, McInerney JO, Burian RM. Evolutionary analyses of non-genealogical bonds produced by introgressive descent. Proc Natl Acad Sci. 2012;109:18266–72.
Bapteste E, van Iersel L, Janke A, Kelchner S, Kelk S, McInerney JO, Morrison DA, Nakhleh L, Steel M, Stougie L, Whitfield J. Networks: expanding evolutionary thinking. Trends Genet. 2013;29:439–41.
Gerlt JA, Bouvier JT, Davidson DB, Imker HJ, Sadkhin B, Slater DR, Whalen KL. Enzyme function initiative-enzyme similarity tool (EFI-EST): a web tool for generating protein sequence similarity networks. Biochim Biophys Acta. 2015;1854:1019–37.
Koufopanou V, Burt A. Degeneration and domestication of a selfish gene in yeast: molecular evolution versus site-directed mutagenesis. Mol Biol Evol. 2005;22:1535–8.
Pietrokovski S. A new intein in cyanobacteria and its significance for the spread of inteins. Trends Genet. 1996;12:287–8.
Archibald JM, Keeling PJ. Recycled plastids: a 'green movement' in eukaryotic evolution. Trends Genet. 2002;18:577–84.
Dorrell RG, Howe CJ. What makes a chloroplast? Reconstructing the establishment of photosynthetic symbioses. J Cell Sci. 2012;125:1865–75.
Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:254–67.
Shah NH, Muir TW. Inteins: nature's gift to protein chemists. Chem Sci. 2014;5:446–61.
Elleuche S, Poggeler S. Inteins, valuable genetic elements in molecular biology and biotechnology. Appl Microbiol Biotechnol. 2010;87:479–89.
Gramespacher JA, Stevens AJ, Nguyen DP, Chin JW, Muir TW. Intein zymogens: conditional assembly and splicing of split inteins via targeted proteolysis. J Am Chem Soc. 2017;139:8074–7.
Callahan BP, Topilina NI, Stanger MJ, Van Roey P, Belfort M. Structure of catalytically competent intein caught in a redox trap with functional and evolutionary implications. Nat Struct Mol Biol. 2011;18:630–3.
Topilina NI, Green CM, Jayachandran P, Kelley DS, Stanger MJ, Piazza CL, Nayak S, Belfort M. SufB intein of Mycobacterium tuberculosis as a sensor for oxidative and nitrosative stresses. Proc Natl Acad Sci. 2015;112:10348–53.
Topilina NI, Novikova O, Stanger M, Banavali NK, Belfort M. Post-translational environmental switch of RadA activity by extein-intein interactions in protein splicing. Nucleic Acids Res. 2015;43:6631–48.
Ciragan A, Aranko AS, Tascon I, Iwai H. Salt-inducible protein splicing in cis and trans by inteins from extremely halophilic archaea as a novel protein-engineering tool. J Mol Biol. 2016;428:4573–88.
Lennon CW, Stanger M, Belfort M. Protein splicing of a recombinase intein induced by ssDNA and DNA damage. Genes Dev. 2016;30:2663–8.
Reitter JN, Cousin CE, Nicastri MC, Jaramillo MV, Mills KV. Salt-dependent conditional protein aplicing of an intein from Halobacterium salinarum. Biochemistry. 2016;55:1279–82.
Zhang L, Zheng Y, Callahan B, Belfort M, Liu Y. Cisplatin inhibits protein splicing, suggesting inteins as therapeutic targets in mycobacteria. J Biol Chem. 2011;286:1277–82.
Chan H, Pearson CS, Green CM, Li Z, Zhang J, Belfort G, Shekhtman A, Li H, Belfort M. Exploring intein inhibition by platinum compounds as an antimicrobial strategy. J Biol Chem. 2016;291:22661–70.
Stop neglecting fungi. Nat Microbiol. 2017;2:17120. https://www.ncbi.nlm.nih.gov/pubmed/?term=stop+neglecting+fungi.
Clerissi C, Grimsley N, Desdevises Y. Genetic exchanges of inteins between prasinoviruses (phycodnaviridae). Evolution. 2013;67:18–33.
Pietrokovski S. Modular organization of inteins and C-terminal autocatalytic domains. Protein Sci 1998;7:64-71.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Stivala A, Wybrow M, Wirth A, Whisstock JC, Stuckey PJ. Automatic generation of protein structure cartoons with pro-origami. Bioinformatics. 2011;27:3315–6.
The authors acknowledge Danielle Kelley and Dr. Valjean R. Bacot-Davis for useful discussions regarding the analysis and figures. They also wish to thank Dr. Nicholas Schiraldi for data analytics and valuable coding.
Work in our laboratory is supposed by National Institutes of Health grants GM39422 and GM44844 to MB.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary files. Additional information is available from the corresponding author upon request.
CMG, ON, and MBconceived the study; CMG performed data mining and CMG and ON performed bioinformatics analyses. CMG, ON, and MB prepared the manuscript and approved the final version. CMG and ON rendered figures.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Tables S3-S13 and Supplementary Figures 1–6. This file contains Supplementary Tables S3–13 and Supplementary Figures 1–6 with legends. (PDF 1447 kb)
Inteins in nDNA and cpDNA. This is a spreadsheet of eukaryotic inteins in nDNA and cpDNA with accession numbers and extein/intein sequences mined in this study. (XLSX 177 kb)
Inteins in vDNA. This is a spreadsheet of eukaryotic inteins in vDNA with accession numbers and extein/intein sequences mined in this study. (XLSX 80 kb)
DnaBi alignment. This file contains an alignment of cpDnaBi, bacterial DnaBi, and cpClpPi for Fig. 6c. (FA 57 kb)
Sequences of all proteins. This file contains all nDNA, cpDNA, and vDNA inteins and HH proteins used in analyses in this study. (TXT 152 kb)
About this article
Cite this article
Green, C.M., Novikova, O. & Belfort, M. The dynamic intein landscape of eukaryotes. Mobile DNA 9, 4 (2018). https://doi.org/10.1186/s13100-018-0111-x