1. Introduction
Serine proteases constitute the major granule content of several hematopoietic cell types. Large amounts of active chymotryptic and tryptic enzymes are stored in cytoplasmic granules of mammalian mast cells for rapid released upon activation. Mammalian neutrophils instead store primarily elastolytic enzymes, as exemplified by neutrophil elastase and proteinase 3, but also a tryptic enzyme, NSP-4 and a chymotryptic enzyme, cathepsin G [
1,
2,
3,
4]. Human cytotoxic T cells and natural killer cells (NK cells) express and store several granzymes with tryptase, asp-ase and met-ase specificities, as exemplified by granzymes A, B and M respectively [
1]. A relatively detailed picture has been established of how these enzymes have appeared and diversified during mammalian evolution [
2,
3,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15]. However, the information concerning the presence and specificity of such enzymes in reptiles, amphibians, and fishes is still only fragmentary [
16,
17,
18,
19,
20,
21,
22,
23]. There has been an increase in the use of non-rodent animal models with much focus on non-mammalian vertebrates including fish species, notably the zebrafish. However, there are few studies on the details regarding the characterisation of hematopoietic serine proteases in other fishes. Here, the (spotted) gar is of particular interest as it represents an basally diverging branch of the ray-finned fish evolutionary tree and has been found to have a genome that has undergone relatively limited number of rearrangements and amplifications compared to many other fish lineages including the zebrafish, which is often used as a model organism [
24]. The gar genome thereby represents a less derived configuration than most other fish genomes and has interestingly been found show major similarities to mammalian genomes [
24]. In order to look closer into early events in the expansion and diversification of these hematopoietic serine protease, this analysis focuses on two hematopoietic serine proteases from the spotted gar.
In mammals, serine proteases play important and diverse roles in a number of physiological processes including blood coagulation, food digestion, fertilisation, immunity and tissue repair [
25]. This large chymotrypsin/trypsin family all share a common mechanism for cleaving peptide bonds, based on their catalytic triad, with three vital residues: His57, Asp102 and Ser195 (chymotrypsinogen numbering) [
26] (
Figure 1A). These key amino acids are located near a substrate binding pocket (termed S1), typically made up of residues 189, 216 and 226 (also chymotrypsinogen numbering) [
26] (
Figure 1B). Together they form the specificity conferring triplet and provide clues as to the primary specificity of the serine proteases (
Figure 1B). For mammalian enzymes and also other tetrapods, these three residues provide a relatively good indication for the primary specificity of an enzyme. However, for the majority of fish proteases, excluding the granzyme A/K members, these three residues give little guidance to their primary specificities. This is primarily due to large sequence divergences making the positioning of the relevant residues difficult to determine based only on the primary sequence. To determine the specificity of the fish proteases we therefore have to rely on direct experimental analysis.
The primary specificity of a protease is the amino acid after which the cleavage occurs and this residue is named the P1 residue (
Figure 1A). The residues N-terminally of the cleavage site are numbered P2, P3, P4 etc and the residues C-terminally of the cleavage site P1´, P2´, P3´ etc (
Figure 1A). The primary specificity determines the main specificity of the enzyme. If the enzyme cleaves after a basic amino acid, such as Arg or Lys, it has tryptic activity, when it cleaves after large aromatic amino acids, such as Phe, Tyr or Trp it has a chymotryptic activity (
Figure 1B). The extended specificity is determined by the selectivity in residues surrounding the cleavage site, most often the P2, P3, P4, P1´, P2´and P3´residues (
Figure 1A).
To determine the cleavage specificity of an enzyme different techniques can be used, including chromogenic substrates, peptide libraries and phage display. Substrate phage display is a method that can provide detailed information concerning both the primary and extended specificity of an enzyme. This method has been used to study the extended specificities of two serine proteases from the spotted gar. These two proteases were extracted from the genomic database depending on the distance of the relationship with other known hematopoietic serine proteases. Both enzymes showed very stringent both primary and extended specificities. One being a Tyr-specific chymase and the other being a Leu-specific Leu-ase. Antisera developed against these two proteases detected cells in Peyer´s patch-like structures of the intestinal region of the gar, with a pattern of expression indicating that the proteases are expressed in cytotoxic T cells or NK cells. However, it is yet not possible to more specifically identify the cell origin as very few reagents are available for studies of immune cells in the gar.
3. Discussion
Hematopoietic serine proteases perform a number of important functions in mammalian immunity. Granzyme B, which is an enzyme found in cytotoxic T cells and NK cells, induces apoptosis in target cells via cleavage of caspases [
1,
28,
29]. Mast cells use the chymases to regulate blood pressure via angiotensin cleavage, to regulate types of immunity by cleaving selective sets of cytokines, to combat blood-feeding parasites such as mosquitos and ticks via the cleavage of anti-coagulant proteins and to inactivate toxins from snakes and scorpions [
4,
30,
31,
32,
33,
34,
35,
36,
37,
38,
39]. Neutrophil proteases help these cells to migrate through tissues to reach the site of infection and cleave the pathogen-associated molecules of the infectious organisms [
1]. However, which functions the corresponding proteases have among fishes and what their targets are, are almost completely unknown. Relatively solid evidence for granzyme A/K equivalents in fish are present. However, the function of mammalian granzyme A and K are have not been well elucidated [
40]. However, they are the most well conserved of all hematopoietic serine proteases and are found in essentially all species from cartilaginous fish to humans but major question still remain concerning their primary targets [
40]. Major similarities are also found among different tetrapods including amphibians, reptiles and birds, and mammals concerning cleavage specificities and cell origin of other hematopoietic serine proteases. However, in ray-finned fishes we know very little concerning this large family of proteases. Almost all of the non-granzyme A/K members of the hematopoietic serine proteases form their own branch in the phylogenetic tree, indicating that they have experienced an independent expansion and diversification compared to tetrapods (
Figure 2). The structure of the ray-finned fish enzymes are also sufficiently different from the mammalian enzymes, so it’s not possible to use the positions of the three residues that form the S1 pocket of the tetrapod enzymes to obtain information concerning the primary specificity of the fish proteases [
1,
2].
What is known about hematopoietic serine proteases in fishes? Granzyme A from the catfish has previously been shown to have cytolytic activity [
23]. However, caution around this interpretation maybe needed as the human enzyme, despite apoptotic activity being reported, when more physiological amounts of the enzyme were used this activity was lost [
40,
41,
42]. A more detailed analysis of this phenomenon maybe needed before definitive conclusions can be drawn. A second catfish enzyme, catfish granzyme like-I, has recently been shown to cleave a sequence in catfish caspase 6 that corresponds to the region in human caspase 3 that is cleaved by human granzyme B [
17]. However, also here no direct proof is yet available for that this cleavage occurs in vivo. The analysis has here been hampered by difficulties in producing recombinant catfish caspase 6, which is needed to enable cleavage analysis with catfish granzyme-like I. However, this enzyme is expressed by catfish NK-like cells indicating that it may be involved in apoptosis induction of cells infected with intracellular parasites. Detailed information concerning the cleavage specificities of three additional fish hematopoietic serine proteases are now available, catfish granzyme-like II, which is also expressed by the same NK-like cells and the two enzymes from the spotted gar presented here [
43]. The catfish granzyme-like II is a highly specific elastase cleaving after Ala with a strong preference for several basic amino acids, primarily Arg, just N-terminally of the cleavage site [
43]. To our knowledge no elastase is expressed by human NK cells, which is why the function of this enzyme may turn out to be a fish-specific immune function. The two gar enzymes are, based on their tissue location and cell shape, expressed in what appears to be NK-like cells or cytotoxic T cells. However, it is not yet clear in what cell these two proteases are expressed due to the lack of suitable reagents to determine the actual cellular origin in gar. Their specificities currently provide little clues to their function. Both enzymes are highly specific. One is a Tyr-specific chymase and the second enzyme a Leu- specific Leu-ase, similar to the mast cell chymases of rabbit and guinea pig, but both appears to be expressed by a different cell type than these latter proteases (
Figure 4). Neither chymases nor Leu-ases seem to be expressed by NK or cytotoxic T cells in mammals. Screening of the total proteome with the consensus sequences for these two gar proteases did not result in any clear candidates for their potential immune functions. A relatively broad panel of different potential targets were obtained for both enzymes with no clear identified immune candidate. A more in-depth analysis involving direct cleavage of a cell extract and possibly 2-D gel analysis of cleavage products could be informative. An alternative explanation could be that these enzymes are not aimed at cleaving host proteins but instead pathogen derived proteins such as bacterial toxins.
The analysis of the intestinal region and the head kidney resulted in an interesting finding concerning the immune organs in this branch of fish evolution. Clear indications for Peyer’s patch-like structures in the intestinal region were observed. Here, the expression of the two gar enzymes was seen, which did not occur in other regions of the fish, including the head kidney, indicating that the proteases are expressed in cells that appear at these sites of immune activity in the intestinal region and not during early expansion of cells in the head kidney.
A clear separation of regions where red blood cells is formed was seen from the gar head kidney sections and other regions within the same organ where most likely the majority of other hematopoietic cells are residing was also an interesting finding that differs from the organization of human bone marrow. The clusters of red blood cells in the regions of other hematopoietic cells may be due to a rich blood supply to support rapidly expanding cell populations in need for oxygen and nutrient supply.
The analyses of the genomic organization of the two gar enzymes also gave insight into the evolution of these proteases. Both of the gar enzymes were located in the met-ase locus (
Figure 9). This locus is known in mammals to encode the majority of neutrophil proteases including N-elastase, proteinase 3, azurocidin and neutrophil protease 4, but also the NK-cell expressed granzyme M and complement factor D. We also observed that independent amplifications had occurred in mammals compared to gar and sterlet sturgeon (
Figure 9). In mammals it is most likely complement factor D that is the origin of the majority of the neutrophil proteases, including N-elastase, proteinase 3 and azurocidin, whereas in the gar it is most likely the PRSS57, encoding NSP4, which has been duplicated forming a small subfamily of six different proteases in the gar, two in sturgeon, and three in the bichir (
Figure 9). These three fish species are representatives of non-teleost branches of the fish evolutionary tree, indicating that the gene duplications of PRSS57 are early expansions in the met-ase locus in fish [
44]. In this regard, the relation and differences between the teleost catfish and zebrafish, and the non-teleost gar, sterlet, and bichir is also interesting. The one or two DDN-like enzymes in one end of the locus found in all three non-teleost fish branches including gar, sterlet and bichir have been lost in both catfish and zebrafish (
Figure 9). The three genes that most likely originate from gene duplications of PRSS57 are still present in the teleosts catfish and zebrafish (
Figure 9). Both the catfish granzyme-like II and gar granzyme G are located as the first gene of these PRSS57 duplicated genes, however they have very different cleavage specificities, where catfish granzyme-like II is an Ala-specific elastase and the gar enzyme is a Tyr- specific chymase. Currently there is no information concerning the other catfish and gar enzymes in this cluster except for catfish granzyme-like I, which is a met-ase. There are striking similarities between catfish and zebrafish as both catfish granzyme-like I and zebrafish AE-like, as well as also catfish granzyme-like II and zebrafish SPA have very similar specificities [
43]. Future studies will hopefully be able to distinguish the similarities and differences in specificities between these enzymes in non-teleost compared to teleost fishes. Such information may guide us into key primary immune-related targets for these enzymes during fish evolution and how this relates to their mammalian counterparts.
4. Materials and Methods
4.1. Phylogenetic Analyses
The phylogenetic analyses aimed at determining the relationship between the two gar enzymes and other hematopoietic serine proteases from fish were performed essentially as described in a previous publication using the same strategy and sequences [
2]. Sequences relating to the two gar enzymes were systematically uncovered by BLASTp searching of all animal NCBI databases. The mature two gar enzymes were used as the query sequence and all novel derived sequences were analysed using the multiple alignment programme MAFFT with G-INS-i strategy and default parameters to determine whether they belonged to the serine protease family. To visualise the relationship between the two gar enzymes with those from other species, a phylogenetic tree using the Bayesian interference of phylogeny algorithm with posterior probabilities in the MRBAYES program was constructed and viewed in FigTree (v1.4). The amino acid sequences for mature proteins of serine proteases branching with the two gar enzymes were aligned using MAFFT.
4.2. Production and Purification of Recombinant Gar Gzm-G and Gar DDN1
The gar granzyme sequences (GenBank accession numbers: Gar Gzm-G-like (XP_015220993) and Gar DDN1-like (XP_006640095)) were designed and ordered from GenScript (Piscataway, NJ, USA). The synthesized construct was cloned in the pU57 cloning vector, containing EcoRI and XhoI sites. The gar granzyme sequences were subsequently transferred to a pCEP-Pu2 vector, used for expression in mammalian cells [
27]. The enzymes were produced as inactive recombinant proteins, with an N-terminal His
6-tag followed by an enterokinase (EK) site. HEK 293 cells were grown to 70% confluency in a 25 cm
3 tissue culture flask (BD VWR) with Dulbecco’s Modified Eagles Medium (DMEM) (GlutaMAX, Invitrogen, Carlsbad, CA, USA) supplemented with 5% fetal bovine serum (FBS) and 50 µg/ml gentamicin. Following transfection with lipofectamine (Invitrogen, Carlsbad, CA, USA), using approximately 25 µg of the gar constructs in pCEP-Pu2, puromycin was added to the DMEM (0.5 µg/mL) to select for cells which had taken up the DNA. Heparin (5 µg/ml) was added to the culture medium to enhance recovery of secreted protein. Cells were expanded and conditioned media collected.
After collecting sufficient amount of media, in general around 750 mL, the conditioned media was filtered (Munktell 00H 150 mm, Falun, Sweden) and 500 µL nickel nitrilotriacetic acid (Ni-NTA) beads were added (Qiagen, Hilden, Germany) to purify the recombinant enzymes. The media with Ni-NTA beads were rotated for 45 mins at 4 °C. Subsequently, the Ni-NTA beads were collected by centrifugation and transferred to a column containing a glass filter (Sartorius, Goettingen, Germany). To remove unbound protein the columns were washed with PBS tween 0.05 % + 10 mM imidazole + 1 M NaCl. Following this wash the recombinant protein was eluted in PBS tween 0.05 % + 100 mM imidazole fractions. The first fraction volume was half the Ni-NTA bead width (200 µL) and further fractions eluted with a full bead width (400 µL). Individual fractions were run on SDS-PAGE gel, their concentrations estimated from a bovine serum albumin standard (BSA) and the most concentrated were pooled and kept at 4 °C.
4.3. Activation of Recombinant Gar Granzymes
To activate the enzymes the initial concentrations of the gar enzymes were first determined by SDS-PAGE and the level of EK (Roche, Mannheim, Germany) adjusted for activation of the enzymes. To 70 µL of the eluted recombinant enzyme we added 1 µL EK and incubated for 3 hrs at 37 °C to activate the enzyme. The activated fractions were stored at 4 °C until use.
4.4. Substrate Phage Display
A T7 phage library containing 5 x 107 variants was used for the phage display analysis to determine the extended cleavage specificity of the two gar enzymes. In this library each phage displays a unique nine amino acid sequence. The nine amino acid region had been inserted into the C-terminal of the capsid 10 protein, followed by a His6-tag. 125 µL Ni-NTA agarose beads via their His6-tags for 1 hr at 4°C with gentle rotation (Qiagen, Hilden, Germany). To these 125 ul Ni-NTA we estimate that approximately 109 plaque forming units (pfu) were bound. Unbound phages were removed by washing ten times with PBS tween 0.05 % + 1 M NaCl, followed by two washes with PBS. The beads were re-suspended in 375 µL PBS and approximately 250 ng of the recombinant gar enzyme was added. This reaction was incubated overnight or for approximately 16 hours at 37 °C with gentle rotation. The enzyme cleavage results in that susceptible phages detach from the Ni-NTA beads. From the supernatant released phages was recovered after centrifugation. Thirty µL was used in a plaque assay to determine the number of released phages. Briefly, ten-fold serial dilutions were made, mixed with (E.coli) BLT5615 (for propagation and visualization of plaques on a bacterial lawn) and plated on LA-Amp (50 µg/mL) plates, incubated for 2.5 hrs at 37 °C and then counted. The remaining supernatant was added to 10 ml BLT5615 bacteria (OD600 0.5) and the culture was incubated at 37 °C for approximately 75 mins until the culture had cleared by phage lysis for phage expansion. From this, 1.5 mL was centrifuged to remove bacterial debris and 800 µL transferred to a microcentrifuge tube containing 100 µL PBS and 100 µL 5M NaCl. This solution was bound to 125 µL Ni-NTA beads directly after centrifugation to start the next selection cycle. The complete process was repeated a further 6 times, constituting 7 selection rounds. Individual plaques were isolated from the final selection round in 100 µL phage buffer before vortexing for 30 mins and stored at 4 °C. The random nine amino acid regions contained in these phages were amplified by PCR (T7Select primers, Novagen, Sacramento, CA, USA) and sequenced by Eurofins (Sequencing centre, Cologne, Germany). The resulting sequences were translated using CLC viewer and aligned using Adobe Illustrator. A parallel control reaction without enzyme (only PBS) was also run under the same conditions and plaque numbers were compared to the enzyme sample.
4.5. Phage Display Sequence Verification using a Two-Thioredoxin (trx) Approach
In order to verify the phage display data, a recombinant trx system was used. This system has been developed in our laboratory and used multiple times to study other proteases with great success. Here, the random nine amino acid cleaved region determined from the phage display was introduced between two adjacent E. coli trx proteins. Originally, a pET21 vector containing a single trx protein was modified to contain a second trx with BamHI and SalI sites in the intervening region. Here the random region was synthesised as oligonucleotides (Sigma, St Louis, USA). These oligos generated overhand suitable for ligation between BamHI and SalI restriction sites of the vector. This resulted in a vector containing a first trx followed by the random cleavable region then a second trx with His6-tag (facilitating purification).
To produce the recombinant protein the construct was expressed in E.coli Rosetta gami (Novagen, Sacramento, CA, USA). Ten ml of an overnight culture was added to 90 mL LB+Amp (50 µg/mL) and 500 µL 20 % glucose. After approximately 1 hr (reaching OD600 0.5), 100 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) was added and the culture placed on a shaker (moderate shaking) at 37 °C for 3 hrs. The bacterial cells of the culture were pelleted by centrifugation at 10,000 rpm for 3 mins and supernatant discarded. The pellet containing the bacteria with the expressed recombinant protein was washed in 10 mL PBS tween 0.05 %, centrifuged and pelleted again followed by resuspension in 1/100th starting volume (i.e., 1 mL) PBS. To obtain the intracellularly expressed protein, the resuspended pellet was sonicated for 6 x 30 secs on ice. The supernatant was transferred to a new microcentrifuge tube after centrifugation at 10,000 rpm for 10 mins at 4 °C. To purify the protein, 125 µL Ni-NTA beads were added and incubated for 45 mins at 4 °C with gentle agitation (Qiagen, Hilden Germany). The solution with Ni-NTA beads was then transferred to a 2 mL column (Terumo, Leuven, Belgium) containing a glass filter (Sartorius, Goettingen, Germany). The column was then first washed with 3 x 2 mL and 2 x 1 mL PBS tween 0.05 % + 10 mM imidazole. To elute the protein, fractions were collected after passing through PBS tween 0.05 % + 100 mM imidazole. The first fraction volume was only half the Ni-NTA bead volume (75 µL) and further fractions eluted with a full bead volume (150 µL). Individual fractions were run on SDS-PAGE gel and their concentrations estimated from a BSA standard (and Bradford assay (Bio-rad, CA, USA)). The most concentrated fractions were then pooled and kept at 4 °C. until use in the assay
For cleavage analysis, approximately 250 ng of recombinant gar enzyme was added to 20 µg of the pooled 2-trx protein (containing different sequences based on the phage display data) and aliquots of 5 µg removed after 0, 15, 45 and 150 mins after enzyme addition. The reactions were run at room temperature and the time point aliquots analysed on SDS-PAGE gel under denaturing conditions using pre-cast 4-12 % Bis-Tris gels (Invitrogen, Carlsbad, CA, USA) and 1x MES buffer (Invitrogen, Carlsbad, CA, USA). Gels were stained with Coomassie colloidal solution to visualize the protein bands [
45].
4.6. Screening for Potential In Vivo Substrates
The derived consensus sequences for the two gar proteases (Ala-Val-Phe-Tyr-Ser-Leu) and (Arg-Arg-Leu-Pro-Leu-Leu), respectively were used in a screening for potential in vivo substrates using a standard protein Basic Local Alignment Search Tool (BLASTp).
4.7. Spotted Gar Husbandry and Tissue Sampling
Spotted gars were raised to juvenile stages from fertilized eggs obtained from hormone-induce spawns of wild-caught broodstock from bayous near Thibodaux, Louisiana. Five unsexed individuals were euthanized at the age of 4 months (average total length: 15.9cm) using MS-222 (250–500 mg/liter water) and fixed overnight in 4% Paraformaldehype/PBS. All animal work was approved by the Michigan State University Institutional Animal Care and Use Committee (PROTO202200298).
4.8. Immunohistochemistry
Approximately 100 µg of the recombinant enzymes were mixed with Freund’s complete adjuvant and injected in one male 20-week Sprague Dawley rat for each enzyme. Three weeks after the first injection the rats were boosted with 20 µg of the recombinant protein in Freunds incomplete adjuvant. After 6 weeks another booster dose was given with 20 µg of recombinant protein in Freund’s incomplete adjuvant. Two weeks after the last booster injection the rats were tail blead to obtain approximately 1 mL of blood. The blood was left overnight to coagulate and thereafter spun at 3000 rpm in an Eppendorf centrifuge and the serum was transferred to a new tube and centrifuged a second time to remove remaining red blood cells. The serum was then aliquoted and kept at -80oC until use in immunohistochemical staining.
The fish were preserved in a formalin solution for several days. Following preservation, the tissue underwent decalcification using 0.5 M ethylenediaminetetraacetic acid at pH 8.0 (Sigma Aldrich, Saint Louis, Missouri, USA) for a duration of one month, during which the solution was replaced on a weekly basis. Subsequently, the tissue was immersed in a 30% sucrose solution in 1x PBS for 72 hours before sectioning. A transverse incision was then performed at the tissue’s midpoint, extracting a section of approximately 1 cm² around the intestinal cavity with the use of a scalpel. The dissected tissue was then embedded in optical cutting compound (OCT, VWR, Radnor, Pennsylvania, USA). Sections of 15 µm thickness were obtained using a cryostat and mounted onto SuperFrost slides (VWR, Radnor, Pennsylvania, USA).
For immunostaining, the tissue sections were initially blocked for a period of 1 hour at room temperature with 1x protein-free blocking buffer (Pierce solution, Thermo Fisher, Waltham, Massachusetts, USA) containing 0.3 % Triton X-100 in 0.1 M PBS. Staining with the rat anti gar-DDN1 and anti-GzmG, (1:100) primary antibody was performed overnight at 4oC. After the primary antibody incubation, the tissue sections underwent triple washing with a PBS solution. Subsequently, the sections were incubated with a donkey anti-rat secondary antibody conjugated to Alexa 488 (1:500, Thermo Fisher, Waltham, Massachusetts, USA) for 2 hours at room temperature. Nuclei were visualized using DAPI (4′,6-diamidino-2-phenylindole)(Thermo Fisher, Waltham, Massachusetts, USA) and mounted using Fluoromount G (Thermo Fisher, Waltham, Massachusetts, USA). Pictures were taken in an epifluorescence Nikon microscope (Nikon Eclipse 90i Upright).