Skip to main content

Features of a novel protein, rusticalin, from the ascidian Styela rustica reveal ancestral horizontal gene transfer event



The transfer of genetic material from non-parent organisms is called horizontal gene transfer (HGT). One of the most conclusive cases of HGT in metazoans was previously described for the cellulose synthase gene in ascidians.


In this study we identified a new protein, rusticalin, from the ascidian Styela rustica and presented evidence for its likely origin by HGT. Discernible homologues of rusticalin were found in placozoans, coral, and basal Chordates. Rusticalin was predicted to consist of two distinct regions, an N-terminal domain and a C-terminal domain. The N-terminal domain comprises two cysteine-rich repeats and shows remote similarity to the tick carboxypeptidase inhibitor. The C-terminal domain shares significant sequence similarity with bacterial MD peptidases and bacteriophage A500 L-alanyl-D-glutamate peptidase. A possible transfer of the C-terminal domain by bacteriophage was confirmed by an analysis of noncoding sequences of C. intestinalis rusticalin-like gene, which was found to contain a sequence similar to the bacteriophage A500 recombination site. Moreover, a sequence similar to the bacteriophage recombination site was found to be adjacent to the cellulose synthase catalytic subunit gene in the genome of Streptomices sp., the donor of ascidian cellulose synthase.


The C-terminal domain of rusticalin and rusticalin-like proteins is likely to be horizontally transferred by the bacteriophage A500. A common mechanism involving bacteriophage mediated gene transfer can be proposed for at least two HGT events in ascidians.


Ascidians are marine benthic animals from the subphylum Tunicata (Urochordata), which is considered the closest living sister group to vertebrates based on genome analysis [1]. The name Tunicata derives from the unique exoskeleton of these animals, the tunic, comprising both proteins and carbohydrates [2]. A remarkable feature of tunicates is biosynthesis and incorporation of cellulose into their tunic. The ascidian life cycle includes a mobile larva possessing a notochord and a sessile filter-feeding adult stage [3]. Ascidians harbor diverse microbiota [4], and their cellulose synthase is thought to have been acquired by horizontal gene transfer (HGT) from the bacterial Streptomyces sp. genome [5, 6]. The adaptive importance of HGT is supported by studies showing that mutants of cellulose synthase exhibit defects in metamorphosis and maintaining a sessile lifestyle, suggesting that it was an acquisition of cellulose synthesizing ability that permitted ascidians to evolve their sessile lifestyle [7].

Most of the described cases of HGT between prokaryotes and eukaryotes are thought to have involved transfer of genes from former to the latter [8, 9]. Possessors of former prokaryotic genes include multicellular animals [10] and, in particular, chordates [11, 12]. The fraction of horizontally acquired genes in a eukaryotic genome can reach 8%, as was described for the bdelloid rotifer Adineta vaga [13]. It has been shown that some of these horizontally transferred genes are expressed and produce functional protein products [14, 15]. Possible mechanisms of HGT between prokaryotes and eukaryotes are widely discussed, with viruses being considered as the most probable vectors of transmission into the genome [16, 17]. The existence of nuclear localization signals in bacteriophage proteins covalently bound to viral DNA lends support to this hypothesis. Facilitation of gene delivery into the eukaryotic nucleus by these signal sequences has been confirmed experimentally [18]. A broad range of gene engineering techniques adopting virus vectors for eukaryotic cells transformation in vitro and in vivo [19, 20] may provide further evidence in support of this hypothesis.

Compelling evidence supports the HGT of the cellulose synthase gene of the ascidian Ciona intestinalis [5]. This gene is expressed in the tunic-producing epidermis [7, 21]. Apart from epidermal cell layer tunic formation involves also blood cells [22, 23]. Several morphotypes of blood cells have been described for ascidians [22, 24,25,26,27], including hyalinocytes. In the blood of a solitary ascidian Styela rustica hyalinocytes and morula cells are two dominating cell groups, with an average abundance of 38 and 56%, respectively [22]. Hyalinocytes are characterized by the presence of numerous small granules. Their density is low, and so they can be separated by density gradient centrifugation [23, 28].

In this work we describe a novel protein, rusticalin, isolated from hyalinocytes of S. rustica and discuss its possible origin by HGT.


cDNA cloning and sequence analysis

Whole blood cells were separated by discontinuous percoll gradient and analyzed by SDS-PAGE (Fig. 1). The upper fraction above 35% percoll containing mainly hyalinocytes showed a major protein band of 23 kDa on SDS-PAGE. This band was subjected to trypsin digestion and MS/MS de novo sequencing, yielding a 7-residue-long peptide GNSYIRC. As a first attempt to find homologous proteins in databases the peptide was queried by tBLASTn against EST DataBase limited to Tunicata or without limitations, but showed a lack of reliable similarity. Therefore, the sequence information was used to design degenerate primers and to amplify full-length rusticalin cDNA through 3′ and 5′ Rapid Amplification of cDNA Ends (RACE) PCR. The rusticalin cDNA was 1002 bp, comprising a 5′-untraslated region of 111 bp, an open reading frame of 690 bp and a 3′-untraslated region of 201 bp. The first ATG at position 94–96 was assigned as the start codon. Two polyadenilation signals (AATAAA) were found 23 and 112 bp upstream of the poly(A) + tail. The ORF encoded a protein of 230 amino acid residues, including predicted signal peptide of 18 amino acid residues (Fig. 2, GenBank accession number MH115429).

Fig. 1

SDS-PAGE of blood cells proteins of ascidian Styela rustica. Lane T: whole blood cells; lane 1: the upper fraction of blood cells from discontinuous percoll gradient (above 35% percoll); lane 2: the intermediate fraction (between 35 and 45% percoll); lane 3: morula cells fraction from discontinuous percoll gradient (between 45 and 60% percoll); lane 4: whole blood cells proteins transferred to PVDF membrane and stained with Ponceau S. The major band of the upper (hyalinocytes) fraction of 23 kDa (black arrow) subjected to further matrix-assisted laser desorption/ionization tandem mass spectrometry MALDI MS/MS analysis. M – Molecular weight markers in kDa

Fig. 2

Deduced amino acid sequence of rusticalin. The predicted signal peptide is printed in red, the internal repeats detected with REPRO and RADAR are underlined, with cysteines shaded yellow, and a stretch of residues with relative solvent accessibility > 50% shaded black. Peptide sequence detected by de novo sequencing is printed in blue

The mature protein comprises 212 amino acid residues with a theoretical molecular mass of 23,309.7 Da. It contains 11 negatively charged residues (Asp + Glu) and 29 positively charged residues (Arg + Lys) yielding a calculated pI of 9.33. N-terminal region of rusticalin was found to contain two repeats, 34 and 33 residues in length, with a common cysteine spacing motif Cx6Cx6-7Cx8Cx7CC (Fig. 2). These protein regions are further referred to as cysteine-rich repeats.

Computational tools were applied to detect putative domains and to predict relative solvent accessibility and disorder in the sequence. Scooby-domain tool predicted two domains with a boundary at Ser95. Prediction of the relative solvent accessibility revealed a stretch of eight highly exposed residues at positions Ser94-Ser102. Protein backbone disorder prediction showed the existence of two rigid regions linked by a short flexible region (Fig. 3). These data suggest that the hydrophilic and flexible region identified may be involved in the formation of a linker about six amino acid residues long connecting N- and C-terminal domains.

Fig. 3

Prediction of disorder. Disordered regions predicted using Disopred3 (red line), SPOT (blue line) and protein backbone dynamics predicted using DynaMine (green line) for Styela rustica rusticalin

Localization of rusticalin mRNA

The localization of rusticalin mRNA in the blood cells was examined with fluorescent in situ hybridization (FISH). Confocal microscopy showed that flattened cells containing numerous small spherical granules were labeled (Fig. 4II, III). These cells were clearly identified as hyalinocytes based on the presence of characteristic granules (Fig. 4Ia), which were absent in morula cells (Fig. 4Ib). Thus the hybridization signal was restricted to the cytoplasmic space of hyalinocytes (Probe 2: Fig. 4IIb). The same results were observed with Probe 1; the negative control showed no hybridization signal (data not shown).

Fig. 4

FISH detection of rusticalin mRNA in the blood cells. I Hyalinocytes (H) and morula cells (Mc), hematoxylin and eosin staining, DIC microscopy. Note hyalinocytes containing numerous small spherical granules (arrows). Scale bar, 5 μm. II Confocal sections showing the distribution of transcripts (b) within the cells with the morphology of hyalinocytes (a). Scale bar, 10 μm. The hybrids with Probe 2 were detected with streptavidin-Alexa594 (red pseudo color). 4′-6-Diamidino-2-phenylindole (DAPI) was used as a general DNA dye (blue pseudo color)

Similarity search

The workflow of rusticalin sequence analysis is shown in Fig. 5. Searches against transcriptome and non-redundant protein databases using tBLASTn for finding close homologues and HHblits for remote similarity search identified discernible rusticalin-like proteins only in tunicates (Oikopleura dioica, Ciona intestinalis, Ciona savignyi, Diplosoma listerianum, and Botryllus schlosseri), cephalochordates (Branchiostoma floridae and Asymmetron lucayanum) and basal multicellular animals (coral Alveopora japonica and placozoan Trichoplax adhaerens) (Table 1). Multiple sequence alignment of all newly identified rusticalin-like proteins was queried against UniProtKB and the NCBI non-redundant protein databases. No other significant hits containing both predicted structural domains and covering more than 90% of a query were found. These results indicate that rusticalin-like proteins are taxonomically restricted to placozoans, corals, and basal chordates.

Fig. 5

Workflow of rusticalin sequence analysis

Table 1 Results of iterative search strategy used to identify rusticalin-like proteins

The alignment of all the proteins showed high sequence similarity and a highly conserved cysteine spacing at the cysteine-rich repeats (Fig. 6). All proteins were found to contain N-terminal signal peptides, and were predicted to be secretory. Notably, rusticalin lacks 40 C-terminal residues, in contrast to its identified homologs.

Fig. 6

Multiple sequence alignment of rusticalin-like proteins. Predicted signal peptides were trimmed out before alignment. Note the perfect match and conservation of cysteine residues highlighted in yellow. N- and C-terminal domains of rusticalin-like proteins are shown in black. Note the lack of 40 C-terminal residues in rusticalin. The tentative linker region between N- and C-terminal domains is shown in gray. Two cysteine-rich repeats inside the N-terminal domain are shown in blue. Amino acid residues above 85% identity threshold are colored according to their physicochemical properties. Styru – Styela rustica, Bosch – Botryllus schlosseri, Cioin – Ciona intestinalis, Ciosa – Ciona savignyi, Dipli – Diplosoma listerianum, Oikdi – Oikopleura dioica, Brafl – Branchiostoma floridae, Alvja – Alveopora japonica, Triad – Trichoplax adhaerens

In order to characterize predicted domains of rusticalin-like proteins we performed a remote similarity search (HHpred) using a multiple sequence alignment separately for each domain as query. The multiple sequence alignment generated for individual cysteine-rich repeats (Fig. 6,) and searched with HHpred against Pfam and SCOP databases revealed similarity with β-defensin family and β-defensin-like fold (SCOP g.9.1), respectively (Fig. 7a). Additionally, the multiple sequence alignment of the N-terminal domains containing a pair of cysteine-rich repeats, searched against PDB database, showed similarity to tick carboxypeptidase inhibitor (PDB ID 1ZLH) (Fig. 7a). Remarkably, tick carboxypeptidase inhibitor is structurally related to β-defensin-like fold and is the only described protein structure comprising two β-defensin repeats. Thus, in silico analysis suggests that the N-terminal domain of rusticalin-like proteins may have a tertiary structure similar to the tick carboxypeptidase inhibitor, acting as a double-headed enzyme inhibitor [29].

Fig. 7

Sequence-structure alignment of N- and C-terminal domains of Botryllus schlosseri rusticalin sequence. a Alignment of the N-terminal domain with β-defensin-like fold of tick carboxypeptidase inhibitor (PDB ID 1ZLH). Conserved cysteine residues are highlighted in black. The sequence identity is 26%. b Alignment of C-terminal domain with Hedgehog/DD-peptidase fold of bacteriophage A500 L-alanyl-D-glutamate peptidase (PDB ID 2VO9). Amino acid residues involved in Zn2+ binding by 2VO9 are highlighted in red, catalytic and substrate-binding residues are highlighted in green and yellow, respectively. The sequence identity is 21%. ss_dssp – template secondary structure as determined by DSSP. The secondary structure is labeled ‘H’ for α-helix, ‘E’ for β-strand, and ‘C’ for coil

In order to determine the nature of the C-terminal domain we queried its multiple sequence alignment against Pfam, SCOP, and PDB databases. The search in Pfam database showed that the C-terminal domain of rusticalin-like proteins share a significant sequence similarity with Peptidase_MD clan (Pfam ID: CL0170). Most of proteins belonging to that clan are bacterial cell-wall degradation enzymes suggesting that the C-terminal domain might originate from a bacterial genome. The search in SCOP database revealed that the C-terminal domain matched structurally with Hedgehog/DD-peptidase fold (SCOP d.65.1). Catalytic, substrate binding, and Zn-binding residues of MD peptidases were conserved (Fig. 7b), thus rusticalin-like proteins are likely to have peptidase activity. However, rusticalin itself appears to lack this activity due to the absence of 40 C-terminal residues. Finally, a high sequence similarity (Fig. 7b, 21%, E-value of 1.5E-16) of C-terminal domain with bacteriophage A500 L-alanyl-D-glutamate peptidase (PDB ID 2VO9,) indicates a possible role of bacteriophage in horizontal transfer of the C-terminal domain coding sequence from the bacterial genome.

Evidence of a horizontal gene transfer (HGT) event

Bacteriophage A500 site-specific recombination involves the 3′ region of the bacterial tRNA gene [30]. Thus all genomes containing rusticalin-like proteins were searched for tRNA genes neighboring rusticalin-like genes. The rusticalin-like gene of C. intestinalis (Gene ID: 100185212) contains seven tRNA genes at antisense orientation situated upstream of the gene and inside the second and third introns from 5′-end (Fig. 8a). Multiple sequence alignment shows that seven tRNA genes of C. intestinalis are highly similar, with the sequence identity from 95 to 100% (Fig. 9). The third intron containing tRNA genes is adjacent to the C-terminal domain. Alignment of tRNA gene (Gene ID: 108950108) lying inside the third intron with bacteriophage A500 recombination site (AttP) showed the presence of a similar sequence (Fig. 8b). Thus in C. intestinalis rusticalin-like protein L-alanyl-D-glutamate peptidase domain is adjacent to the intron containing a region resembling the bacteriophage recombination site, confirming the domain’s horizontal transfer by means of a viral genome. We also conducted a nucleotide BLAST of the bacteriophage AttP site against all Tunicata genomic sequences, which gave a hit with B. schlosseri contig89252 (Fig. 10). We can conclude that a sequence similar to bacteriophage A500 AttP is present in Tunicata genomes. Nucleotide BLAST against T. adhaerensis genome gave no positive results.

Fig. 8

Rusticalin-like gene from Ciona intestinalis contains bacteriophage A500 recombination site inside the non-coding region. a Position and antisense orientation of seven tRNA-Arg genes inside the non-coding regions of C. intestinalis rusticalin-like gene (Gene ID: 100185212). The first tRNA gene is situated upstream of the protein coding sequence, four of tRNA genes are inside the second intron and two are inside the third intron neighboring the C-terminal domain coding region. b Alignment of bacteriophage A500 recombination site (AttP) with Ciona intestinalis tRNA gene (Gene ID: 108950122) situated inside the third intron. The sequence identity and score are, respectively, 65.8% and 73 (calculated by EMBOSS Matcher)

Fig. 9

Alignment of seven Ciona intestinalis tRNA genes neighboring rusticalin-like gene coding sequence. Asterisks indicate conservative positions. Three of the tRNA genes differ from the other tRNA genes at four nucleotide positions (shaded grey). The sequence identity is 95–100%

Fig. 10

Result of nucleotide BLAST of bacteriophage A500 recombination site (AttP) against Tunicata genome sequences. Alignment of AttP site with Botryllus schlosseri genome sequence - contig89252. The sequence identity is 88%, E-value: 0.17

We analyzed the genome of Streptomices sp., the prokaryote donor of ascidian cellulose synthase gene [5]. Cellulose synthase catalytic subunit gene (bcsA) was found to be adjacent to tRNA-Lys gene in this genome (Fig. 11a). Pairwise alignment of the tRNA gene with bacteriophage A500 AttP showed the presence of a highly similar sequence (Fig. 11b). This result suggests the involvement of tRNA gene in HGT of cellulose synthase into the tunicates genome.

Fig. 11

Donor of ascidian cellulose synthase Streptomyces sp. contains bacteriophage recombination site adjacent to cellulose synthase gene. a Streptomyces sp. cellulose synthase catalytic subunit gene (bcsA) is adjacent to tRNA-lys gene. b Alignment of bacteriophage A500 recombination site (AttP) with the tRNA-Lys gene (APS67_000733). The sequence identity and score are, respectively, 91.2% and 143 (calculated by EMBOSS Matcher)

Fig. 12

Alignment of C-terminal domains of rusticalin-like proteins with designated positions of introns. Blue vertical bars represent the positions of introns in the corresponding DNA sequence of: Cioin_3 – Ciona intestinalis XP_002122335.1; Cioin_1 – Ciona intestinalis XP_002128942.1; Brafl – Branchiostoma floridae XP_002588042.1; Triad – Trichoplax adhaerens XP_002117795.1

Taking into account that bacteriophage A500 genes do not contain introns [30] we used the presence of introns and their positions to predict the number of independent HGT events. Information about the positions of introns was available for C. intestinalis, B. floridae, and T. adhaerens rusticalin-like genes. The positions were mapped on the corresponding protein sequences. Two introns were found to be located inside the C-terminal domain coding region (Fig. 12), and their positions were strictly conserved in the sequences analyzed. This fact suggests that the C-terminal domain was formed as a result of a single gene transfer event of L-alanyl-D-glutamate peptidase. Synonymous distances counted between bacteriophage A500 enzyme and C-terminal domain of those four proteins indicated that the shortest distance of 34 substitutions is in bacteriophage A500 and C. intestinalis (Cioin_1) comparison. Based on this data we speculate that the first acceptor of a foreign gene belonged to the Tunicata lineage.


Specific expression of rusticalin in hyalinocytes

As previously shown, percoll gradients are suitable for isolation of cell populations in marine invertebrates. They have been successfully used for identification of cell-type-specific proteins through antibody (AB) production [23, 31] or by MALDI MS\MS analysis with subsequent RACE PCR [32, 33]. About 40% of the blood cells in the ascidian Styela rustica are represented by hyalinocytes [22]. Hyalinocytes or their equivalents in other ascidian species perform functions such as phagocytosis [22, 27], cytokine synthesis [34], and protease release upon LPS induction [28]. In order to isolate the rusticalin protein of hyalinocytes and describe its gene, we conducted MALDI and RACE. DNA-RNA FISH of the newly identified gene confirmed its specific expression in hyalinocytes (Fig. 4). The rusticalin gene with the deduced amino acid sequence was compared to other genome and transcriptome sequences from many species using both BLAST search and methods specialized for remote similarity search – Hhblits and HHpred. This approach allowed us to characterize a new protein, rusticalin, and predict properties for rusticalin as well as for group of homologous rusticalin-like Proteins. Rusticalin-like proteins are present in basal chordates and, also in primitive multicellular animals: coral A. japonica and placozoan T. adhaerens.

Putative function of rusticalin-like proteins

Prediction of protein disorder and solvent accessibility for rusticalin showed the existence of two distinct structural domains (Fig. 3). N-terminal domain contained two cysteine-rich repeats. Querying of sequence and predicted structure of cysteine-rich repeats in protein databases showed that they resembled β-defensins, antimicrobial proteins responsible for the lysis of pathogens [35,36,37] by disrupting their membranes [38]. On the other hand, the C-terminal domain of rusticalin-like proteins was identified as a Peptidase_MD clan member and, more specifically, as being close to L-alanyl-D-glutamate peptidase. Catalytic, substrate binding, and Zn-binding sites of the enzyme [39] were conserved in all rusticalin-like proteins suggesting that they may have peptidase activity. Other members of MD peptidases are bacterial cell-wall digesting enzymes [39,40,41,42]. Though the precise function of rusticalin-like proteins cannot be identified yet, we may venture a guess that the N-terminal domain perforates bacterial cell walls while the C-terminal domain digests them. Accordingly, all rusticalin-like proteins are predicted by TargetP to be secretory. The fact that rusticalin is specific to hyalinocytes does not contradict its putative immune function since at least hyaline amoebocytes are also known to be capable of phagocytosis [27]. Another protein previously characterized as Zn-dependent metallo-protease from the ascidian Halocynthia roretzi hemocytes is activated by lipopolysaccharide (LPS) [28, 43, 44] and hence might also be involved in immune reactions.

At the same time, a pair of the cysteine-rich repeats analyzed separately showed a significant similarity with carboxypeptidase inhibitor of the tick Rhipicephalus bursa. This protein is also related to β-defensin-like fold [29] but its function is to inhibit carboxypeptidase-A/B of mammalian blood [45]. Based on this finding, we propose an alternative scenario for the interaction of the N- and C-terminal domains, where the N-terminal domain exerts no bactericidal function but acts as a regulatory subunit. This mode of interaction has been described for carboxypeptidases A/B (M14) [46], for zinc-dependent matrix metalloproteases (MMPs) [47], and also for LytM [48], which is related to MD peptidases [49]. Thus, the ancestral state of the N-terminal domain’s function might have been the perforation of the bacterial membrane. Whatever the case, putative functions of the newly described protein should be verified experimentally by production of recombinant protein. Rusticalin of S. rustica is 40 amino acids shorter and lacks a part of the active site. This means that it cannot perform an enzymatic function but might still be involved in the signaling pathways of the immune reaction [50, 51], similarly to the Hedgehog signaling molecule, another member of peptidase MD family [52].

Possible horizontal gene transfer (HGT)

Cellulose synthase of the ascidian C. intestinalis provides one of the clearest examples of HGT [5]. In the present study we described another ascidian protein, rusticalin, whose C-terminal domain probably originated by means of HGT from a bacterial cell-wall digesting enzyme. Moreover similarity with bacteriophage A500 L-alanyl-D-glutamate peptidase suggests a possible involvement of a bacteriophage as a vector. This hypothesis is supported by the fact that the C-terminal domain belongs to bacterial MD peptidases (Pfam ID CL0170) and at the same time shows significant sequence similarity with bacteriophage protein (E-value 1.5e-16). It is further confirmed by an analysis of noncoding regions of C. intestinalis rusticalin-like gene, which contained a sequence similar to the bacteriophage A500 recombination site [30]. While many cases of HGT are described based on sequence similarity alone [15, 53,54,55,56,57,58,59], in the case of rusticalin we also demonstrated strong evidence of the mechanism of transfer by identifying the recombination site.

Rusticalin-like proteins are also present in a primitive multicellular animal Trichoplax adhaerens [60, 61] and the coral Alveopora japonica. However, no remains of bacteriophage A500 recombination sites were found in the T. adhaerens or A. japonica nucleotide sequences. The signatures of the bacteriophage gene transfer might have been erased from the T. adhaerens genome as a result of intron shortening [61] (Table 2) but preserved in the C. intestinalis genome, possibly, due to the possession of functioning tRNA genes inside the introns (Fig. 7). We also found that Streptomices sp., the prokaryote donor of the ascidian cellulose synthase gene [5], contained tRNA-Lys gene and a sequence similar to the bacteriophage recombination site (AttP) adjacent to the cellulose synthase catalytic subunit gene (bcsA). This fact supports the hypothesis that viral recombination with tRNA genes was involved in HGT events and suggests a common mechanism for at least these two cases of HGT.

Table 2 Intron length in rusticalin-like genes

Since T. adhaerens, A. japonica, and Chordata are distant animal relatives [62], it can’t be ruled out that HGT events for C-terminal domains of their rusticalin-like proteins were independent. Still, the position of the fourth intron inside the C-terminal domain coding region is identical for the placozoan T. adhaerens, the ascidian C. intestinalis, and the cephalochordate B. floridae. Given that the genome of the bacteriophage A500 contains no introns [30], they must have been introduced right after the gene transfer to the eukaryote genome [63]. It seems improbable that the identical intron positions are the result of an independent intron gain. Thus, we assume that the fourth intron appeared as a result of a single event of intron insertion into the C-terminal domain coding region. This means, in turn, that the C-terminal domain of rusticalin and rusticalin-like proteins emerged as a result of a single HGT event of L-alanyl-D-glutamate peptidase, inserted by the bacteriophage into the eukaryote genome. We performed a synonymous distance analysis between the bacteriophage A500 enzyme and the C-terminal domains of four rusticalin-like proteins that possess the identical intron positions. The C. intestinalis gene (Cioin_1) appeared to have shortest synonymous distance to the bacteriophage enzyme. The same gene contains a tRNA and a sequence similar to the AttP site inside its introns. Thus, this supports the hypothesis that the first HGT event mediated by a bacteriophage happened in the Tunicata lineage.


We described a new protein, rusticalin, from the hyalinocytes of the ascidian Styela rustica and predicted its features based on the sequence analysis. Discernible homologues of rusticalin were found only in basal chordates, coral, and placozoans. Sequence similarity and the presence of a putative bacteriophage recombination site support the hypothesis of transfer of the C-terminal domain from a bacteriophage genome. A similar mechanism involving bacteriophage as a vector can be proposed for the cellulose synthase catalytic subunit gene.



Ascidians Styela rustica Linnaeus (1767) were collected off Fettakh Island near the Biological Station of the Zoological Institute of the Russian Academy of Sciences at Cape Kartesh (Kandalaksha Bay, the White Sea) in June–August of 2013–2017. The ascidians were kept in cages at a depth of 3–4 m throughout the experimental period.

Collection of hemocytes

All manipulations with ascidians were carried out in a temperature-controlled room at 10 °C. Before bleeding, the animal was washed with sea water and dried with absorbent paper. Then the sampling area was sterilized with 70% ethanol and the ascidian body wall was cut with a razor blade to the muscular layer without injuring the internal organs. Hemolymph was collected from the cut with a micropipette and transferred into a tube containing an anticoagulant solution (AS) (0.3 M NaCl, 20 mM KCl, 15 mM EDTA, 10 mM HEPES pH 7.6) [23].

Discontinuous percoll gradient for hemocytes fractionation of hemocytes

Percoll solution (Sigma) was mixed with appropriate volumes of AS to obtain final concentrations of 60, 45, and 35%. Three milliliters of each mixture was overlaid sequentially into a glass centrifuge tube. The blood sample was made by pooling blood from four animals and mixing it with AS (1:1). Three milliliters of the blood sample was layered onto the percoll gradient and the tube was centrifuged in a swing rotor at 800 g for 30 min. Cells from the density boundary were collected by gentle aspiration and washed thrice in AS. The cell composition of fractions was determined by phase-contrast microscopy. The protein composition of the fractions was analyzed by SDS-PAGE.


Protein samples for SDS-PAGE were prepared out of whole blood cells or cell fractions after separation in percoll gradient. Cells were centrifuged at 800 g for 10 min, resuspended in 7 mM EDTA, 1 mM PMSF, 10% β-mercaptoetanol, and frozen (− 20 °C). After thawing the suspension was mixed with 2x loading buffer (0.3 Tris-HCl pH 6.8; 20% glycerol; 4% SDS; 5% β-mercaptoetanol) and boiled for 5 min. SDS-PAGE was performed on 15% gels with Mini-Protean II electrophoretic cell (Bio-Rad). Unstained Protein MW marker (Thermo Scientific) was used as a size standard. To visualize proteins, the gel slabs were stained with Coomassie BB R-250 (Biolot, Russia).

Protein sequencing and tandem mass spectrometry

After SDS-PAGE of whole blood cells proteins were transferred to PVDF membrane and stained with Ponceau S (Fig. 1, line 4). A protein band of apparent molecular mass 23 kDa was excised and subjected to Edman degradation (Alta Bioscience, interior code of sample: S6269, Birmingham, UK). This method provided no accurate amino acid sequence. Therefore, an equal protein band was excised from polyacrylamide gel and subjected to digestion with Proteomics Grade Trypsin (Sigma). Tryptic fragments were further extracted from the gel matrix and analyzed by MALDI MS/MS at PostGenome analysis center (http://xn--h1aaoah.xn--p1ai/services-and-rates/mass-spectrometry.html, Moscow). The resulting partial amino acid sequence was used to create nested degenerate oligonucleotide primers designed with iCODEHOP [64].

Cloning and sequencing of rusticalin cDNA

Total RNA was extracted from blood cells of S. rustica using TRI Reagent (Sigma) and reverse-transcribed with MINT cDNA synthesis kit (Evrogen) according to the manufacturer’s instructions. MINT RACE cDNA Amplification Set (Evrogen) was used for 3′ and 5’ RACE. For 3’RACE, nested degenerate oligonucleotide primers were designed using the iCODEHOP algorithm [64] on the basis of the determined amino acid sequence (Table 3; #1, 2). Primers for 5′RACE (Table 3; #3, 4) were based on the DNA sequence obtained in 3’RACE. Both 3′ and 5’PCR products were cloned in pAL2-T vector, using Quick-TA kit (Evrogen, Russia), and Sanger sequenced in Evrogen.

Table 3 Oligonucleotide primers used in the study

DNA-RNA fish

Two synthetic 26–27-mer 5′-end biotin-labeled DNA probes were used for DNA-RNA FISH. Probe 1 (/Biotin/CAGTTGTTGCTCATAACCGGCGATGC-3′) was complementary to 113–138 nucleotide region corresponding to N-terminal domain of rusticalin, while Probe 2 (/Biotin/GGCGACTCGAATTACCTTGCCCTGATA-3′) was complementary to 400–426 nucleotide region corresponding to the C-terminal domain of rusticalin. Hybridization without probe served as negative control.

Ascidian blood was collected as described above. Blood drops were transferred from the cut in the body wall directly onto a glass slide (Superfrost Plus, Menzel) and left for 20 min at 10 °C for cell attachment. The cells were fixed with 4% PFA in AS for 10 min at 10 °C and washed successively in AS, distilled water, and methanol. The slides were dried and stored frozen (− 20 °C) until use. For morphological control several slides with spread cells were resolved and stained with hematoxylin and eosin, dehydrated, and embedded in Dammar resin. Images were taken on Leica DM6000 with DIC (Nomarsky optics).

Before FISH the excessive PFA was washed off with PBT (1 × PBS, 0.1% Tween 20). Cells were pretreated with 2 μg·ml− 1 proteinase K (Thermo Scientific), 0.1% SDS in PBS for 2 min. The proteinase K was then inactivated by incubation with 200 μM PMSF. Cells were postfixed in 4% PFA and washed again with 200 μM PMSF. Excessive PFA was washed off with PBT. Endogenous biotin was blocked as described by Miller and Kubier [65]. The cells were then washed thrice for 10 min with PBS and postfixed in 4% PFA. Excessive PFA was washed with PBT.

To perform DNA-RNA FISH the cells were rinsed in 4 × SSC and prehybridized in hybridization buffer (1% dextran sulfate, 50% formamide, 1 mg·ml− 1 salmon sperm DNA in 4 × SSC) for 15 min at 36 °C. Hybridization was performed with 0.5 μM of probe in hybridization buffer for 17 h at 36 °C. After hybridization the samples were washed in 50% formamide, 4 × SSC at 36 °C and then in 0.2 × SSC, 0.1% Tween 20 at 45 °C. After blocking in 1× In Situ Hybridization Blocking solution (Vector laboratories) in PBT at 37 °C for 60 min, the probe was detected using strepavidin-Alexa594 (1:500, Life technologies) at 37 °C for 120 min. The samples were washed thrice at 37 °C in PBT, counterstained with 3 μg/ml DAPI and mounted in 80% glycerol in 1 × PBS. Fluorescent images were taken with the use of confocal laser microscope LEICA TCS SP5 MP.

Sequence analysis and database searches

The workflow of sequence analysis and database searches is shown in Fig. 5. The average molecular mass and isoelectric point of rusticalins were calculated with ProtParam [66] on the ExPASy server ( Signal peptides were predicted with Phobius at EMBL-EBI [67] and SignalP [68]. Subcellular location was predicted with SCL-Epred [69] and TargetP [70]. Globular domains were predicted with Scooby-domain. Internal repeats were identified with REPRO [71] and RADAR [72] algorithms. Relative solvent accessibility was predicted with PaleAle [73]. Disordered regions were predicted with Disopred3 [74] and SPOT-disorder [75], and protein backbone dynamics was predicted with DynaMine [76]. All secondary structure predictions were made after removal of the signal peptide.

The initial tBLASTn searches were performed against transcriptome database (EST) available at NCBI server. HHblits [77] was used to search in UniProtKB and the NCBI non-redundant protein databases. Obtained hits showing both conservation of cysteine residues and more than 90% sequence coverage were trimmed to remove putative signal peptide and aligned using MSAProbs [78]. The aligned sequences were filtered to 90% identity and subjected to remote similarity searches using HHpred [79] in PDB, SCOP, and Pfam 30.0 protein databases. Multiple sequence alignment was visualized with CHROMA software [80].

Genomic sequences and gene structure

tRNA genes positions in genomic sequences were retrieved from whole-genome shotgun sequences of Ciona intestinalis (GCA_000224145.2) and Streptomices sp. AVP053U2 isolated from Styela clava (LMTQ02000003.1) [81]. Sequences of seven C. intestinalis tRNA genes (Gene ID: 108950112, 108,950,111, 108,950,110, 108,950,109, 108,950,108, 108,950,122, 108,950,121) were obtained from NC_020179.2 genome region (Chromosome 14). Sequences of Streptomices sp. tRNA gene (APS67_000733) were obtained from region 156,485–156,560 of contig000003. tRNA genes were aligned using the Clustal Omega multiple sequence alignment program [82]. Pairwise alignment of bacteriophage recombination site sequence and tRNA genes was made in EMBOSS Matcher [83]. Database searches restricted to Tunicata were performed using BLASTn against GenBank nucleotide collections: nr/nt database, expressed sequence tags (EST), and whole-genome shotgun contigs (WGS).

Information about gene structure was available in GenBank for four rusticalin-like sequences: two Ciona intestinalis genes GeneID:100181995, XM_002122299.4, XP_002122335.1 and GeneID:100185212, XM_002128906.4, XP_002128942.1; Branchiostoma floridae gene GeneID:7231622, XM_002587996.1, XP_002588042.1 and Trichoplax adhaerens gene GeneID:6759007, XM_002117759.1, XP_002117795.1. Intron positions were mapped on the corresponding amino acid sequences preserving alignment.

The same gene sequences with addition of bacteriophage A500 gene (GeneID:5601386) were used to calculate synonymous distances with SNAP v2.1.1 [84]. Distances were calculated based on codon-alignment preserving alignment of amino acid sequences.



Antibodies, LPS – lipopolysaccharide


Anticoagulant solution




Fluorescent in situ hybridization


Horizontal gene transfer


Matrix-assisted laser desorption/ionization tandem mass spectrometry


Matrix metalloprotease


1×PBS, 0.1% Tween 20


Rapid Amplification of cDNA Ends


  1. 1.

    Delsuc F, Brinkmann H, Chourrout D, Philippe H. Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature. 2006;439:965.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Daele Y, Revol J, Gaill F, Goffinet G. Characterization and supramolecular architecture of the cellulose-protein fibrils in the tunic of the sea peach (Halocynthia papillosa, Ascidiacea, Urochordata). Biol Cell. 1992;76:87–96.

    Article  Google Scholar 

  3. 3.

    Shenkar N, Swalla BJ. Global diversity of ascidiacea. PLoS One. 2011;6(6):e20657.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Schreiber L, Kjeldsen KU, Funch P, Jensen J, Obst M, López-Legentil S, et al. Endozoicomonas are specific, facultative symbionts of sea squirts. Front Microbiol. 2016;7:1042.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Nakashima K, Yamada L, Satou Y, Azuma J, Satoh N. The evolutionary origin of animal cellulose synthase. Dev Genes Evol. 2004;214(2):81–8.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Sagane Y, Zech K, Bouquet J-M, Schmid M, Bal U, Thompson EM. Functional specialization of cellulose synthase genes of prokaryotic origin in chordate larvaceans. Development. 2010;137(9):1483 LP–1492.

    Article  CAS  Google Scholar 

  7. 7.

    Sasakura Y, Ogura Y, Treen N, Yokomori R, Park S-J, Nakai K, et al. Transcriptional regulation of a horizontally transferred gene from bacterium to chordate. Proc R Soc B Biol Sci. 2016;283(1845).

  8. 8.

    Andersson JO. Gene transfer and diversification of microbial eukaryotes. Annu Rev Microbiol. 2009;63(1):177–93.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Tucker RP. Horizontal gene transfer in choanoflagellates. J Exp Zool Part B Mol Dev Evol. 2012;320(1):1–9.

    Article  CAS  Google Scholar 

  10. 10.

    Boto L. Horizontal gene transfer in the acquisition of novel traits by metazoans. Proc R Soc B Biol Sci. 2014;281(1777):20132450.

    Article  Google Scholar 

  11. 11.

    Graham LA, Lougheed SC, Ewart KV, Davies PL. Lateral transfer of a lectin-like antifreeze protein gene in fishes. PLoS One. 2008;3(7):e2616.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  12. 12.

    Riley DR, Sieber KB, Robinson KM, White JR, Ganesan A, Nourbakhsh S, et al. Bacteria-human somatic cell lateral gene transfer is enriched in cancer samples. Eisen JA, editor. PLoS Comput Biol. 2013;9(6):e1003107.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Flot J-F, Hespeels B, Li X, Noel B, Arkhipova I, Danchin EGJ, et al. Genomic evidence for ameiotic evolution in the bdelloid rotifer Adineta vaga. Nature. 2013;500:453.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  14. 14.

    Boschetti C, Carr A, Crisp A, Eyres I, Wang-Koh Y, Lubzens E, et al. Biochemical diversification through foreign gene expression in bdelloid rotifers. PLoS Genet. 2012;8(11):e1003035.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Ren Q, Wang C, Jin M, Lan J, Ye T, Hui K, et al. Co-option of bacteriophage lysozyme genes by bivalve genomes. Open Biol. 2017;7(1):160285.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  16. 16.

    Ryan F. The mysterious world of the human genome. Amherst: Prometheus Books; 2016. p. 300.

    Google Scholar 

  17. 17.

    Gilbert C, Peccoud J, Chateigner A, Moumen B, Cordaux R, Herniou EA. Continuous influx of genetic material from host to virus populations. PLoS Genet. 2016;12(2):e1005838.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  18. 18.

    Redrejo-Rodríguez M, Muñoz-Espín D, Holguera I, Mencía M, Salas M. Functional eukaryotic nuclear localization signals are widespread in terminal proteins of bacteriophages. Proc Natl Acad Sci. 2012;109(45):18482 LP–18487.

    Article  Google Scholar 

  19. 19.

    Chira S, Jackson CS, Oprea I, Ozturk F, Pepper MS, Diaconu I, et al. Progresses towards safe and efficient gene therapy vectors. Oncotarget. 2015;6(31):30675–703.

    PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Chatterjee S, Sullivan HA, MacLennan BJ, Xu R, Hou Y, Lavin TK, et al. Nontoxic, double-deletion-mutant rabies viral vectors for retrograde targeting of projection neurons. Nat Neurosci. 2018;21(4):638–46.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Matthysse AG, Deschet K, Williams M, Marry M, White AR, Smith WC. A functional cellulose synthase from ascidian epidermis. Proc Natl Acad Sci U S A 2004;101(4):986 LP-991.

  22. 22.

    Chaga OY. Blood cells in the ascidian Styela (Goniocarpa) rustica. I. Histological analysis. Tsitol. 1998;40:31–44.

    Google Scholar 

  23. 23.

    Podgornaya OI, Shaposhnikova TG. Antibodies with the cell-type specificity to the morula cells of the solitary ascidians Styela rustica and Bolteni echinata. Cell Struct Funct. 1998;23(6):349–55.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  24. 24.

    Radford JL, Hutchinson AE, Burandt M, Raftos DA. A Hemocyte classification scheme for the tunicate Styela plicata. Acta Zool. 1998;79(4):344–50.

    Article  Google Scholar 

  25. 25.

    Hirose E, Shirae M, Saito Y. Ultrastructures and classification of circulating hemocytes in 9 botryllid ascidians (chordata: ascidiacea). Zool Sci. 2003;20(5):647–56.

    PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Ballarin L, Kawamura K. The hemocytes of Polyandrocarpa mysakiensis : morphology and immune-related activities. ISJ. 2009;6:154–61.

    Google Scholar 

  27. 27.

    Cima F, Peronato A, Ballarin L. The haemocytes of the colonial aplousobranch ascidian Diplosoma listerianum: structural, cytochemical and functional analyses. Micron. 2017;102:51–64.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Azumi K, Satoh N, Yokosawa H. Functional and structural characterization of hemocytes of the solitary ascidian, Halocynthia roretzi. J Exp Zool. 1993;265:309–16.

    CAS  Article  Google Scholar 

  29. 29.

    Arolas JL, Popowicz GM, Lorenzo J, Sommerhoff CP, Huber R, Aviles FX, et al. The three-dimensional structures of tick carboxypeptidase inhibitor in complex with A/B carboxypeptidases reveal a novel double-headed binding mode. J Mol Biol. 2005;350(3):489–98.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Dorscht J, Klumpp J, Bielmann R, Schmelcher M, Born Y, Zimmer M, et al. Comparative genome analysis of listeria bacteriophages reveals extensive mosaicism, programmed translational frameshifting, and a novel prophage insertion site. J Bacteriol. 2009;191(23):7206–15.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Mukhina YI, Kumeiko VV, Podgornaya OI, Efremova SM. The fate of larval flagellated cells during metamorphosis of the sponge Halisarca dujardini. Int J Dev Biol. 2006;5:533–41.

    Google Scholar 

  32. 32.

    Shaposhnikova T, Matveev I, Napara T, Podgornaya O. Mesogleal cells of the jellyfish Aurelia aurita are involved in the formation of mesogleal fibres. Cell Biol Int. 2005;29:952–8.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Matveev I, Shaposhnikova T, Podgornaya O. A novel Aurelia aurita protein mesoglein contains DSL and ZP domains. Gene. 2007;399:20–5.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Parrinello N. Focusing on Ciona intestinalis (Tunicata) innate immune system. Evolutionary implications. Invertebr Surviv J. 2009;6(1):S46–57.

    Google Scholar 

  35. 35.

    White SH, Wimley WC, Selsted ME. Structure, function, and membrane integration of defensins. Curr Opin Struct Biol. 1995;5(4):521–7.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Ding J, Chou Y-Y, Chang TL. Defensins in viral infections. J Innate Immun. 2009;1(5):413–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Wilson SS, Wiens ME, Smith JG. Antiviral mechanisms of human defensins. J Mol Biol. 2013;425(24):4965–80.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Sahl HG, Pag U, Bonness S, Wagner S, Antcheva N, Tossi A. Mammalian defensins:structures and mechanism of antibiotic activity. J Leukoc Biol. 2005;77(4):466–75.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Korndörfer IP, Kanitz A, Danzer J, Zimmer M, Loessner MJ, Skerra A. Structural analysis of the l-alanoyl-d-glutamate endopeptidase domain of Listeria bacteriophage endolysin Ply500 reveals a new member of the LAS peptidase family. Acta Crystallogr Sect D. 2008;64(6):644–50.

    Article  CAS  Google Scholar 

  40. 40.

    Loessner MJ, Wendlinger G, Scherer S. Heterogeneous endolysins in Listeria monocytogenes bacteriophages: a new class of enzymes and evidence for conserved holin genes within the siphoviral lysis cassettes. Mol Microbiol. 1995;16:1231–41.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Loessner MJ, Kramer K, Ebel F, Scherer S. C-terminal domains of Listeria monocytogenes bacteriophage murein hydrolases determine specific recognition and high-affinity binding to bacterial cell wall carbohydrates. Mol Microbiol. 2002;44(2):335–49.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Fukushima T, Yao Y, Kitajima T, Yamamoto H, Sekiguchi J. Characterization of new L, D-endopeptidase gene product CwlK (previous YcdD) that hydrolyzes peptidoglycan in Bacillus subtilis. Mol Genet Genomics. 2007;278:371–83.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Azumi K, Yokosawa H. Characterization of novel Metallo-proteases released from ascidian Hemocytes by treatment with calcium Ionophore. Zool Sci. 1996;13(3):365–70.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Azumi K, Yokosawa H. Characterization of protease-releasing factors isolated from hemocytes of the solitary ascidian, Halocynthia roretzi. Zool Sci. 1997;14(3):391–5.

    CAS  Article  Google Scholar 

  45. 45.

    Arolas JL, Bronsoms S, Ventura S, Avilés F, Calvete J. Characterizing the tick carboxypeptidase inhibitor - molecular basis for its two-domain nature. J Biol Chem. 2006;281:22906–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Guasch A, Coll M, Avilés FX, Huber R. Three-dimensional structure of porcine pancreatic procarboxypeptidase A. A comparison of the A and B zymogens and their determinants for inhibition and activation. J Mol Biol. 1992;224(1):141–57.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. 47.

    Van Wart HE, Birkedal-Hansen H. The cysteine switch: a principle of regulation of metalloproteinase activity with potential applicability to the entire matrix metalloproteinase gene family. Proc Natl Acad Sci U S A. 1990;87(14):5578–82.

    PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Odintsov SG, Sabala I, Marcyjaniak M, Bochtler M. Latent LytM at 1.3Å resolution. J Mol Biol. 2004;335(3):775–85.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. 49.

    Bochtler M, Odintsov SG, Marcyjaniak M, Sabala I. Similar active sites in lysostaphins and D-ala-D-ala metallopeptidases. Protein Sci. 2009;13(4):854–61.

    Article  CAS  Google Scholar 

  50. 50.

    Otsuka A, Dreier J, Cheng PF, Nägeli M, Lehmann H, Felderer L, et al. Hedgehog pathway inhibitors promote adaptive immune responses in basal cell carcinoma. Clin Cancer Res. 2015;21(6):1289 LP–1297.

    Article  CAS  Google Scholar 

  51. 51.

    Westendorp BF, Büller NVJA, Karpus ON, van Dop WA, Koster J, Versteeg R, et al. Indian hedgehog suppresses a stromal cell–driven intestinal immune response. Cell Mol Gastroenterol Hepatol. 2018;5(1):67–82.e1.

    PubMed  Article  PubMed Central  Google Scholar 

  52. 52.

    Fuse N, Maiti T, Wang B, Porter JA, Hall TM, Leahy DJ, et al. Sonic hedgehog protein signals not as a hydrolytic enzyme but as an apparent ligand for patched. Proc Natl Acad Sci U S A. 1999;96(20):10992–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Naranjo-Ortíz MA, Brock M, Brunke S, Hube B, Marcet-Houben M, Gabaldón T. Widespread inter- and intra-domain horizontal gene transfer of D-amino acid metabolism enzymes in eukaryotes. Front Microbiol. 2016;7:2001.

    PubMed  PubMed Central  Article  Google Scholar 

  54. 54.

    Andersson JO. Evolution of patchily distributed proteins shared between eukaryotes and prokaryotes: Dictyostelium as a case study. J Mol Microbiol Biotechnol. 2011;20(2):83–95.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  55. 55.

    Jackson DJ, Macis L, Reitner J, Wörheide G. A horizontal gene transfer supported the evolution of an early metazoan biomineralization strategy. BMC Evol Biol. 2011;11(1):238.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  56. 56.

    Syvanen M. Evolutionary implications of horizontal gene transfer. Annu Rev Genet. 2012;46:341–58.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  57. 57.

    Crisp A, Boschetti C, Perry M, Tunnacliffe A, Micklem G. Expression of multiple horizontally acquired genes is a hallmark of both vertebrate and invertebrate genomes. Genome Biol. 2015;16(1):50.

    PubMed  PubMed Central  Article  Google Scholar 

  58. 58.

    Grau-Bové X, Ruiz-Trillo I, Rodriguez-Pascual F. Origin and evolution of lysyl oxidases. Sci Rep. 2015;5:10568.

    PubMed  PubMed Central  Article  Google Scholar 

  59. 59.

    Davín AA, Tannier E, Williams TA, Boussau B, Daubin V, Szöllősi GJ. Gene transfers can date the tree of life. Nat Ecol Evol. 2018;2(5):904–9.

    PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Dellaporta SL, Xu A, Sagasser S, Jakob W, Moreno MA, Buss LW, et al. Mitochondrial genome of Trichoplax adhaerens supports Placozoa as the basal lower metazoan phylum. Proc Natl Acad Sci. 2006;103(23):8751 LP–8756.

    Article  CAS  Google Scholar 

  61. 61.

    Srivastava M, Begovic E, Chapman J, Putnam NH, Hellsten U, Kawashima T, et al. The Trichoplax genome and the nature of placozoans. Nature. 2008;454:955.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Schierwater B, Eitel M, DeSalle R. World Placozoa Database. Trichoplax Schultze, 1883. World Register of Marine Species. 2018. on 2018-07-18. Accessed 28 May 2018.

    Google Scholar 

  63. 63.

    Jo B-S, Choi SS. Introns: the functional benefits of introns in genomes. Genomics Inform. 2015;13(4):112–8.

    PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Rose TM, Schultz ER, Henikoff JG, Pietrokovski S, McCallum CM, Henikoff S. Consensus-degenerate hybrid oligonucleotide primers for amplification of distantly related sequences. Nucleic Acids Res. 1998;26(7):1628–35.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Miller RT, Kubier P, Reynolds B, Henry T, Turnbow H. Blocking of endogenous avidin-binding activity in immunohistochemistry: the use of skim milk as an economical and effective substitute for commercial biotin solutions. Appl Immunohistochem Mol Morphol. 1999;7:63–5.

    Google Scholar 

  66. 66.

    Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, et al. In: Walker JM, editor. Protein identification and analysis tools on the ExPASy server BT - the proteomics protocols handbook. Totowa: Humana Press; 2005. p. 571–607.

    Google Scholar 

  67. 67.

    Käll L, Krogh A, Sonnhammer ELL. A combined transmembrane topology and signal peptide prediction method. J Mol Biol. 2004;338(5):1027–36.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  68. 68.

    Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785.

    CAS  Article  Google Scholar 

  69. 69.

    Mooney C, Cessieux A, Shields DC, Pollastri G. SCL-Epred: a generalised de novo eukaryotic protein subcellular localisation predictor. Amino Acids. 2013;45(2):291–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  70. 70.

    Emanuelsson O, Brunak S, von Heijne G, Nielsen H. Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. 2007;2:953.

    CAS  PubMed  Article  Google Scholar 

  71. 71.

    George RA, Heringa J. The REPRO server: finding protein internal sequence repeats through the web. Trends Biochem Sci. 2000;25(10):515–7.

    CAS  PubMed  Article  Google Scholar 

  72. 72.

    Heger A, Holm L. Rapid automatic detection and alignment of repeats in protein sequences. Proteins. 2000;41(2):224–37.

    CAS  PubMed  Article  Google Scholar 

  73. 73.

    Mirabello C, Pollastri G. Porter, PaleAle 4.0: high-accuracy prediction of protein secondary structure and relative solvent accessibility. Bioinformatics. 2013;29(16):2056–8.

    CAS  PubMed  Article  Google Scholar 

  74. 74.

    Jones DT, Cozzetto D. DISOPRED3: precise disordered region predictions with annotated protein-binding activity. Bioinformatics. 2015;31(6):857–63.

    CAS  PubMed  Article  Google Scholar 

  75. 75.

    Hanson J, Yang Y, Paliwal K, Zhou Y. Improving protein disorder prediction by deep bidirectional long short-term memory recurrent neural networks. Bioinformatics. 2017;33(5):685–92.

    PubMed  Google Scholar 

  76. 76.

    Cilia E, Pancsa R, Tompa P, Lenaerts T, Vranken WF. The DynaMine webserver: predicting protein dynamics from sequence. Nucleic Acids Res. 2014;42(W1):W264–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  77. 77.

    Remmert M, Biegert A, Hauser A, Söding J. HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat Methods. 2011;9:173.

    PubMed  Article  CAS  Google Scholar 

  78. 78.

    González-Domínguez J, Liu Y, Touriño J, Schmidt B. MSAProbs-MPI: parallel multiple sequence aligner for distributed-memory systems. Bioinformatics. 2016;32(24):3826–8.

    PubMed  Article  CAS  Google Scholar 

  79. 79.

    Söding J, Biegert A, Lupas AN. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005;33:W244–8.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  80. 80.

    Goodstadt L, Ponting CP. CHROMA: consensus-based colouring of multiple alignments for publication. Bioinformatics. 2001;17(9):845–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  81. 81.

    deMayo JA, Maas KR, Klassen JL, Balunas MJ. Draft genome sequence of Streptomyces sp. AVP053U2 isolated from Styela clava, a tunicate collected in long island sound. Genome Announc. 2016;4(5):e00874–16.

    PubMed  PubMed Central  Article  Google Scholar 

  82. 82.

    Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal omega. Mol Syst Biol. 2011;7:539.

    PubMed  PubMed Central  Article  Google Scholar 

  83. 83.

    Rice P, Longden I, Bleasby A. EMBOSS: the european molecular biology open software suite. Trends Genet. 2000;16(6):276–7.

    CAS  Article  Google Scholar 

  84. 84.

    Korber B. HIV signature and sequence variation analysis. In: Rodrigo AG, Learn GH, editors. Computational analysis of HIV molecular sequences. Dordrecht: Kluwer Academic Publishers; 2000. p. 55–72.

    Google Scholar 

Download references


The authors greatly appreciate the help received at the Kartesh White Sea Biological Station of the Zoological Institute of the Russian Academy of Sciences. We used the core facilities of the Research Park of St. Petersburg State University: Center for Molecular and Cell Technologies, Center for Microscopy and Microanalysis, and Observatory of Environmental Safety Center. We would also like to thank Alexey Gurevich for help with bioinformatics analysis and Laurel Sky Hiebert for help with text editing.


This work was supported by the “Molecular and Cell Biology” program of the Presidium of the Russian Academy of Sciences (grant no. 01.2.01457147) and the Russian Foundation for Basic Research (grant no. 15–04-06008-а).

Availability of data and materials

The datasets used and/or analysed during the current study are available at public databases GenBank (, UniProtKB (, PDB (, SCOP (, and Pfam 30.0 ( or included in this published article (and its supplementary information files).

Author information




MD, SS, AS, TS and LA performed the experiments. OP designed the study. All authors discussed the results and contributed to the final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Maria A. Daugavet.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Daugavet, M.A., Shabelnikov, S., Shumeev, A. et al. Features of a novel protein, rusticalin, from the ascidian Styela rustica reveal ancestral horizontal gene transfer event. Mobile DNA 10, 4 (2019).

Download citation


  • Ascidians
  • Bacteriophage
  • Hemocytes
  • Horizontal gene transfer
  • L-alanyl-D-glutamate peptidase
  • Trichoplax
  • tRNA