Casposons: a new superfamily of self-synthesizing DNA transposons at the origin of prokaryotic CRISPR-Cas immunity
© Krupovic et al.; licensee BioMed Central Ltd. 2014
Received: 18 April 2014
Accepted: 13 May 2014
Published: 19 May 2014
Diverse transposable elements are abundant in genomes of cellular organisms from all three domains of life. Although transposons are often regarded as junk DNA, a growing body of evidence indicates that they are behind some of the major evolutionary innovations. With the growth in the number and diversity of sequenced genomes, previously unnoticed mobile elements continue to be discovered.
We describe a new superfamily of archaeal and bacterial mobile elements which we denote casposons because they encode Cas1 endonuclease, a key enzyme of the CRISPR-Cas adaptive immunity systems of archaea and bacteria. The casposons share several features with self-synthesizing eukaryotic DNA transposons of the Polinton/Maverick class, including terminal inverted repeats and genes for B family DNA polymerases. However, unlike any other known mobile elements, the casposons are predicted to rely on Cas1 for integration and excision, via a mechanism similar to the integration of new spacers into CRISPR loci. We identify three distinct families of casposons that differ in their gene repertoires and evolutionary provenance of the DNA polymerases. Deep branching of the casposon-encoded endonuclease in the Cas1 phylogeny suggests that casposons played a pivotal role in the emergence of CRISPR-Cas immunity.
The casposons are a novel superfamily of mobile elements, the first family of putative self-synthesizing transposons discovered in prokaryotes. The likely contribution of capsosons to the evolution of CRISPR-Cas parallels the involvement of the RAG1 transposase in vertebrate immunoglobulin gene rearrangement, suggesting that recruitment of endonucleases from mobile elements as ready-made tools for genome manipulation is a general route of evolution of adaptive immunity.
Cellular organisms in the three domains of life are under constant onslaught of invading mobile genetic elements (MGE), such as transposons, viruses and plasmids. Many, if not most, of these diverse selfish elements insert into the chromosomes of the cellular hosts, either as an obligate part of their life cycles or at least occasionally, and in multicellular eukaryotes constitute a substantial proportion of the host genome. For example, sequencing of the human genome has shown that transposons or relics thereof amount to 35% to 50% of the genome [1, 2] and subsequent analyses have only revised these estimates upward [3, 4]. Even more strikingly, in some green plants, MGE-derived DNA seems to represent more than 70% of the genome [5, 6]. Although not as abundant as in eukaryotes, proviruses and other MGE constitute up to 30% of some bacterial genomes [7, 8]. The effects of MGE integration vary from beneficial (gain of new phenotypic traits, such as antibiotic resistance or toxin production) to deleterious (disruption or inactivation of essential cellular genes upon MGE insertion) [7, 9–11]. For most prokaryotic plasmids and viruses, the circular form of the MGE genome is inserted into specific loci (site-specific integration) of the cellular chromosome with the aid of MGE-encoded enzymes known as integrases . The integrases are grouped into two major families on the basis of sequence conservation and mechanistic relatedness: (1) tyrosine recombinases use a catalytic tyrosine residue which attacks the DNA and becomes covalently linked to it during strand exchange; and (2) serine recombinases for the same purpose use a nucleophilic serine residue .
Transposons are DNA segments that move from one location in the host genome to another. Although several classification schemes have been proposed, based on the nature of the transposition intermediate, transposons can be generally grouped into two classes [10, 13] or types . Class I (or type 2) elements—also known as retrotransposons—transpose via an RNA intermediate which prior to integration is copied back to the DNA form by the element-encoded reverse transcriptase. Class II (or type 1) DNA transposons move in the genome by the so-called ‘cut-and-paste’ mechanism whereby the transposon is excised from its initial location and inserted into a new genomic locus. Most of the class II transposons possess characteristic terminal inverted repeats (TIR) but differ widely in terms of the transposases they encode, the specific mechanisms of transposition, the element size and gene content [10, 13, 15]. Although most transposases belong to the DDE superfamily (named after two aspartate and one glutamate residues that form the catalytic triad of these enzymes) [10, 13, 16], some transposons encode transposases homologous to the rolling-circle replication initiation endonucleases found in numerous viruses and plasmids [17–19], to phage integrase-like tyrosine recombinases [20, 21] or to the serine integrases/invertases . Furthermore, some bacterial and eukaryotic viruses encode transposases that are involved in the integration of the viral genome into the host chromosome, thereby partially blurring the distinction between different MGE types [23–25].
A distinct group of MGE consists of large (15 to 20 kb), self-synthesizing DNA transposons, called Mavericks or Polintons [26, 27]. The defining feature of Polintons/Mavericks is that they encode their own protein-primed type B DNA polymerase which is most likely involved in the transposon replication (hence ‘self-synthesizing’ transposons) . In addition, these transposons encode several hallmark viral proteins, the genome packaging ATPase and protease. Recently, we have shown that Polintons/Mavericks also encode major and minor capsid proteins, suggesting that these elements combine features of bona fide viruses and transposons . Polintons/Mavericks are widespread in diverse unicellular and multicellular eukaryotes [26, 27]. In contrast, no such self-synthesizing DNA transposons have been described in prokaryotes.
To survive the proliferation of various MGE and to maintain genetic integrity, cellular organisms have evolved numerous defense lines, including a variety of innate and adaptive immunity mechanisms [28–30]. Although once considered to be characteristic exclusively of animals, adaptive immunity has been recently discovered in bacteria and archaea [30–33]. This system consists of arrays of clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated proteins (Cas) and elicits interference against foreign nucleic acids by degrading them in a sequence-specific fashion. The specificity is ensured by the unique spacers homologous to viral or plasmid DNA and integrated into the CRISPR loci. The action of the CRISPR-Cas system can be divided into three stages. The first stage, called adaptation, involves insertion of foreign DNA spacers into the CRISPR repeats. This step is mediated by the two most conserved core proteins of the CRISPR-Cas system, Cas1 and Cas2 [34–36]. Although the mechanistic details of adaptation remain poorly understood, it has been demonstrated that Cas1 is the endonuclease responsible for the excision of the protospacer from the foreign DNA and its insertion into the CRISPR cassette [36–40]. During the second stage, expression and processing, the CRISPR locus containing the arrays of spacers is transcribed, producing a long pre-crRNA (CRISPR RNA), which is subsequently processed by Cas proteins into short guide crRNAs. The final stage is called interference and involves degradation of the alien DNA or RNA by the Cas enzymatic machinery guided by the bound crRNA [30, 31, 35]. Phylogenomic analyses of the Cas proteins from diverse archaea and bacteria yielded a wealth of information on the diversity and evolution of the CRISPR-Cas immunity . However, it remains unclear how this sophisticated defense system emerged in the first place.
Here, we describe the discovery and characterization of a new superfamily of MGE that possess several features resembling the eukaryotic self-synthesizing DNA transposons but are integrated in the genomes of various archaea and some bacteria. Along with family B DNA polymerases (PolB), that are related either to viral protein-primed PolBs or to typical archaeal PolBs, these elements, which we denote ‘Casposons’, encode Cas1 proteins of a distinct subfamily. We propose that, different from other known MGE, casposons utilize Cas1 endonucleases for integration into the host genomes via a mechanism resembling that of spacer integration by CRISPR-Cas systems. Given that Cas1 is a key enzyme of the CRISPR-Cas immunity and considering the deep branching of casposon homologs in the Cas1 phylogeny, casposons appear to have played a pivotal role in the origin of the adaptive immune system in prokaryotes.
Genomic islands containing stand-alone cas1genes
A recent comparative genomic survey of cas genes revealed two distinct groups of cas1 genes that are not associated with CRISPR loci and form two distinct clades in the Cas1 phylogeny (hereinafter ‘Cas1-solo’) . The first Cas1-solo group was exclusively found in members of the archaeal order Methanomicrobiales and did not show any evidence of horizontal gene transfer (HGT) whereas the second group displayed a more patchy distribution. Most of the group 2 members were from the euryarchaeal class Methanomicrobia; however, several representatives were also detected in members of Thaumarchaeota as well as in the hyperthermophilic euryarchaeon Aciduliprofundum boonei affiliated with the order Thermoplasmatales. We hypothesized that Cas1-solo might be ancestral to the Cas1 proteins found in CRISPR-Cas systems and set out to investigate their provenance and potential function.
Previous site-directed mutagenesis study has identified four conserved residues constituting the active site of Cas1 endonucleases (E141, H208, D218 and D221 in Escherichia coli Cas1); alanine substitutions at any of these positions abolished the nuclease activity of the Cas1 against all substrates tested . Examination of the multiple alignment of Cas1-solo protein sequences showed that none of the group 1 members had the full complement of active site residues [see Additional file 1: Figure S1a], indicating that these proteins are unlikely to be active endonucleases (or less likely, rely on a different set of catalytic residues) and evolved under different constraints than the functional Cas1 proteins. In contrast, the four catalytic residues are strictly conserved in all Cas1-solo proteins from group 2 [see Additional file 1: Figure S1b]. Therefore, for further analysis we focused on group 2 of Cas1-solo.
Discovery of casposons
Gene content analysis showed that Cas1-solo-encoding genomic islands from Thaumarchaeota contain PolBs that belong to the group of protein-primed polymerases. These polymerases are encoded by various viruses and eukaryotic self-synthesizing transposons of the Polinton/Maverick family [26, 27] but generally not by cellular organisms. Thus, we hypothesized that these islands represent integrated MGE, analogous to the eukaryotic self-synthesizing transposons. DNA transposons typically possess TIR and upon integration into the genome often contain a specific mark, the target site duplication (TSD) which flanks the transposon [13, 15]. We investigated the Cas1 and PolB-containing genomic islands for the presence of these features and found that in nearly all cases these loci were flanked by TIRs and direct repeats which correspond to TSD [see Additional file 1: Table S1]. None of these elements contained identifiable genes for serine or tyrosine recombinases nor did they carry conserved transposase genes (see also below). The only enzyme that is consistently present in all these elements and, judged by its experimentally characterized activity, is capable of mediating the integration of the elements into the host genome is Cas1. Accordingly, we denote this new group of transposon-like elements ‘Casposons’. The conservation of polB genes places casposons as a new (super) family into the class of self-synthesizing large DNA transposons .
TIRs, TSDs and integration sites
The unique casposon TIRs are highly variable in length (25 to 602 nucleotides, median of 56) and could be identified in all casposons, except for the three closely related elements (MetBur-C1 to C3) in the genome of Methanococcoides burtonii DSM 6242 [see Additional file 1: Figure S2]. Some of the TIRs contain internal palindromic sequences [see Additional file 1: Figure S2}.
The TSDs result from the fill-in repair of staggered nicks introduced by transposases at the target site upon insertion of MGE [15, 41]. The length of the TSD differs depending on the transposase involved but in addition varies within as well as between transposon families [13, 15]. The great majority of casposons are flanked by perfect direct repeats corresponding to TSD and ranging in length from 1 to 27 nucleotides (median of 15; Additional file 1: Table S1). In a substantial fraction of the identified casposons, one or both TIRs partially overlap with the TSDs [see Additional file 1: Table S1 and Figure S2], suggesting that, prior to integration, these casposons contained short terminal overhangs that were partially complementary to the staggered ends of the nicked target site. By contrast, the casposons in which the overlaps between the TSDs and TIRs were not present likely had blunt termini prior to integration.
Most transposons do not display strong target site preference but some are known to integrate site-specifically [42, 43]. Ten casposons were found to be inserted into intergenic loci whereas for eight others, the target sites overlapped with coding sequences. The target sites of the three complete thaumarchaeal casposons were located within the 3′-distal region of the gene encoding the translation elongation factor aEF-2, whereas in five euryarchaeal casposons the target site overlapped four to seven 3′-distal nucleotides of different tRNA genes [see Additional file 1: Table S1]. Notably, eukaryotic transposons of the recently described DADA superfamily integrate site-specifically into snRNA and tRNA genes . However, unlike DADA transposons, which integrate close to the anticodon loop of tRNA genes, casposons do not alter the sequence of their target genes (either tRNA or aEF-2) and are located proximal to these genes. This pattern of integration is reminiscent of the bacterial Tn7 transposon which recognizes the 3′-distal region of the highly conserved glutamine synthetase gene (glmS) but inserts downstream of the glmS coding region, preserving the integrity of the latter . Such a strategy ensures that the integration of Tn7 and casposons does not disrupt genes essential for host viability, thereby ensuring successful propagation of both the host and the respective MGE.
In most cases, when complete genome sequences are available, casposons are present in one copy per genome, consistent with their site-specific integration. However, M. burtonii DSM 6242 encompasses three closely related casposons (MetBur-C1 to -C3) which are adjacent on the genome [see Additional file 1: Figure S3a], suggesting recent activity of casposons in this archaeon. Notably, MetBur-C3 appears to be inactivated because two of the genes in this element contain amber mutations. Again, a parallel with the Tn7 transposon can be drawn. Tn7 is usually present in a single copy per genome. However, with lower efficiency, additional copies can be integrated into the same target site, forming islands of tandem transposons, some of which become inactivated with time . M. psychrophilus R15 is another organism in which remnants of a second, adjacent casposon are present [see Additional file 1: Figure S3b]. The TIRs of the latter element could not be identified, suggesting that it might be in the process of deterioration. The patchy taxonomic distribution of casposons in archaeal and bacterial genomes as well as their amplification in certain organisms suggests that they are active mobile elements. However, the possibility that the amplification of casposons in M. burtonii DSM 6242 and M. psychrophilus R15 genomes is a result of segmental duplication cannot be ruled out. Experimental study of casposon integration and excision, as well as analysis of many more complete archaeal genomes, is necessary to provide definitive answers regarding the mobility of these MGE.
Classification of casposons
Family 2 casposons are present in diverse members of the archaeal phylum Euryarchaeota [see Additional file 1: Table S1], including the unclassified human gut-associated methanogen Methanomassiliicoccus luminyensis B10  as well as the hyperthermoacidophile A. boonei T469 . PolBs encoded by casposons of family 2 are related to the PolB3 family of typical archaeal RNA-primed DNA polymerases [see Additional file 1: Figure S4] [52, 53]. The casposon polymerases form a sister group to a small clade of PolBs from Thermoproteales (phylum Crenarchaeota). Notably, in the latter clade the polB gene of Ignisphaera aggregans DSM 17230 is located within an integrated mobile element which is unrelated to casposons and carries a gene for a tyrosine integrase (KSM, MK, EVK, unpublished work). This observation suggests that, as with the family 1 casposons, there could have been exchange of PolB genes between family 2 casposons and other types of MGE. Cas1 proteins of all family 2 casposons contain a conserved C-terminal fusion of an HTH domain which is not found in any other Cas1 proteins. Notably, a similar HTH domain is also found in the C-termini of PolBs of family 2 casposons. Although family 2 casposons vary considerably in size (6 to 16 kb) and gene content, most of them share a core of five genes which encode Cas1, PolB, an HNH endonuclease and two distinct HTH proteins (Figure 2). One of the conserved HTH proteins contains a C-terminal HEAT repeat domain (PF02985); HEAT repeats form rod-like helical structures that mediate various protein-protein interactions . The conserved HTH proteins and the HNH endonuclease might be involved in the target site recognition and subsequent casposon integration, in concert with Cas1. This mechanism of integration would resemble integration of the site-specific transposon Tn7 mentioned above. The heterotrimeric Tn7 transposase TnsABC binds the termini of the transposon whereas targeting to the specific region of the host glmS gene is mediated by the sequence-specific DNA-binding protein TnsD .
Family 3 casposons are present in the genomes of different bacteria, including an uncultivated thermophilic bacterium Candidatus ‘Acetothermum autotrophicum’ [see Additional file 1: Table S1]. In the Cas1 phylogeny, Family 3 casposons form a distinct clade that is a sister group to the rest of the casposons (Figure 1). By contrast, in the PolB tree, the Family 3 clade emerges from within the Family 2 casposons [see Additional file 1: Figure S4], compatible with the possibility that casposons emerged in archaea and were horizontally transferred to bacteria subsequent to the divergence of the casposon families 1 and 2. Cas1 protein of NitAL212-C1 contains a zinc-binding domain (ZBD) and a HTH domain fused to the N- and C-termini of the Cas1 domain, respectively, whereas in AceAut-C1 both ZBD and HTH are fused to the C-terminus of the Cas1 domain (Figure 2). The Cas1 proteins of StrAlb-C1 and HenMar-C1 do not contain any additional domains, similar to the Cas1 of Family 1 casposons. Three of the four group 3 casposons contain genes for the HNH endonuclease and a conserved HTH protein shared with the group 2 casposons (Figure 2). StrAlb-C1 encodes only a homolog of the HNH endonuclease although a gene for an unrelated HTH protein, which might be functionally equivalent, was also identified (Figure 2). Thus, both the PolB phylogeny and the comparison of the sets of predicted genes point to an affinity between the casposon families 2 and 3.
Casposon gene repertoire
Casposons vary greatly in terms of gene content, both within and between the three families described above, and carry many lineage-specific genes. Virtually all of the genes for which functions could be inferred are predicted to be involved in various DNA manipulations. Three consistent themes could be discerned among the products of the casposon genes (Figure 2 and Additional file 1: Table S2).
The first group of proteins includes predicted nucleases potentially involved in casposon integration/excision. In addition to the Cas1 endonuclease, the hallmark casposon protein, this group includes the HNH endonuclease that is present in all Family 2 and Family 3 casposons and is likely to cooperate with Cas1 in the integration and excision processes. MetLum-C1 and MetPsy-C1 contain genes for apurinic/apyrimidinic (AP) endonucleases, whereas MetPsy-C1 also encodes a GIY-YIG nuclease and a uracil-DNA glycosylase (UDG) that might be involved in the repair of the termini following casposon integration. MetTin-C1 encodes an exonuclease which could contribute to the processing of the casposon termini.
The second group of casposon proteins is implicated in DNA replication. Besides the two types of PolB genes, many casposons carry genes for various helicases which might assist during the replication of the casposon DNA. Notably, AceAut-C1 encodes not only a HerA-like helicase but also a putative DnaC-like helicase loader as well as an additional protein containing the AAA + ATPase domain, which is found in many helicases, including MCM . HenMar-C1 encodes a, so far, unique fusion protein containing an N-terminal nuclease domain related to the Cas4-like proteins of the CRISPR-Cas systems and a C-terminal superfamily I helicase domain.
The third category consists of various small DNA-binding proteins containing HTH, ZBD or ribbon-helix-helix (RHH) domains. Various combinations of these proteins are encoded in most casposons (Figure 2), and it cannot be ruled out that some of the uncharacterized small proteins contain highly derived versions of DNA-binding domains. These proteins could contribute to both integration/excision and replication of the casposons, and in addition, some of them might regulate expression of casposon and/or host genes.
In addition, casposons encode various enzymes some of which are associated with genome, RNA or chromatin modification whereas others are implicated in defense functions and metabolic processes [see Additional file 1: Table S2]. Three casposons encode highly divergent Sir2-like proteins; in Sulfolobus, Sir2 has been shown to deacetylate the major archaeal chromatin protein, Alba, thereby modulating the chromatin structure . Although the casposon-encoded Sir2-like proteins might also function as deacetylases, alternative enzymatic activities of these derived Sir2 homologs cannot be ruled out. By contrast, AmiAut-C1 encodes a unique fusion protein containing an N-terminal GCN5 acetyltransferase domain and a C-terminal queuine tRNA-ribosyltransferase domain. In addition, AciBoo-C1 and HenMar-C1 encode N6 and C5 DNA methyltransferases, respectively, whereas MetLum-C1 encodes a putative restriction endonuclease. Other notable proteins encoded by casposons include a KAP family P-loop ATPase  (MetMah-C1), Ser/Thr kinase (HenMar-C1) and a Sm-like RNA-binding protein (AmiAut-C1).
Two casposons, MetArv-C1 and NitAL212-C1, carry insertion sequence (IS) elements of the families ISNCY and IS1595, respectively . In both cases, the IS transposase genes are flanked by typical short inverted repeats and TSDs, indicating that the IS elements parasitize casposons rather than participate in their propagation. The sporadic conservation of functionally diverse genes in distinct casposons, even those that belong to the same family, indicates that, similar to viruses, casposons can horizontally acquire genes from various sources.
Transposons as a type of MGE are polyphyletic with respect to the enzymes mediating their transposition [10, 13, 15]. Here, we described a new type of mobile elements, the casposons, which appear to rely on Cas1-like endonucleases for genome integration. Structures of several Cas1 proteins have been solved [37, 39, 40] showing that Cas1 proteins adopt a novel structural fold, unrelated to the folds of any of the transposases described to date. Nevertheless, casposons share a number of features with known DNA transposons. On the one hand, they appear to behave as site-specific transposons, akin to the bacterial transposon Tn7 [43, 44]. On the other hand, the molecular structure of casposons is highly reminiscent of the eukaryotic self-synthesizing DNA transposons of the Polinton/Maverick superfamily [26, 27]. Similar to the Polintons/Mavericks, casposons possess TIRs and encode their own DNA polymerase genes which in family 1 casposons belong to the same, protein-primed class as the polymerase of Polintons/Mavericks.
Given the presence of a PolB gene in all casposons, we propose to classify these elements as the second superfamily within the class of self-synthesizing DNA transposons, in addition to the Polintons/Mavericks . Most of the polintons show a distinct virus-like character, with two conserved genes encoding major and minor capsid proteins, suggesting that these elements form virus particles under some circumstances and prompting their proposed re-classification as polintoviruses . The casposons, however, do not encode any detectable homologs of capsid proteins and accordingly are likely to adhere to the transposon lifestyle. We define casposons as self-synthesizing MGE which rely on Cas1 for integration. Under this definition, casposons may encode distinct DNA polymerases (as is indeed the case for family 1 compared to families 2 and 3) but should possess the ability to self-synthesize.
The Cas1 endonuclease is the key player in the adaptation step of the CRISPR-Cas immunity. Consistent with its importance, Cas1 is the most stable and conserved component of functional CRISPR-Cas systems and is considered a signature gene for these defense systems [32, 33, 35]. A model of spacer acquisition, which integrates the available experimental data, has been proposed [30, 34–36]. This model helps to predict the specific path of the Cas1-mediated integration of casposons which is also consistent with the detailed analysis of the terminal sequences of integrated casposons [see Additional file 1: Figure S2]. Figures 4b and 4c depict the parallel flows of events underlying the insertion of new CRISPR spacers and casposons, respectively. In both cases, staggered nicks are introduced into the target sequence which in the case of CRISPR-Cas corresponds to the first repeat proximal to the leader sequence. In the case of some casposons, the TIR-containing termini are processed to produce short overhangs complementary to the tips of the nicked target site (see above). In the next step, the ends of the protospacer/casposon are joined to those of the nicked target site. The observation that in CRISPR-Cas systems Cas1 is the only protein whose enzymatic activity is essential for integration of new spacers , a process that involves cutting and rejoining of the cellular DNA within the CRISPR repeat arrays, suggests that the casposon Cas1 also possesses both DNA cutting and joining activities. However, the latter activity of Cas1 remains to be demonstrated experimentally. Finally, the target site is fill-in repaired, completing the casposon/spacer insertion and resulting in the TSD (for casposons) or repeat duplication (in CRISPR-Cas) (Figure 4b, c).
The discovery of casposons has important evolutionary implications. The deep branching of casposon Cas1 homologs within the global Cas1 phylogeny (Figure 1) is compatible with the possibility that the Cas1 family of endonucleases emerged in the context of mobile elements and only later was adapted for cellular defense. Consequently, we propose that casposons played a pivotal role in the origin of prokaryotic CRISPR-Cas immunity. The origin of Cas1 appears not to be the only contribution of transposable elements to the evolution of CRISPR-Cas. Indeed, recent comparative genomic analysis of the type II CRISPR-Cas systems has shown that Cas9, the key protein of these systems involved in the RNA processing and interference stages, most likely, also evolved from a distinct class of transposon proteins .
It has been previously hypothesized that the CRISPR-Cas system originated in archaea , and the present observations on the likely archaeal origin of casposons appear compatible with this hypothesis. However, it cannot be ruled out that similar to some other MGE, casposons are even more ancient and antedate advanced cellular life forms .
Strikingly, transposons can also be placed at the root of adaptive immunity in eukaryotes. The RAG1 protein, which plays a central role during the V(D)J recombination, was derived from the DDE transposase of Transib transposons . The parallel contribution of transposons to the origin of adaptive immunity in prokaryotes and eukaryotes emphasizes that MGE are the molecular architects behind some of the major evolutionary innovations of their hosts, in particular, the cellular defense systems [23, 65]. More specifically, given the mechanistic similarity between MGE transposition and integration, on the one hand, and insertion of spacers by the CRISPR-Cas system and immunoglobulin gene rearrangement, on the other hand, integrases and transposases appear to be ready-made tools that can be recruited and utilized by adaptive immunity systems.
The diversity of MGE is astounding and is far from being fully explored. This state of affairs is well illustrated by the discovery of casposons described here. Casposons constitute the second superfamily of self-synthesizing transposon-like MGE, beside the eukaryotic Polinton/Maverick transposons, and are the first representatives of this class of elements in prokaryotes. Different MGE have evolved a number of unrelated molecular mechanisms to perform similar tasks that ensure their propagation within the host cells. The casposons, so far, are unique as the only group of MGE that apparently rely on Cas1 endonucleases, a key component of the prokaryotic CRISPR-Cas defense system, for insertion into and excision from the host genome. The perennial arms race between cellular organisms and various MGE seems to be one of the major driving forces underlying the evolution of both interacting parties and it is becoming increasingly clear that cells and MGE exchange molecular inventions that emerge in the process of this evolutionary struggle. The adaptive immunity of both prokaryotes and eukaryotes apparently evolved via recruitment of recombinases from distinct MGE, the casposons and the Transib family transposons, respectively. Although this route of evolution seems paradoxical given that MGE are the primary targets of the immunity systems, it is becoming clear that throughout the course of evolution, MGE have served as a rich source of naturally evolved tools for cellular genome engineering that had a major impact on the diversification of cellular organisms.
Casposon protein sequences were analyzed using PSI-BLAST , CD-Search , and HHpred . Inverted and direct repeats flanking the casposons were analyzed using Unipro UGENE . The palindromic repeats within the casposon TIR sequences were identified using Mfold . Insertion sequences were analyzed using ISfinfer . Multiple sequences alignments were built using Promals3D  and Muscle . The Polinton/Maverick PolB sequences were recovered from the Repbase Update database . For phylogenetic analysis, gapped columns (more than 30% of gaps) and columns with low information content were removed from the alignment . Phylogenetic analysis was carried out by using FastTree , with the Jones–Taylor–Thornton model of amino acid evolution and γ-CAT estimation of evolutionary rates across sites. The trees were visualized using MEGA6 . For the Cas1 phylogeny, Cas1 protein sequences representing all major types and subtypes of the CRISPR-Cas systems were obtained from  and supplemented with the casposon-encoded Cas1 protein sequences. The Cas1 sequence alignment used to generate the tree is provided in Additional file 2.
MK was partly supported by the European Molecular Biology Organization (ASTF 82–2014). EVK and KSM are supported by intramural funds of the US Department of Health and Human Services (to the National Library of Medicine). PF is supported by the European Union’s Seventh Framework Program (FP/2007-2013) / Project EVOMOBIL - ERC Grant Agreement no.340440. DP is supported by the Agence nationale de la recherche (ANR) program BLANC, project EXAVIR.
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, Funke R, Gage D, Harris K, Heaford A, Howland J, Kann L, Lehoczky J, LeVine R, McEwan P, McKernan K, Meldrim J, Mesirov JP, Miranda C, Morris W, Naylor J, Raymond C, Rosetti M, Santos R, Sheridan A, Sougnez C, et al: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.PubMedView ArticleGoogle Scholar
- Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, Holt RA, Gocayne JD, Amanatides P, Ballew RM, Huson DH, Wortman JR, Zhang Q, Kodira CD, Zheng XH, Chen L, Skupski M, Subramanian G, Thomas PD, Zhang J, Gabor Miklos GL, Nelson C, Broder S, Clark AG, Nadeau J, McKusick VA, Zinder N, et al: The sequence of the human genome. Science. 2001, 291: 1304-1351. 10.1126/science.1058040.PubMedView ArticleGoogle Scholar
- Huda A, Jordan IK: Epigenetic regulation of mammalian genomes by transposable elements. Ann N Y Acad Sci. 2009, 1178: 276-284. 10.1111/j.1749-6632.2009.05007.x.PubMedView ArticleGoogle Scholar
- López-Flores I, Garrido-Ramos MA: The repetitive DNA content of eukaryotic genomes. Genome Dyn. 2012, 7: 1-28.PubMedView ArticleGoogle Scholar
- Defraia C, Slotkin RK: Analysis of retrotransposon activity in plants. Methods Mol Biol. 2014, 1112: 195-210. 10.1007/978-1-62703-773-0_13.PubMedView ArticleGoogle Scholar
- Zedek F, Smerda J, Smarda P, Bures P: Correlated evolution of LTR retrotransposons and genome size in the genus Eleocharis. BMC Plant Biol. 2010, 10: 265-10.1186/1471-2229-10-265.PubMed CentralPubMedView ArticleGoogle Scholar
- Carle P, Saillard C, Carrere N, Carrere S, Duret S, Eveillard S, Gaurivaud P, Gourgues G, Gouzy J, Salar P, Verdin E, Breton M, Blanchard A, Laigret F, Bové JM, Renaudin J, Foissac X: Partial chromosome sequence of Spiroplasma citri reveals extensive viral invasion and important gene decay. Appl Environ Microbiol. 2010, 76: 3420-3426. 10.1128/AEM.02954-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Casjens S: Prophages and bacterial genomics: what have we learned so far?. Mol Microbiol. 2003, 49: 277-300. 10.1046/j.1365-2958.2003.03580.x.PubMedView ArticleGoogle Scholar
- Cambray G, Guerout AM, Mazel D: Integrons. Annu Rev Genet. 2010, 44: 141-166. 10.1146/annurev-genet-102209-163504.PubMedView ArticleGoogle Scholar
- Hua-Van A, Le Rouzic A, Boutin TS, Filee J, Capy P: The struggle for life of the genome’s selfish architects. Biol Direct. 2011, 6: 19-10.1186/1745-6150-6-19.PubMed CentralPubMedView ArticleGoogle Scholar
- Roberts AP, Mullany P: A modular master on the move: the Tn916 family of mobile genetic elements. Trends Microbiol. 2009, 17: 251-258. 10.1016/j.tim.2009.03.002.PubMedView ArticleGoogle Scholar
- Grindley ND, Whiteson KL, Rice PA: Mechanisms of site-specific recombination. Annu Rev Biochem. 2006, 75: 567-605. 10.1146/annurev.biochem.73.011303.073908.PubMedView ArticleGoogle Scholar
- Wicker T, Sabot F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, Flavell A, Leroy P, Morgante M, Panaud O, Paux E, SanMiguel P, Schulman AH: A unified classification system for eukaryotic transposable elements. Nat Rev Genet. 2007, 8: 973-982. 10.1038/nrg2165.PubMedView ArticleGoogle Scholar
- Kapitonov VV, Jurka J: A universal classification of eukaryotic transposable elements implemented in Repbase. Nat Rev Genet. 2008, 9: 411-412. 10.1038/nrg2165-c1. author reply 414PubMedView ArticleGoogle Scholar
- Jurka J, Kapitonov VV, Kohany O, Jurka MV: Repetitive sequences in complex genomes: structure and evolution. Annu Rev Genomics Hum Genet. 2007, 8: 241-259. 10.1146/annurev.genom.8.080706.092416.PubMedView ArticleGoogle Scholar
- Rice PA, Baker TA: Comparative architecture of transposase and integrase complexes. Nat Struct Biol. 2001, 8: 302-307. 10.1038/86166.View ArticleGoogle Scholar
- Chandler M, de la Cruz F, Dyda F, Hickman AB, Moncalian G, Ton-Hoang B: Breaking and joining single-stranded DNA: the HUH endonuclease superfamily. Nat Rev Microbiol. 2013, 11: 525-538. 10.1038/nrmicro3067.PubMedView ArticleGoogle Scholar
- Ilyina TV, Koonin EV: Conserved sequence motifs in the initiator proteins for rolling circle DNA replication encoded by diverse replicons from eubacteria, eucaryotes and archaebacteria. Nucleic Acids Res. 1992, 20: 3279-3285. 10.1093/nar/20.13.3279.PubMed CentralPubMedView ArticleGoogle Scholar
- Krupovic M: Networks of evolutionary interactions underlying the polyphyletic origin of ssDNA viruses. Curr Opin Virol. 2013, 3: 578-586. 10.1016/j.coviro.2013.06.010.PubMedView ArticleGoogle Scholar
- Goodwin TJ, Butler MI, Poulter RT: Cryptons: a group of tyrosine-recombinase-encoding DNA transposons from pathogenic fungi. Microbiology. 2003, 149: 3099-3109. 10.1099/mic.0.26529-0.PubMedView ArticleGoogle Scholar
- Goodwin TJ, Poulter RT: A new group of tyrosine recombinase-encoding retrotransposons. Mol Biol Evol. 2004, 21: 746-759. 10.1093/molbev/msh072.PubMedView ArticleGoogle Scholar
- Boocock MR, Rice PA: A proposed mechanism for IS607-family serine transposases. Mob DNA. 2013, 4: 24-10.1186/1759-8753-4-24.PubMed CentralPubMedView ArticleGoogle Scholar
- Koonin EV, Dolja VV: A virocentric perspective on the evolution of life. Curr Opin Virol. 2013, 3: 546-557. 10.1016/j.coviro.2013.06.008.PubMed CentralPubMedView ArticleGoogle Scholar
- Krupovic M, Bamford DH, Koonin EV: Conservation of major and minor jelly-roll capsid proteins in Polinton (Maverick) transposons suggests that they are bona fide viruses. Biol Direct. 2014, 9: 6-10.1186/1745-6150-9-6.PubMed CentralPubMedView ArticleGoogle Scholar
- Malik HS, Henikoff S, Eickbush TH: Poised for contagion: evolutionary origins of the infectious abilities of invertebrate retroviruses. Genome Res. 2000, 10: 1307-1318. 10.1101/gr.145000.PubMedView ArticleGoogle Scholar
- Kapitonov VV, Jurka J: Self-synthesizing DNA transposons in eukaryotes. Proc Natl Acad Sci U S A. 2006, 103: 4540-4545. 10.1073/pnas.0600833103.PubMed CentralPubMedView ArticleGoogle Scholar
- Pritham EJ, Putliwala T, Feschotte C: Mavericks, a novel class of giant transposable elements widespread in eukaryotes and related to DNA viruses. Gene. 2007, 390: 3-17. 10.1016/j.gene.2006.08.008.PubMedView ArticleGoogle Scholar
- Labrie SJ, Samson JE, Moineau S: Bacteriophage resistance mechanisms. Nat Rev Microbiol. 2010, 8: 317-327. 10.1038/nrmicro2315.PubMedView ArticleGoogle Scholar
- Makarova KS, Wolf YI, Koonin EV: Comparative genomics of defense systems in archaea and bacteria. Nucleic Acids Res. 2013, 41: 4360-4377. 10.1093/nar/gkt157.PubMed CentralPubMedView ArticleGoogle Scholar
- Westra ER, Swarts DC, Staals RH, Jore MM, Brouns SJ, van der Oost J: The CRISPRs, they are a-changin’: how prokaryotes generate adaptive immunity. Annu Rev Genet. 2012, 46: 311-339. 10.1146/annurev-genet-110711-155447.PubMedView ArticleGoogle Scholar
- Barrangou R: CRISPR-Cas systems and RNA-guided interference. Wiley Interdiscip Rev RNA. 2013, 4: 267-278. 10.1002/wrna.1159.PubMedView ArticleGoogle Scholar
- Makarova KS, Aravind L, Wolf YI, Koonin EV: Unification of Cas protein families and a simple scenario for the origin and evolution of CRISPR-Cas systems. Biol Direct. 2011, 6: 38-10.1186/1745-6150-6-38.PubMed CentralPubMedView ArticleGoogle Scholar
- Makarova KS, Wolf YI, Koonin EV: The basic building blocks and evolution of CRISPR-cas systems. Biochem Soc Trans. 2013, 41: 1392-1400.PubMedView ArticleGoogle Scholar
- Datsenko KA, Pougach K, Tikhonov A, Wanner BL, Severinov K, Semenova E: Molecular memory of prior infections activates the CRISPR/Cas adaptive bacterial immunity system. Nat Commun. 2012, 3: 945-PubMedView ArticleGoogle Scholar
- Makarova KS, Haft DH, Barrangou R, Brouns SJ, Charpentier E, Horvath P, Moineau S, Mojica FJ, Wolf YI, Yakunin AF, van der Oost J, Koonin EV: Evolution and classification of the CRISPR-Cas systems. Nat Rev Microbiol. 2011, 9: 467-477. 10.1038/nrmicro2577.PubMedView ArticleGoogle Scholar
- Yosef I, Goren MG, Qimron U: Proteins and DNA elements essential for the CRISPR adaptation process in Escherichia coli. Nucleic Acids Res. 2012, 40: 5569-5576. 10.1093/nar/gks216.PubMed CentralPubMedView ArticleGoogle Scholar
- Babu M, Beloglazova N, Flick R, Graham C, Skarina T, Nocek B, Gagarinova A, Pogoutse O, Brown G, Binkowski A, Phanse S, Joachimiak A, Koonin EV, Savchenko A, Emili A, Greenblatt J, Edwards AM, Yakunin AF: A dual function of the CRISPR-Cas system in bacterial antivirus immunity and DNA repair. Mol Microbiol. 2011, 79: 484-502. 10.1111/j.1365-2958.2010.07465.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Díez-Villaseñor C, Guzmán NM, Almendros C, García-Martínez J, Mojica FJ: CRISPR-spacer integration reporter plasmids reveal distinct genuine acquisition specificities among CRISPR-Cas I-E variants of Escherichia coli. RNA Biol. 2013, 10: 792-802. 10.4161/rna.24023.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim TY, Shin M, Huynh Thi Yen L, Kim JS: Crystal structure of Cas1 from Archaeoglobus fulgidus and characterization of its nucleolytic activity. Biochem Biophys Res Commun. 2013, 441: 720-725. 10.1016/j.bbrc.2013.10.122.PubMedView ArticleGoogle Scholar
- Wiedenheft B, Zhou K, Jinek M, Coyle SM, Ma W, Doudna JA: Structural basis for DNase activity of a conserved protein implicated in CRISPR-mediated genome defense. Structure. 2009, 17: 904-912. 10.1016/j.str.2009.03.019.PubMedView ArticleGoogle Scholar
- Mahillon J, Chandler M: Insertion sequences. Microbiol Mol Biol Rev. 1998, 62: 725-774.PubMed CentralPubMedGoogle Scholar
- Kojima KK, Jurka J: A superfamily of DNA transposons targeting multicopy small RNA genes. PLoS One. 2013, 8: e68260-10.1371/journal.pone.0068260.PubMed CentralPubMedView ArticleGoogle Scholar
- Peters JE, Craig NL: Tn7: smarter than we thought. Nat Rev Mol Cell Biol. 2001, 2: 806-814. 10.1038/35099006.PubMedView ArticleGoogle Scholar
- Parks AR, Peters JE: Tn7 elements: engendering diversity from chromosomes to episomes. Plasmid. 2009, 61: 1-14. 10.1016/j.plasmid.2008.09.008.PubMed CentralPubMedView ArticleGoogle Scholar
- Könneke M, Bernhard AE, de la Torre JR, Walker CB, Waterbury JB, Stahl DA: Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature. 2005, 437: 543-546. 10.1038/nature03911.PubMedView ArticleGoogle Scholar
- Park SJ, Kim JG, Jung MY, Kim SJ, Cha IT, Ghai R, Martin-Cuadrado AB, Rodriguez-Valera F, Rhee SK: Draft genome sequence of an ammonia-oxidizing archaeon, “Candidatus Nitrosopumilus sediminis” AR2, from Svalbard in the Arctic Circle. J Bacteriol. 2012, 194: 6948-6949. 10.1128/JB.01869-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Bath C, Cukalac T, Porter K, Dyall-Smith ML: His1 and His2 are distantly related, spindle-shaped haloviruses belonging to the novel virus group, Salterprovirus. Virology. 2006, 350: 228-239. 10.1016/j.virol.2006.02.005.PubMedView ArticleGoogle Scholar
- Krupovic M, Quemin ER, Bamford DH, Forterre P, Prangishvili D: Unification of the globally distributed spindle-shaped viruses of the Archaea. J Virol. 2014, 88: 2354-2358. 10.1128/JVI.02941-13.PubMed CentralPubMedView ArticleGoogle Scholar
- Pietilä MK, Atanasova NS, Manole V, Liljeroos L, Butcher SJ, Oksanen HM, Bamford DH: Virion architecture unifies globally distributed pleolipoviruses infecting halophilic archaea. J Virol. 2012, 86: 5067-5079. 10.1128/JVI.06915-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Gorlas A, Robert C, Gimenez G, Drancourt M, Raoult D: Complete genome sequence of Methanomassiliicoccus luminyensis, the largest genome of a human-associated Archaea species. J Bacteriol. 2012, 194: 4745-10.1128/JB.00956-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Reysenbach AL, Liu Y, Banta AB, Beveridge TJ, Kirshtein JD, Schouten S, Tivey MK, Von Damm KL, Voytek MA: A ubiquitous thermoacidophilic archaeon from deep-sea hydrothermal vents. Nature. 2006, 442: 444-447. 10.1038/nature04921.PubMedView ArticleGoogle Scholar
- Filée J, Forterre P, Sen-Lin T, Laurent J: Evolution of DNA polymerase families: evidences for multiple gene exchange between cellular and viral proteins. J Mol Evol. 2002, 54: 763-773. 10.1007/s00239-001-0078-x.PubMedView ArticleGoogle Scholar
- Tahirov TH, Makarova KS, Rogozin IB, Pavlov YI, Koonin EV: Evolution of DNA polymerases: an inactivated polymerase-exonuclease module in Pol epsilon and a chimeric origin of eukaryotic polymerases from two classes of archaeal ancestors. Biol Direct. 2009, 4: 11-10.1186/1745-6150-4-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Andrade MA, Petosa C, O’Donoghue SI, Muller CW, Bork P: Comparison of ARM and HEAT protein repeats. J Mol Biol. 2001, 309: 1-18. 10.1006/jmbi.2001.4624.PubMedView ArticleGoogle Scholar
- Iyer LM, Leipe DD, Koonin EV, Aravind L: Evolutionary history and higher order classification of AAA + ATPases. J Struct Biol. 2004, 146: 11-31. 10.1016/j.jsb.2003.10.010.PubMedView ArticleGoogle Scholar
- Bell SD, Botting CH, Wardleworth BN, Jackson SP, White MF: The interaction of Alba, a conserved archaeal chromatin protein, with Sir2 and its regulation by acetylation. Science. 2002, 296: 148-151. 10.1126/science.1070506.PubMedView ArticleGoogle Scholar
- Aravind L, Iyer LM, Leipe DD, Koonin EV: A novel family of P-loop NTPases with an unusual phyletic distribution and transmembrane segments inserted within the NTPase domain. Genome Biol. 2004, 5: R30-10.1186/gb-2004-5-5-r30.PubMed CentralPubMedView ArticleGoogle Scholar
- Filée J, Siguier P, Chandler M: Insertion sequence diversity in archaea. Microbiol Mol Biol Rev. 2007, 71: 121-157. 10.1128/MMBR.00031-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Klassen R, Meinhardt F: Linear protein-primed replicating plasmids in eukaryotic microbes. Microbiol Monogr. 2007, 7: 188-216.Google Scholar
- Salas M: Protein-priming of DNA replication. Annu Rev Biochem. 1991, 60: 39-71. 10.1146/annurev.bi.60.070191.000351.PubMedView ArticleGoogle Scholar
- Krupovic M, Koonin EV: Evolution of eukaryotic single-stranded DNA viruses of the Bidnaviridae family from genes of four other groups of widely different viruses. Sci Rep. 2014, in pressGoogle Scholar
- Nunez JK, Kranzusch PJ, Noeske J, Wright AV, Davies CW, Doudna JA: Cas1-Cas2 complex formation mediates spacer acquisition during CRISPR-Cas adaptive immunity. Nat Struct Mol Biol. 2014, doi:10.1038/nsmb.2820Google Scholar
- Chylinski K, Makarova KS, Charpentier E, Koonin EV: Classification and evolution of the protein and RNA components of type II CRISPR-Cas systems. Nucleic Acids Res. in pressGoogle Scholar
- Kapitonov VV, Jurka J: RAG1 core and V(D)J recombination signal sequences were derived from Transib transposons. PLoS Biol. 2005, 3: e181-10.1371/journal.pbio.0030181.PubMed CentralPubMedView ArticleGoogle Scholar
- Forterre P, Prangishvili D: The major role of viruses in cellular evolution: facts and hypotheses. Curr Opin Virol. 2013, 3: 558-565. 10.1016/j.coviro.2013.06.013.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView ArticleGoogle Scholar
- Marchler-Bauer A, Bryant SH: CD-Search: protein domain annotations on the fly. Nucleic Acids Res. 2004, 32: W327-W331. 10.1093/nar/gkh454.PubMed CentralPubMedView ArticleGoogle Scholar
- Söding J: Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005, 21: 951-960. 10.1093/bioinformatics/bti125.PubMedView ArticleGoogle Scholar
- Okonechnikov K, Golosova O, Fursov M: Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012, 28: 1166-1167. 10.1093/bioinformatics/bts091.PubMedView ArticleGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMed CentralPubMedView ArticleGoogle Scholar
- Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M: ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006, 34: D32-D36. 10.1093/nar/gkj014.PubMed CentralPubMedView ArticleGoogle Scholar
- Pei J, Kim BH, Grishin NV: PROMALS3D: a tool for multiple protein sequence and structure alignments. Nucleic Acids Res. 2008, 36: 2295-2300. 10.1093/nar/gkn072.PubMed CentralPubMedView ArticleGoogle Scholar
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 113-10.1186/1471-2105-5-113.PubMed CentralPubMedView ArticleGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.PubMedView ArticleGoogle Scholar
- Yutin N, Makarova KS, Mekhedov SL, Wolf YI, Koonin EV: The deep archaeal roots of eukaryotes. Mol Biol Evol. 2008, 25: 1619-1630. 10.1093/molbev/msn108.PubMed CentralPubMedView ArticleGoogle Scholar
- Price MN, Dehal PS, Arkin AP: FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009, 26: 1641-1650. 10.1093/molbev/msp077.PubMed CentralPubMedView ArticleGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013, 30: 2725-2729. 10.1093/molbev/mst197.PubMed CentralPubMedView 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/4.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.