Deiminated proteins in extracellular vesicles and serum of llama (Lama glama)-Novel insights into camelid immunity.

Peptidylarginine deiminases (PADs) are phylogenetically conserved calcium-dependent enzymes which post-translationally convert arginine into citrulline in target proteins in an irreversible manner, causing functional and structural changes in target proteins. Protein deimination causes generation of neo-epitopes, affects gene regulation and also allows for protein moonlighting. Furthermore, PADs have been found to be a phylogenetically conserved regulator for extracellular vesicle (EVs) release. EVs are found in most body fluids and participate in cellular communication via transfer of cargo proteins and genetic material. In this study, post-translationally deiminated proteins in serum and serum-EVs are described for the first time in camelids, using the llama (Lama glama L. 1758) as a model animal. We report a poly-dispersed population of llama serum EVs, positive for phylogenetically conserved EV-specific markers and characterised by TEM. In serum, 103 deiminated proteins were overall identified, including key immune and metabolic mediators including complement components, immunoglobulin-based nanobodies, adiponectin and heat shock proteins. In serum, 60 deiminated proteins were identified that were not in EVs, and 25 deiminated proteins were found to be unique to EVs, with 43 shared deiminated protein hits between both serum and EVs. Deiminated histone H3, a marker of neutrophil extracellular trap formation, was also detected in llama serum. PAD homologues were identified in llama serum by Western blotting, via cross reaction with human PAD antibodies, and detected at an expected 70 kDa size. This is the first report of deiminated proteins in serum and EVs of a camelid species, highlighting a hitherto unrecognized post-translational modification in key immune and metabolic proteins in camelids, which may be translatable to and inform a range of human metabolic and inflammatory pathologies.

Introduction

Lamoids, or llamas, belong to a family of camelids which are economically important livestock and have developed complex features and immunological traits related to their habitat (Wu et al., 2014; Saadeldin et al., 2018a). The llama (Lama glama), Bactrian camel (Camelus bactrianus), dromedary (Camelus dromedarius) and alpaca (Vicugna pacos) differ in their habitat. The Bactrian camel and dromedary are adapted to arid-desert-adapted environments, alpacas to plateaus, and the llama to higher altitudes. Llamas are historically found in the Andean highlands, specifically the Altiplano of southeast Perú and western Bolivia, as well as in Chile and Argentina, which has the third highest population of llamas (Wu et al., 2014). The domesticated llama is closely related to two extant wild South American camelids, the vicuña (Vicugna vicugna) and guanaco (Lama guanicoe). Previous genomic studies have revealed a range of specific adaptions in camelids relating to fat and water metabolism, osmoregulation, blood glucose level regulation, stress responses to heat, aridity, as well as to intense ultraviolet radiation and dust (Wu et al., 2014). Furthermore, a particularly important feature in camelid immunity is the production of small, homodimeric heavy chain-only, antibodies (HCAbs) which are of great value for the biomedical industry (Henry et al., 2019).

Peptidylarginine deiminases (PADs) are phylogenetically conserved calcium-dependent enzymes which post-translationally convert arginine into citrulline in target proteins in an irreversible manner. This can cause functional and structural changes in target proteins (Vossenaar, 2003; György et al., 2006; Wang and Wang, 2013; Bicker and Thompson, 2013). Structures most prone to deimination are beta-sheets and intrinsically disordered proteins, while the position of the arginine is also important; arginines sitting next to aspartic acid residues are most prone to citrullination, arginines next to glutamic acid residues are rarely citrullinated and those flanked by proline are poorly citrullinated (Nomura, 1992; Tarcsa et al., 1996; György et al., 2006). Protein deimination can affect gene regulation, cause generation of neoepitopes (Witalison et al., 2015; Lange et al., 2017) and may also allow for protein moonlighting, an evolutionary acquired phenomenon facilitating proteins to exhibit several physiologically relevant functions within one polypeptide chain (Henderson and Martin, 2014; Jeffrey, 2018; Magnadóttir et al., 2018a). PADs have been identified throughout phylogeny from bacteria to mammals, with 5 tissue specific PAD isozymes in mammals, 3 in chicken, 1 in bony fish and arginine deiminase homologues in bacteria (Vossenaar et al., 2003; Rebl et al., 2010; Magnadóttir et al., 2018a, a; Kosgodage et al., 2019). While studies on PADs in relation to human pathophysiology, including cancer, autoimmune and neurodegenerative diseases (Wang and Wang, 2013; Witalison et al., 2015; Lange et al., 2017; Kosgodage et al., 2017 & 2018) and CNS regeneration (Lange et al., 2011 and 2014) exist, relatively little phylogenetic research has been carried out on PADs in relation to normal physiology and evolutionary acquired adaptions of the immune system. Recent comparative studies focussing on roles for PADs in teleost fish have identified post-translational deimination in key proteins of innate, adaptive and mucosal immunity (Magnadóttir et al., 2018a, b; Magnadottir et al., 2019a; Magnadóttir et al., 2019b). A recent study in shark also revealed novel insights into this post-translational modification in relation to key immune factors, including shark immunoglobulins (Criscitiello et al., 2019). As the camelid family has developed unusual small immunoglobulins similar as to shark through convergent evolution, indicating common factors between shark and camelid immunity, we felt that an investigation of post-translationally deiminated proteins in camelids was warranted.

As PADs have been identified to be a key regulator of extracellular (EV)-release, a mechanism that has been found to be phylogenetically conserved from bacteria to mammals (Kholia et al., 2015; Kosgodage et al., 2017, 2018; Gavinho et al., 2019; Kosgodage et al., 2019), the characterisation of EVs in camelids is of further interest. Extracellular vesicles (EVs) are found in most body fluids and participate in cellular communication via transfer of cargo proteins and genetic material (Inal et al., 2013; Colombo et al., 2014; Lange et al., 2017; Turchinovich et al., 2019; Vagner et al., 2019). EVs in body fluids, including serum, can also be useful biomarkers to reflect health status (Hessvik and Llorente, 2018; Ramirez et al., 2018). Previous work on EVs has hitherto mainly been in the context of human pathologies, while recently comparative studies are growing (Iliev et al., 2018; Yang et al., 2015; Magnadóttir et al., 2019b; Criscitiello et al., 2019; Gavinho et al., 2019; Kosgodage et al., 2019). Few studies have been performed on EVs in camelids but therapeutic effects of EVs isolated from camel milk have been identified in halting cancer progression (Badawy et al., 2018). A recent study in shark identified for the first time deiminated small immunoglobulin proteins as part of EV cargo (Criscitiello et al., 2019). Due to the link between camelid and shark immunity through convergent evolution including the unusual immunoglobulin structure of small heavy chain-only Ig’s, a comparative study on EVs and deiminated EV cargo in a camelid species may provide further insights into shared immunological traits. Camelids have an unusual Ig repertoire and a large diversity of functional nanobodies has been identified in the llama (Harmsen and De Haard, 2007; Deschaght et al., 2017). This has made camelids an important source for small immunoglobulins that can be used for immunotherapy purposes, including for tumour targeting (van Lith et al., 2016), as well as for assessing cancer metastasis (Ramos-Gomes et al., 2018). As these nanobodies can also penetrate the blood-brain barrier (Širochmanová et al., 2018) they are of great value for a range of therapeutic treatment applications, including for brain cancers (Iqbal et al., 2010). Furthermore, as the camelid family has acquired unique metabolic features, they are also of interest as a model species for informing metabolic diseases.

In the current study we assessed post-translationally deiminated proteins in llama serum and serum-derived EVs, and report for the first time EV-mediated export of deiminated key immune, metabolic and nuclear proteins in serum of a camelid species.

Materials and methods

Animals and sampling

Llama (Lama glama L. 1758) serum was shared from excess blood collected in routine health checks of a resident male llama at the Texas A&M Winnie Carter Wildlife Center. Blood collected from the jugular vein of this 21 year old llama was allowed to clot at room temperature for 2 h before serum was collected by centrifuging at 300 g for 10 min. Serum was aliquoted and immediately frozen at −80 °C until further use.

Extracellular vesicle (EV) isolation and nanoparticle tracking analysis (NTA)

EVs were isolated by step-wise centrifugation according to established protocols using ultracentrifugation and the recommendations of MISEV2018 (the minimal information for studies of extracellular vesicles 2018; Théry et al., 2018). Llama serum was diluted 1:5 in ultrafiltered (using a 0.22 μm filter) Dulbecco’s PBS (DPBS, 100 μl serum added to 400 μl DPBS) and then centrifuged at 4000 g for 30 min at 4 °C for removal of cells and cell debris. The supernatant was collected and centrifuged at 100,000 g for 1 h at 4 °C. The pellet was then resuspended in DPBS and washed again at 100,000 g for 1 h at 4 °C. The resulting EV-enriched pellet was resuspended in 100 μl DPBS, diluted 1/100 in DPBS and analysed by NTA, based on Brownian motion of particles in suspension, using the NanoSight NS300 system (Malvern, UK). The NanoSight was used in conjunction with a syringe pump to ensure continuous flow of the sample, with approximately 40–60 particles per frame and videos were recorded for 5 x 90 s. The replicate histograms generated from the recordings were averaged.

Transmission electron microscopy (TEM)

EVs were isolated from serum as described above, the EV pellets were fixed with 2.5 % glutaraldehyde in 100 mM sodium cacodylate buffer (pH 7.0) for 1 h at 4 °C, resuspended in 100 mM sodium cacodylate buffer (pH 7.0), placed on to a grid with a glow discharged carbon support film, stained with 2 % aqueous uranyl acetate (Sigma-Aldrich) and thereafter viewed in TEM. Imaging was performed using a JEOL JEM 1400 transmission electron microscope (JEOL, Japan) operated at 80 kV at a magnification of 80,000–100,000. Digital images were recorded using an AMT XR60 CCD camera (Deben, UK).

Western blotting

Llama serum and EV isolates (an EV pellet derived from 100 μl serum, reconstituted in 100 μl DPBS after isolation and purification) were diluted 1:1 in 2x Laemmli sample buffer, boiled for 5 min at 100 °C and separated by SDS-PAGE on 4–20 % gradient TGX gels (BioRad UK). Approximately 5 μg protein was loaded per lane and transferred to nitrocellulose membranes using semi-dry Western blotting. Blocking of membranes was performed in 5 % BSA in TBS-T for 1 h at room temperature (RT) and incubation with primary antibodies, diluted in TBS-T, was carried out at 4 °C overnight (F95 MABN328, Merck, 1/1000; PAD2 ab50257, Abcam, 1/1000; PAD3 ab50246, 1/1000; PAD4 ab50247, 1/1000; citH3 ab5103, 1/1000; CD63 ab216130, 1/1000; Flot-1 ab41927, 1/2000). The membranes were washed in TBS-T for 3 × 10 min at RT and thereafter incubated in the corresponding secondary antibody (anti-rabbit IgG BioRad or anti-mouse IgM BioRad, diluted 1/4000 in TBS-T) for 1 h at RT. Membranes were washed for 6 × 10 min in TBS-T and visualisation performed using electrochemiluminescence (ECL) and the UVP BioDoc-ITTM System (Thermo Fisher Scientific, UK).

Immunoprecipitation and identification of deiminated proteins in llama serum and EVs

For isolation of total deiminated proteins from llama serum and serum derived EVs, the Catch and Release immunoprecipitation kit (Merck, UK) was used together with the F95 pan-deimination antibody (MABN328, Merck), which has been developed against a deca-citrullinated peptide and specifically detects proteins modified by citrullination (Nicholas and Whitaker, 2002). For F95 enrichment, 50 μl serum was used according to the manufacturer’s instructions (Merck). For EVs, total protein was first extracted from EV-enriched pellets derived from 100 μl serum, using 100 μl radioimmunoprecipitation assay (RIPA) buffer, containing protease inhibitor cocktail (P8340, Sigma, UK), and shaken gently on ice for 2 h. Thereafter proteins were isolated from the EVs by centrifugation at 16,000 g for 30 min, collecting the supernatant containing the proteins. Immunoprecipitation was carried out according to the manufacturer’s instructions (Merck), using a rotating platform overnight at 4 °C. The F95 bound proteins were eluted using denaturing elution buffer, according to the manufacturer’s instructions (Merck). The F95 enriched eluates were then either analysed by Western blotting or by LC—MS/MS (Cambridge Proteomics, Cambridge, UK). For LC—MS/MS, the F95-enriched eluates were run 1 cm into a SDS-PAGE gel and the whole F95-enriched eluate was cut out as one band, whereafter it was processed for proteomic analysis (carried out by Cambridge Proteomics). Peak files were submitted to in-house Mascot (Matrix Science; Cambridge Proteomics). Databases used for protein identification (in house, Cambridge Proteomics) were as follows: Camelidae_family_20190613 (21429 sequences; 9086806 residues) and also specifically for llama: Lama_glama_20190613 (234 sequences; 52757 residues).

Results

EV analysis in llama serum

EVs from llama serum were characterised, following step-wise ultracentrifugation, by size exclusion using NTA, by morphological analysis using TEM and by Western blotting using EV-specific markers (Fig. 1 ). A poly-dispersed population of EVs in the size range of 30–576 nm, with main peaks at 38, 119, 167, 237, 323 and 403 nm was identified by NTA analysis (Fig. 1A). Western blotting confirmed that the llama serum EVs were positive for the EV-specific markers CD63 and Flotillin-1 (Fig. 1B). TEM analysis confirmed a poly-dispersed EVs population (Fig. 1C).

Fig. 1

Extracellular vesicles (EVs) isolated from llama serum. A. Nanoparticle tracking analysis shows a poly-dispersed population of EVs in the size range of 30–576 nm, with main peaks at 38, 119, 167, 237, 323 and 403 nm. B. Llama serum EVs are positive for the EV-specific markers CD63 and Flotillin-1 (Flot-1). C. Transmission electron microscopy (TEM) imaging of EVs isolated from llama serum shows a polydispersed population; scale bar represents 100 nm.

PAD and deiminated proteins in llama serum

A cross-reaction with human PAD2, 3 and 4 isozyme-specific antibodies was observed in llama serum by Western blotting, at an approximate 70–75 kDa size range as expected for mammalian PADs (Fig. 2 A). Deiminated histone H3 was also detected in llama serum by Western blotting at the expected approximate 20 kDa size (Fig. 2A). Total deiminated proteins in llama serum-EVs were detected by Western blotting using the F95 pan-deimination antibody, revealing a range of proteins between 25–100 kDa in size (Fig. 2B). Deiminated proteins were also assessed by Western blotting after F95 enrichment, using immunoprecipitation (IP) from llama serum and serum-derived EVs (Fig. 2C). Deiminated proteins from the F95 enriched eluates were further identified by LC–MS/MS analysis. In llama serum, 103 hits for camelid proteins were identified, as listed in Table 1 and Supplementary Table 1. Overall, 43 of these deiminated protein hits were common to whole serum and serum-derived EVs, while 60 hits were specific for whole serum (Fig. 2D).

Fig. 2

Western blotting of deiminated proteins and PAD in llama serum. A. Llama PAD homologues were identified at the expected size of approximately 70–75 kDa using the anti-human PAD2, PAD3 and PAD4 isozyme specific antibodies respectively. Deiminated histone H3 (citH3), representative of neutrophil extracellular traps (NETs), was verified in llama serum. B. Total deiminated proteins were assessed by Western blotting in llama serum EVs, using the F95 pan-deimination specific antibody. C. Immunoprecipitated (IP) deiminated proteins after F95-enrichment were assessed both in serum-EVs and whole serum of llama by Western blotting. D. The Venn diagram represents deiminated proteins identified in total serum and serum-derived EVs by F95 enrichment and LC–MS/MS analysis. Overall, 43 proteins were identified in common with both samples, while 60 proteins were found deiminated in serum only and 25 proteins were identified as deiminated in EVs only.

Table 1

Deiminated proteins identified by F95 enrichment and LC–MS/MS in total serum of llama (Lama glama). Deiminated proteins were isolated by immunoprecipitation using the pan-deimination F95 antibody. The F95 enriched eluate was analysed by LC–MS/MS and peak list files were submitted to in-house Mascot (Cambridge Proteomics). Peptide sequence hits scoring with Lama glama (LAMGL) are included as well as hits with other camelids (CAMFR=Camelus ferus; CAMDR=Camelus dromedaries; LAMGU=Lama guanicoe; VICPA=Vicugna pacos (Alpaca)). Hits with uncharacterized proteins are omitted in the list. For a full list of peptide sequences and m/z values see Supplementary Table 1. An asterix (*) indicates that the protein hit is specific to whole serum only.

Protein name (*unique for serum)	Number of peptide sequences identified	Total score (p<0.05)ⱡ
*O97643_LAMGL	6	282
Fibrinogen A-alpha chain
*P01973|HBA_LAMGL	5	189
Hemoglobin subunit alpha
P68226|HBB_LAMGL	4	166
Hemoglobin subunit beta
Q865W8_LAMGL	4	109
Beta actin
A0A1W5VKM5_LAMGL	2	117
Anti-RON nanobody
*A0A1W5VKM7_LAMGL	2	111
Anti-RON nanobody
*A0A1W5VKR8_LAMGL	2	96
Anti-RON nanobody
*A0A1W5VKQ9_LAMGL	1	66
Anti-RON nanobody
S9XDK9_CAMFR	34	1575
Complement C3-like protein
S9WI87_CAMFR	20	972
Serum albumin
S9XAP9_CAMFR	9	522
Keratin, type I cytoskeletal 14-like protein
S9Y6J1_CAMFR	11	502
Keratin, type II cytoskeletal 5 isoform 13-like protein
*S9 × 688_CAMFR	10	462
Keratin 6A-like protein
*S9YD43_CAMFR	9	443
Complement component 4A-like protein
*S9Y253_CAMFR	8	380
Kininogen-2 isoform I
*T0MII3_CAMFR	7	371
Alpha-2-macroglobulin-like protein
S9XI90_CAMFR	8	357
Keratin, type II cytoskeletal 75-like isoform
S9 × 494_CAMFR	6	285
Keratin, type I cytoskeletal 42
S9XBS9_CAMFR	6	263
Ig gamma-3 chain C region
*S9XYY2_CAMFR	5	246
Hemopexin
A0A075T9L1_CAMDR	4	174
Dipeptidylpeptidase 4
S9XXW2_CAMFR	3	172
Fibrinogen beta chain
*S9WDV3_CAMFR	3	169
Fibrinogen gamma chain isoform gamma-B
*A0A1K0GY87_VICPA	4	166
Globin A1
S9WB99_CAMFR	3	146
Histone H2B
S9XNF8_CAMFR	3	142
Xaa-Pro dipeptidase
*S9YS49_CAMFR	3	132
Putative E3 ubiquitin-protein ligase Roquin
*S9YGW7_CAMFR	3	123
Heparin cofactor 2
S9WPM4_CAMFR	3	115
Adiponectin
A2V743_CAMDR	4	109
Beta actin
*S9XP08_CAMFR	1	109
Inter-alpha-trypsin inhibitor heavy chain H1	2
S9XM15_CAMFR	2
Ferritin
T0NNK2_CAMFR	1	93
L-lactate dehydrogenase
S9 × 3E8_CAMFR	1	90
Ig kappa chain V-II region RPMI 6410-like protein
A0A0A0PAR2_CAMDR	2	88
Heat shock protein 90
*S9XYF2_CAMFR	2	87
Heat shock cognate protein HSP 90-beta-like isoform 3
S9XHZ4_CAMFR	1	83
Phosphotriesterase-related protein
S9XM68_CAMFR	2	73
Xaa-Pro dipeptidase isoform 3
S9WT57_CAMFR	2	73
Tubulin beta chain
*S9YV02_CAMFR	1	69
Non-specific protein-tyrosine kinase
*S9XC57_CAMFR	1	63
Plasminogen
*S9YL21_CAMFR	2	62
Apolipoprotein A-I
*S9WKZ8_CAMFR	1	62
Inter-alpha-trypsin inhibitor heavy chain H4
S9WIA5_CAMFR	1	61
Glutathione synthetase
T0MHN9_CAMFR	1	60
Pyruvate kinase
*S9Y4U4_CAMFR	1	58
Complement C1q subcomponent subunit C isoform 2
S9WAX5_CAMFR	2	57
Unconventional myosin-Va isoform 2
S9WVY1_CAMFR	1	56
Actinin, alpha 1 isoform 6-like protein
S9XA40_CAMFR	1	55
Heat shock cognate protein
*S9YFM0_CAMFR	1	52
Keratin, type II cytoskeletal 71
A0A0U2KTX5_CAMDR	1	52
VHH5 (Fragment)
S9XR87_CAMFR	1	50
Ig lambda chain C regions isoform 19-like protein
S9WGH8_CAMFR	1	49
Lysozyme
S9WMX2_CAMFR	2	49
Dystonin
*S9 × 8K9_CAMFR	1	46
Transaldolase
*T0NM23_CAMFR	2	45
Rootletin
S9W6I0_CAMFR	1	45
Ferritin
S9WF34_CAMFR	1	43
Tubulin alpha chain
*S9YSI7_CAMFR	1	43
Triosephosphate isomerase
*S9XSQ6_CAMFR	1	40
Vitamin d-binding protein-like protein
*S9YMC0_CAMFR	2	40
Transcription factor 20 isoform 1
*S9Y636_CAMFR	1	39
Receptor-type tyrosine-protein phosphatase-like N
*S9 × 6M4_CAMFR	1	37
Dyslexia-associated protein
T0MH94_CAMFR	1	36
Rabenosyn-5-like protein
S9WJW3_CAMFR	1	35
N6-adenosine-methyltransferase subunit
*A8IY99_LAMGU	1	35
Gamma-fibrinogen
S9WRI7_CAMFR	1	35
Nuclear receptor coactivator 5 isoform 3-like protein
S9W421_CAMFR	1	35
Hemoglobin, epsilon 1
*S9W711_CAMFR	1	34
Charged multivesicular body protein 4c
*S9 × 089_CAMFR	1	34
Ig lambda chain V-III region LOI-like protein
S9WUC8_CAMFR	1	32
Ig kappa chain V-II region RPMI 6410-like protein
*S9WVI6_CAMFR	1	32
Complement C1q subcomponent subunit A
*S9Y5S1_CAMFR	1	32
Transcriptional repressor NF-X1
*S9YC53_CAMFR	1	32
Alpha-1-antitrypsin-like protein
*S9Y3F6_CAMFR	1	31
Dual specificity testis-specific protein kinase 1
S9WKI8_CAMFR	1	31
HEAT repeat-containing protein 7B1
*S9WVS9_CAMFR	1	29
Peroxisome proliferator-activated receptor gamma coactivator-related protein 1
*S9XVK5_CAMFR	1	29
Transthyretin
*S9Y967_CAMFR	1	28
General transcription factor II, i isoform 4 isoform 1-like protein
*T0MC04_CAMFR	1	28
Spermatogenesis-associated protein 2-like protein
*S9YSZ6_CAMFR	1	28
Centromere protein J

Ions score is -10*Log(P), where P is the probability that the observed match is a random event. Individual ions scores > 20 indicated identity or extensive homology (p < 0.05). Protein scores were derived from ions scores as a non-probabilistic basis for ranking protein hits. Cut-off was set at Ions score 20.

Identification of deiminated proteins in EVs from llama serum

Llama serum-derived EVs showed positive for deiminated proteins by Western blotting, using the pan-deimination F95 antibody (Fig. 2B). Deiminated proteins were further identified by F95 enrichment and LC–MS/MS analysis, revealing 68 deiminated protein hits in total for EVs, with 25 hits unique to EVs (not identified from serum). Peptide sequences of hits with camelid proteins and m/z values are listed in Table 2 and Supplementary Table 2. Overlap with deiminated proteins identified in whole llama serum and EVs is represented in the Venn diagram in Fig. 2D.

Table 2

Deiminated proteins identified by F95 enrichment and LC–MS/MS in EVs isolated from serum of llama (Lama glama). Deiminated proteins were isolated by immunoprecipitation using the pan-deimination F95 antibody, the F95 enriched eluate was analysed by LC–MS/MS and peak list files were submitted to Mascot. Peptide sequence hits scoring with L. Lama glama (LAMGL) are presented as well as hits with other camelids (CAMFR=Camelus ferus; CAMDR=Camelus dromedaries; LAMGU=Lama guanicoe). Hits with uncharacterised proteins are not listed. For a full list of peptide sequences and m/z values see Supplementary Table 2. An asterix (*) indicates that the protein hit is unique for EVs only.

Protein name (*unique for EVs)	Number of peptide sequences identified	Total score (p<0.05)ⱡ
A0A1W5VKM5_LAMGL	2	164
Anti-RON nanobody
Q865W8_LAMGL	2	85
Beta actin
*S9XAP9_CAMFR	9	554
Keratin, type I cytoskeletal 14-like protein
*S9 × 688_CAMFR	10	496
Keratin 6A-like protein
S9Y6J1_CAMFR	10	438
Keratin, type II cytoskeletal 5 isoform 13-like protein
S9WI87_CAMFR	9	430
Serum albumin
*S9YN99_CAMFR	6	417
Keratin, type I cytoskeletal 17-like isoform
*S9XI90_CAMFR	6	318
Keratin, type II cytoskeletal 75-like isoform
S9 × 494_CAMFR	5	269
Keratin, type I cytoskeletal 42
S9XBS9_CAMFR	4	162
Ig gamma-3 chain C region
A0A075T9L1_CAMDR	3	153
Dipeptidylpeptidase 4
S9 × 684_CAMFR	2	136
Keratin, type II cytoskeletal 8
S9WB99_CAMFR	3	133
Histone H2B
S9YQ51_CAMFR	3	114
Tubulin beta chain
*S9WX81_CAMFR	3	89
Histone 1, H2ai isoform 3-like protein
*S9 × 8G9_CAMFR	2	88
Desmoplakin
A2V743_CAMDR	2	85
Beta actin
A0A0A0PAR2_CAMDR	2	85
Heat shock protein 90
S9XA40_CAMFR	1	85
Heat shock cognate protein
T0NNK2_CAMFR	1	81
L-lactate dehydrogenase
S9XNF8_CAMFR	2	74
Xaa-Pro dipeptidase
T0MHN9_CAMFR	1	71
Pyruvate kinase
S9WVY1_CAMFR	1	60
Actinin, alpha 1 isoform 6-like protein
*A0A0E3Z5I3_CAMDR	1	59
Superoxide dismutase
S9XHZ4_CAMFR	1	57
Phosphotriesterase-related protein
S9W9Y4_CAMFR	1	57
Ferritin
S9XR87_CAMFR	1	56
Ig lambda chain C regions isoform 19-like protein
S9 × 3E8_CAMFR	1	50
Ig kappa chain V-II region RPMI 6410-like protein
S9WAX5_CAMFR	2	49
Unconventional myosin-Va isoform 2
S9WF34_CAMFR	1	44
Tubulin alpha chain
*S9W806_CAMFR	2	43
Filamin-A isoform 1
*S9 × 6X3_CAMFR	2	42
Scaffold attachment factor B-like protein
T0MH94_CAMFR	1	40
Rabenosyn-5-like protein
*S9Y0S0_CAMFR	1	39
DNA-directed RNA polymerase subunit beta
S9WJW3_CAMFR	1	38
N6-adenosine-methyltransferase subunit
S9XM68_CAMFR	1	38
Xaa-Pro dipeptidase isoform 3
*S9YV02_CAMFR	1	38
Non-specific protein-tyrosine kinase
S9WGH8_CAMFR	1	32
Lysozyme
S9WRI7_CAMFR	1	32
Nuclear receptor coactivator 5 isoform 3-like protein
S9W421_CAMFR	1	32
Hemoglobin, epsilon 1
*S9WI71_CAMFR	1	32
Metabotropic glutamate receptor 3
*S9WB50_CAMFR	1	31
TSC22 domain family protein 3-like protein
S9WMX2_CAMFR	1	31
Dystonin
S9WKI8_CAMFR	1	31
HEAT repeat-containing protein 7B1
S9WIA5_CAMFR	1	31
Glutathione synthetase
*S9XC05_CAMFR	1	31
Telomere-associated protein RIF1 isoform 1
A0A0U2KTX5_CAMDR	1	30
VHH5
*T0MGG7_CAMFR	1	30
Nucleoredoxin
*S9XET3_CAMFR	1	30
Rac GTPase-activating protein 1
S9XDK9_CAMFR	1	30
Complement C3-like protein
*S9XMI2_CAMFR	1	29
Pseudopodium-enriched atypical kinase 1
*S9Y3S9_CAMFR	1	29
Core histone macro-H2A.1 isoform 2
S9XXW2_CAMFR	1	29
Fibrinogen beta chain
*S9YGX6_CAMFR	1	29
PAS domain-containing serine/threonine-protein kinase
*T0NMU1_CAMFR	1	28
SH2 domain-containing protein 7
*T0MIT6_CAMFR	1	28
Serine-tRNA ligase, mitochondrial
*S9W449_CAMFR	1	28
Fc receptor-like protein 5

Ions score is -10*Log(P), where P is the probability that the observed match is a random event. Individual ions scores > 22 indicated identity or extensive homology (p < 0.05). Protein scores were derived from ions scores as a non-probabilistic basis for ranking protein hits. Cut-off was set at Ions score 20.

Discussion

For the first time deiminated proteins are described in a camelid, using the llama (Lama glama) as a model species. Post-translational deimination was identified in key immune, nuclear and metabolic proteins. A llama PAD homologue was identified at an expected 70–75 kDa size similar as for human PADs by Western blotting via cross reaction with anti-human PAD2, PAD3 and PAD4 antibodies. PAD2 is known to be the phylogenetically most conserved PAD form (Vossenaar et al., 2003; Magnadóttir et al., 2018a) and has also been seen in shark (Criscitiello et al., 2019). Deiminated histone H3, a marker of neutrophil extracellular trap formation (NETosis), was also detected in llama serum by Western blotting and is described for the first time in a camelid species but was recently described in shark (Criscitiello et al., 2019). NETosis is driven by PADs (Li et al., 2010), is conserved throughout phylogeny and is important in innate immune defences against a range of pathogens including bacteria, viruses and helminths (Brinkmann et al., 2004; Branzk et al., 2014; Schönrich and Raftery, 2016). NETosis has also been associated with clearance of apoptotic cells and during tissue remodelling in teleosts (Magnadóttir et al., 2018a, 2019a). Furthermore, NETosis is linked to a range of autoimmune diseases, due to NET activation via neo-epitopes (O’Neil and Kaplan, 2019), which can also lead to organ damage (Lee et al., 2017). NETosis is also associated with cancer (Gonzalez-Aparicio and Alfaro, 2019) and neurodegenerative diseases (Pietronigro et al., 2017).

Further deiminated proteins identified in llama serum and serum-derived EVs by F95 enrichment and LC—MS/MS analysis included key proteins of camelid innate and adaptive immunity, nuclear proteins, as well as proteins involved in metabolic function. Using STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) analysis (https://string-db.org/) for protein-protein interaction, the PPI enrichment value for deiminated proteins in whole llama serum was found to be p <  1.0e-16 for 43 deiminated proteins out of 103 identified in serum, which indicates that these proteins have more interactions among themselves than what would be expected for a random set of proteins of similar size, drawn from the genome (STRING analysis, see Fig. 3 A). For deiminated proteins in 29 out of the 68 deiminated EV cargo proteins, the PPI enrichment value was found to be p = 0.0193 (STRING analysis, see Fig. 4 A). Such enrichment indicates that the proteins are at least partially biologically connected, as a group. As the camelid proteins are not present in the STRING database, corresponding protein homologues for human were chosen to create the protein-protein interaction networks shown in Fig. 3, Fig. 4. Out of the 103 camelid proteins identified as deiminated in serum, protein IDs for 43 were present for Homo sapiens to perform the assessment of protein-protein interactions and identification of main biological GO pathways (response to stress, response to wounding, oxygen transport, small molecule metabolic process, vesicle mediated transport and regulated exocytosis; Fig. 3B). For llama EVs, out of 68 proteins identified as deiminated, 29 corresponding protein homologues were found in the STRING database for Homo sapiens and were analysed highlighting main biological GO pathways (response to stress, cytoskeleton organisation and vesicle-mediated transport; Fig. 4B).

Fig. 3

Protein-protein interaction networks of deiminated protein hits identified in whole llama (Lama glama) serum. Reconstruction of protein-protein interactions based on known and predicted interactions using STRING analysis. Due to annotations for camelids not being present in STRING, proteins are based on corresponding human protein identifiers. A. Coloured nodes represent query proteins and first shell of interactors; white nodes are second shell of interactors. B. Biological GO processes are highlighted for the same protein network as follows: red = response to stress; blue = response to wounding; green = vesicle mediated transport; yellow = oxygen transport; purple = regulated exocytosis; dark green = small molecule metabolic process. Coloured lines indicate whether protein interactions are identified via known interactions (curated databases, experimentally determined), predicted interactions (gene neighbourhood, gene fusion, gene co-occurrence) or via text mining, co-expression or protein homology (see colour key for connective lines in A).

Fig. 4

Protein-protein interaction networks of deiminated protein hits identified in EVs of llama (Lama glama) serum. Reconstruction of protein-protein interactions based on known and predicted interactions using STRING analysis. Due to annotations for camelids not being present in STRING, proteins are based on corresponding human protein identifiers. A. Coloured nodes represent query proteins and first shell of interactors. B. Biological GO processes are highlighted as follows: red = response to stress; blue = cytoskeleton organisation; green = vesicle mediated transport. Coloured lines indicate whether protein interactions are identified via known interactions (curated databases, experimentally determined), predicted interactions (gene neighbourhood, gene fusion, gene co-occurrence) or via text mining, co-expression or protein homology (see colour key for connective lines in A).

Deimination protein candidates identified here in llama serum and EVs, which are involved in immune, nuclear and metabolic functions, are further discussed below, including where appropriate in a comparative context with relevant human diseases. Proteins that have previously been identified as deiminated in other species are listed first.

are based on immunoglobulin single variable domains, derived from the variable domains of heavy chain-only antibodies, which occur naturally in camelids (Deschaght et al., 2017). In the llama the heavy chain-only antibodies are comprised of at least two subclasses (Henry et al., 2019). These types of homodimeric, light chain-less antibodies have evolved through convergent evolution (Brooks et al., 2018), in at least two groups (camelids and cartilaginous fish), and their variable binding domains (VHHS) are of great value for therapeutic and diagnostic applications (Muyldermans, 2013; Criscitiello, 2014; Steeland et al., 2016; Könning et al., 2017; Henry et al., 2019). Nanobodies can act against challenging targets such as small molecules and toxins (Wesolowski et al., 2009; Bever et al., 2016), viruses (Wei et al., 2011; Hassiki et al., 2016; Vanlandschoot et al., 2011; Cohen, 2018), enzymes (Muyldermans, 2013), ion channels (Wei et al., 2011; Danquah et al., 2016) and G protein-coupled receptors GPCRs (Cromie et al., 2015). Llama nanobodies have been shown to tether to early endosomes and to mitochondria (Traub, 2019), be used for diagnostics (Shriver-Lake et al., 2018), be used for design of cancer immunotherapeutics (Hussack et al., 2018; Bannas and Koch-Nolte, 2018; Rossotti et al., 2019) and have been approved for passive immunotherapy (Sheridan, 2019). Our current finding, that llama nanobodies can be post-translationally deiminated may shed some light on their observed structural variation which still remains to be fully explained as sequence alignment does not fully elucidate their diversity (Mitchell and Colwell, 2018). As a structurally analogous immunoglobulin in shark, new antigen receptor (NAR) (Greenberg et al., 1995; Barelle et al., 2009; Flajnik and Dooley, 2009; De Silva et al., 2019) was recently also found to be deiminated (Criscitiello et al., 2019), our current finding may provide novel insights into function of these immune proteins and be useful for refinement in therapeutic nanobody development. Llama single-chain antibodies were here found to be deimination candidates both in llama whole serum and in serum derived EVs, highlighting also their EV-mediated export.

were identified here as being deiminated in llama serum and in serum-derived EVs, scoring with Ig components from other camelids. Immunoglobulins (Ig) are key molecules in adaptive immunity but post-translational deimination of Ig’s has hitherto received little attention. Deimination of the IgG Fc region in patients with bronchiectasis and RA has been identified (Hutchinson et al., 2017). Furthermore, deimination of Ig’s in teleost fish was recently described (Magnadóttir et al., 2019a), as well as in shark (Criscitiello et al., 2019). Given the increased focus on understanding of Ig diversity throughout phylogeny (Smith et al., 2012; Zhang et al., 2013; de los Rios et al., 2015; Zhang et al., 2016, 2017; Stanfield et al., 2018) and the unique features of camelid immunoglobulins (Plasil et al., 2019), our current finding of deimination of llama Ig’s highlights a novel concept that may further understanding of Ig diversity throughout evolution. Of additional interest is the finding that some deiminated Ig proteins were also found to be exported in EVs.

identified to be deiminated in llama serum included C3, C4 and C1q, while only C3 was identified to be deiminated in serum derived EVs. Complement component C3 plays a central role in all pathways of complement activation and can also be directly activated by self- and non-self surfaces via the alternative pathway without a recognition molecule (Dodds and Law, 1998; Dodds, 2002). In camelids, the alternative and classical pathway haemolytic activity of serum has been assessed in the dromedary camel with respect to age and gender (Olaho-Mukani et al., 1995a; and 1995b). Hitherto little is known about roles for post-translational deimination of complement components throughout phylogeny. In teleost fish, C3 has been identified in serum in deiminated form (Magnadóttir et al., 2019a) and also found to be deiminated in mucosal and serum EVs of teleost fish (Magnadóttir et al., 2019b, 2019c), while post-translationally deiminated C3 was identified in shark total serum but not EVs (Criscitiello et al., 2019). Other complement components, including C4 and C1q, which belong to the classical pathway of complement activation were here identified as deiminated in llama whole serum and some of those have also recently been reported to be deiminated in teleost fish (Magnadóttir et al., 2019a). The C1q subcomponent can bind to the Fc region of immunoglobulins that are bound to antigen and activate the classical part of the complement pathway (Reid et al., 2002; Reid, 2018). Interestingly, an essential role for arginine in C1q has previously been suggested for C1q-IgG interaction (Kojouharova et al., 2004). C1q also serves as a potent pattern recognition molecule which recognises self, non-self and altered self-signals (Nayak et al., 2012; Reid, 2018) and may therefore also bind to deiminated neo-epitopes (Magnadóttir et al., 2019a). The complement system has multifaceted roles. It forms part of the first line of immune defence against invading pathogens, acting in clearance of necrotic or apoptotic cells (Dodds and Law, 1998; Sunyer and Lambris, 1998; Fishelson et al., 2001; Carroll and Sim, 2011). Complement also has roles in regeneration (Del Rio-Tsonis et al., 1998; Haynes et al., 2013) and tissue remodelling (Lange et al., 2004a, 2004b; Lange et al., 2005, 2006, 2019). Furthermore, C1 is also implicated in multiple non-complement functions including binding of apoptotic cells, cleavage of nuclear antigens and cleavage of MHC class I (Lu and Kishore, 2017). Post-translational deimination of complement components may possibly influence their function including cleavage ability, binding, deposition and generation of the convertase.

is primarily involved in lipid metabolism where conformational plasticity and flexibility are regarded as key structural features (Arciello et al., 2016). Apo A–I is associated with regulation of mitochondrial function and bioenergetics (White et al., 2017). Furthermore, Apo A–I has been shown to have a regulatory role in the complement system by affecting membrane attach complex (MAC) assembly and thus the final lytic pathway (Hamilton et al., 1993; Jenne et al., 1991; French et al., 1994). Given the diverse roles of Apo A–I, the current finding of deiminated forms in serum may be of quite some relevance and has previously been identified in teleost fish (Magnadóttir et al., 2019a). Apo A–I was here found to be deiminated in whole llama serum only.

is a major acidic plasma protein in vertebrates and serves as a transport molecule for fatty acids, bilirubin, steroids, amino acids and copper, as well as having roles in maintaining the colloid osmotic pressure of blood (Peters, 1996). In camelids, total albumin levels have been assessed, for example in relation to reproductive efficiency (El-Malky et al., 2018) and as biomarkers of oxidative stress (El-Deeb and Buczinski, 2015). While albumin has been identified as a glycoprotein in some species (Metcalf et al., 2007) investigation of post-translational deimination has been limited, but was recently identified in teleost fish (Magnadóttir et al., 2019a). Serum albumin was identified as deiminated in both whole llama serum and EVs.

is a scavenger protein of haemoglobin and a predominant heme binding protein, which contributes to heme homeostasis (Smith and McCulloh, 2015; Immenschuh et al., 2017). Hemopexin also associates with high density lipoproteins (HDL), influencing their inflammatory properties (Mehta and Reddy, 2015). Hemopexin is a plasma glycoprotein that has been previously identified as a deimination candidate in teleost fish (Magnadóttir et al., 2019a) and shark (Criscitiello et al., 2019). Here, hemopexin was found deiminated in whole llama serum only.

(heavy chain H1) belongs to the serpin family of proteins, which have protease-inhibitory functions and are involved in diverse physiological and pathological processes including fertilisation, ovulation, coagulation, inflammation, as well as tumorigenesis, metastasis and dementia (Zhuo and Kimata, 2008; Weidle et al., 2018). Inter-alpha-trypsin inhibitor is synthesised in the liver, circulates in the blood and has two chains, a light and heavy chain, whereof the heavy-chain (ITIH) includes a von Willebrand domain and can interact with the extracellular matrix (Bost et al., 1998). ITIH is downregulated in tumours via methylation and ITIH2 is strongly reduced in invasive cancers (Hamm et al., 2008). While ITIH was recently identified as a deimination candidate in teleost fish (Magnadóttir et al., 2019a) further studies on the regulation of ITIH2 via post-translational deimination have not been carried out. ITIH was here identified as deiminated in whole llama serum only.

is a glycoprotein, synthesised in liver (Tennent et al., 2007) and forms part of the acute phase response as part of the coagulation cascade (Tiscia and Margaglione, 2018). In camelids, fibrinogen is a biomarker for stress and infection (Greunz et al., 2018; El-Bahr and El-Deeb, 2016; El-Deeb and Buczinski, 2015). Impaired mechanism of fibrinogen formation and fibrin polymerization are implicated with various pathologies including coagulopathies and ischemic stroke (Weisel and Litvinov, 2013), while acquired fibrinogen disorders can be associated with cancer, liver disease or post-translational modifications (Besser and MacDonald, 2016). Fibrinogen is indeed a known deimination candidate and this post-translational modification contributes to its antigenicity in autoimmune diseases (Hida et al., 2004; Muller and Radic, 2015; Blachère et al., 2017). Fibrinogen was here identified as deiminated both in whole llama serum and EVs.

participates in cytoskeletal rearrangement and its deimination has previously been linked to EV release (Kholia et al., 2015). Deimination of tubulin may therefore be crucial for facilitating EV-mediated processes in homeostasis, immune responses and in relation to pathologies. Here, tubulin was identified as deiminated in both whole llama serum and EVs.

was identified as being deiminated in both llama serum and EVs and in addition, Histone 1 (H2ai isoform 3-like protein) and core histone macro-H2A.1 isoform 2 were identified as deiminated in EVs only. Histones undergo various posttranslational modifications that affect gene regulation and can also act in concert (Latham and Dent, 2007; Bird, 2007). In addition to acetylation, phosphorylation and ubiquitination, histones are known to undergo deimination, including H2B (Sohn et al., 2015) and H2A (Hagiwara et al., 2005), as identified in this study in llama. Other histones that are known to undergo deimination include H3 and H4 (Chen et al., 2014; Kosgodage et al., 2018).

(Hsp90) was here found to be deiminated both in llama whole serum and EVs. HSP90 has been described in camelid (Saeed et al., 2015). Hsp90 is a phylogenetically highly conserved chaperone protein involved in protein folding, stabilisation of proteins against heat stress, and aids in protein degradation (Buchner, 1999; Picard, 2002). Hsp90 also stabilizes a number of proteins required for tumour growth and is therefore important in anti-cancer drug investigations (Goetz et al., 2003). Hsp90 is responsible for most of the ATPase activity of the proteasome (Imai et al., 2003) and has an ATP binding region, which also is the main binding site of drugs (Chiosis et al., 2006). In camelids, Hsp90 has been related to adaptive tolerance of camel somatic cells to acute and chronic heat shock (up to 20 h at 45 degrees Celsius), which is lethal to many mammalian cells (Saadeldin et al., 2018b). Hsp90 has previously been described to be post-translationally deiminated in rheumatoid arthritis, allowing deimination-induced shifts in protein structure to generate cryptic epitopes capable of bypassing B cell tolerance (Travers et al., 2016). It is of some interest to find that Hsp90 is also found deiminated in llama serum, further highlighting protein moonlighting functions of Hsp90 in physiological and pathophysiological context via post-translational deimination. In addition, finding post-translational deimination of the same protein throughout phylogeny also supports translational value between species to further understanding of this post-translational modification in human pathologies.

The following deimination protein candidates identified in the current study in llama serum and serum-EVs have to our knowledge not been previously reported as deiminated:

is closely related to other thioester containing proteins, such as complement proteins C3, C4 and C5 (Sottrup-Jensen et al., 1985; Davies and Sim, 1981). Alpha-2-M is phylogenetically conserved from arthropods to mammals and found at high levels in mammalian plasma. Alpha-2-M forms part of the innate immune system and clears active proteases from tissue fluids (Armstrong and Quigley, 1999). Here, Alpha-2-M was found deiminated in whole llama serum only and has not reported as deiminated before, to our knowledge.

is the most abundant secreted adipokine with pleiotrophic roles in physiological and pathophysiological processes (Fiaschi, 2019). It has received considerable interest in the field of metabolic and obesity research (Frankenberg et al., 2017; Spracklen et al., 2019), as well as in diabetes (Yamauchi et al., 2003), due to its key function in regulating glucose (Yamauchi et al., 2002). Adiponectin is furthermore linked to longevity (Chen et al., 2019a, 2019b), regenerative functions (Fiaschi et al., 2014) and has roles in myopathies, such as Duchenne muscular dystrophy and collagen VI-related myopathies (Gamberi et al., 2019). Adiponectin is also implicated in a range of cancers, often in relation to obesity (Parida et al., 2019). Furthermore, adiponectin plays roles in reproduction, embryo pre-implantation and embryonic development (Barbe et al., 2019). Due to the range of functions in relation to key pathophysiologies there is great interest in drug development for modulating adiponectin signalling (Fiaschi, 2019). Given the unique metabolic adaptive features of camelids, the identification of post-translational deimination of this key metabolic protein may be of some interest. Recent studies in rheumatoid arthritis made a correlation between inflammation, autoantibodies and adiponectin levels (Hughes-Austin et al., 2018; Liu et al., 2019), while adiponectin itself has not been previously identified to be deiminated to our knowledge. Post-translational deimination may be a hitherto unrecognized mechanism for adiponectin, also in humans, to adapt moonlighting functions via changes in protein folding and therefore interaction with other proteins. Adiponectin is a small 244 aa protein (NP_001171271.1) in humans and contains 2 unfolded regions and 7 arginine sites, while camel adiponectin has 8 arginine sites. These could be subjected to PAD-mediated deimination and therefore modulate adiponectin folding and function, depending on which arginine is deiminated. Here, deiminated adiponectin was identified in whole llama serum only.

is a plakin-family adhesion junction plaque protein and was here identified as deiminated in both llama serum and EVs. Dystonin, along with XVII collagen, form hemidesmosomes and both proteins are autoantigens believed to be responsible for the Type II hypersensitivy pathologic in the pruritic skin disease bullous pemphigoid (Bağcı et al., 2017; Basseri et al., 2018). Dystonin has also been linked to Sjögren’s syndrome and linked to a hypermethylation status (Gonzalez et al., 2011), but post-translational deimination of dystonin has not been previously described to our knowledge, although deimination is associated with a range of autoimmune diseases, including Sjögren’s syndrome (Konsta et al., 2014; Selmi et al., 2015).

also known as prolidase, is an enzyme that in humans is encoded by the PEPD gene. Post-translational modifications of prolidase regulate its enzymatic abilities. Deficiency in prolidase leads to a rare, severe autosomal recessive disorder (prolidase deficiency) that causes many chronic, debilitating health conditions in humans (Viglio et al., 2006). These phenotypic symptoms vary and may include skin ulcerations, mental retardation, splenomegaly, recurrent infections, photosensitivity, hyperkeratosis, and unusual facial appearance. Furthermore, prolidase activity was found to be abnormal compared to healthy levels in various medical conditions including: bipolar disorder, breast cancer, endometrial cancer, keloid scar formation, erectile dysfunction, liver disease, lung cancer, hypertension, melanoma, and chronic pancreatitis (Kitchener and Grunden, 2012). In some cancers with increased levels of prolidase activity, such as melanoma, the differential expression of prolidase and its substrate specificity for dipeptides with proline at the carboxyl end suggests the potential of prolidase in becoming a viable, selective endogenous enzyme target for proline prodrugs (Mittal et al., 2005). Serum prolidase enzyme activity is also currently being explored as a biomarker for diseases including chronic hepatitis B and liver fibrosis (Duygu et al., 2013; Şen et al., 2014; Stanfliet et al., 2015). Phosphorylation of prolidase has been shown to increase its activity while dephosphorylation leads to a decrease in enzyme activity. Post-translational deimination of prolidase has not been described before to our knowledge and may add to understanding of how this enzyme is regulated. Prolidase was here identified as deiminated in both llama serum and EVs.

(DPP4, also known as CD26) was here identified to be deiminated in both llama whole serum and EVs. DPP4 controls glucose homeostasis and has complex roles in inflammation and homeostasis, including in liver cytokine expression, while its activity in plasma has been shown to correlate with body weight and fat mass (Varin et al., 2019). Interestingly, in camel milk, DDP4 inhibitory peptides have been identified and suggested to play roles in the regulation of glycaemia in humans (Nongonierma et al., 2018). Furthermore, roles for DDP4 in cancer have been found to relate to its post-translational processing of chemokines, thereby limiting lymphocyte migration to sites of inflammation and tumours (Barreira da Silva et al., 2015). As the success of antitumour immune responses depends on the infiltration of solid tumours by effector T cells, a process which is guided by chemokines, DPP4 inhibitors have been suggested as a strategy to enhance tumour immunotherapy (Barreira da Silva et al., 2015). Furthermore, serum DDP4 activity levels in primary HIV infection were found to be significantly decreased and to correlate with inflammation and HIV-induced intestinal damage (Ploquin et al., 2018). Middle East respiratory syndrome coronavirus (MERS-CoV) has been found to utilize dipeptidyl peptidase 4 (DPP4) as an entry receptor, via glycosylated sites (Peck et al., 2017). Therefore the identification of DPP4 as a deimination candidate may be of some relevance as such post-translational modification can affect DPP4 structure and function, allowing for moonlighting functions which may vary in pathological compared to pathophysiological milieus.

, belongs to the Roquins which are a family of highly conserved RNA-binding proteins involved in ubiquination, are crucial for T-cell-dependent B-cell responses (Athanasopoulos et al., 2016) and play important roles in modulating T-cell activity (Akef and Muljo, 2018). Roquins repress constitutive decay elements containing mRNAs and play a critical role in RNA homeostasis and immunological self-tolerance (Zhang et al., 2015; Essig et al., 2018; Mino and Takeuchi, 2018). Roquin plays multifaceted roles both in the generation of a homeostatic immune response, as well as during chronic inflammation and autoimmunity (Schaefer and Klein, 2016; Lee et al., 2019). While roquin causes post-translational ubiquination of target proteins, post-translational deimination of roquin itself may modulate is function and is described here for the first time. Roquin was here identified as deiminated in whole llama serum only.

was here identified as a deimination candidate protein in llama EVs only. It has been linked to HUPRA syndrome which is a rare mitochondrial disease characterized by hyperuricemia, pulmonary hypertension, renal failure in infancy and alkalosis (Rivera et al., 2013). It has furthermore been linked to progressive spastic paresis (Linnankivi et al., 2016). Post-translational deimination of this mitochondrial protein has not been described before to our knowledge.

(Nrx) was here identified to be deiminated in EVs only. It is an oxidoreductase of the thioredoxin family of proteins with numerous functions in the redox regulation of metabolic pathways, cellular morphology, and signal transduction (Urbainsky et al., 2018). Nrx has been shown to inhibit Wnt-beta-catenin signalling (Funato et al., 2006) and is linked to Ca2+-mediated mitochondrial reactive oxygen species metabolism (Rharass et al., 2014). Nrx has furthermore been identified as an epigenetic cancer marker related to the oxidative status of human blood (Schöttker et al., 2015), but hitherto not been described as deiminated.

was here identified to be deiminated in total serum of llama. It is associated with developmental dyslexia (Levecque et al., 2009) and has been shown to have roles in cell-cell interactions and signalling and neuronal migration (Velayos-Baeza et al., 2008), as well as in axon guidance (Poon et al., 2011). Dyslexia-associated protein has previously been found to follow the classical clathrin-mediated endocytosis pathway and its surface expression seems regulated by endocytosis, indicating that the internalization and recycling of the protein may be involved in fine-tuning its role in neuronal migration (Levecque et al., 2009). It has been descried to be highly glycosylated in different mammalian cell lines (Velayos-Baeza et al., 2008), while deimination has not been described before.

was identified here as deiminated in llama EVs only. It has, like other components of the glutamatergic system, a widespread distribution outside the central nervous system (CNS) and has been linked to regulation of the brain-gut axis (Julio-Pieper et al., 2013). It has also been linked to pain (Acher and Goudet, 2015) and psychotic disorders (Matrisciano et al., 2016). The group III mGlu receptors have been described in human stomach and colon, revealing a huge potential for these receptors in the treatment of peripheral disorders, including gastrointestinal dysfunction (Julio-Pieper et al., 2013). As post-translational deimination has not been described before in this receptor, but may well affect is tertiary structure, our current finding may be of some relevance in relation to strategies for developing antagonistic probes (Wenthur et al., 2019). Such pharmacological tools originally designed for mGlu receptors in the CNS may also be directed towards new disease targets in the periphery. Ulcerative colitis and Crohn’s disease are potential targets, as irritable bowel diseases can be co-morbid with anxiety and depression (Julio-Pieper et al., 2013).

(GSS) was here identified to be deiminated in both llama whole serum and EVs. GSS is the second enzyme in the glutathione (GSH) biosynthesis pathway involved in homeostasis and cellular maintenance and also acts as a potent antioxidant (Njålsson and Norgren, 2005). In camelids, glutathione peroxidase, which also belongs to the GSH biosynthesis pathway and is a potent anti-oxidant, has been identified as a seminal plasma fertility biomarker (Waheed et al., 2015). Furthermore, camel milk has been shown to boost glutathione and total anti-oxidant capacity in sera and exudates in animal studies (Arab et al., 2017). In a diabetic mouse model, the activities of GSH, alongside glucose insulin and ROS levels, were restored after camel whey protein treatment (Sayed et al., 2017). In humans, GSS deficiency is an autosomal recessive disorder with varyingly severe clinical manifestations that include metabolic acidosis, hemolytic anemia, hyperbilirubinemia, neurological disorders and sepsis (Guney Varal et al., 2019). Deimination of GSS identified here has not been assessed before and may possibly affect the GSH biosynthesis pathway and such post-translational regulation remains to be further investigated.

, also known as ciliary rootlet coiled-coil protein (CROCC), is a protein that is required for centrosome cohesion and is therefore important in mitosis (Bahe et al., 2005; Graser et al., 2007). It was identified here as deiminated in whole llama serum only. Rootletin has been shown to be phosphorylated and to have the ability to form centriole-associated fibers, suggesting a dynamic model for centrosome cohesion based on entangling filaments (Bahe et al., 2005). Deletion of rootletin in mouse models causes photoreceptor degeneration and impaired mucociliary clearance, supporting its key function in rootlet structures (Yang et al., 2005). Post-translational deimination of rootletin, as identified here, may possibly facilitate its dynamic functions.

is is a highly conserved and ubiquitiously expressed centrosomal protein, involved in microtubule disassembly and plays a structural role in the maintenance of centrosome integrity, genome stability and normal spindle morphology during cell division (Tang et al., 2009; McIntyre, 2012). It was here identified as deiminated in whole llama serum only. Knockout mouse models targeting the gene encoding this protein have phenotypes of impaired glucose tolerance, hypoalbuminemia and increased micronuclei, indicative of genomic instability (Gerdin, 2010). In humans, it is associated with primary autosomal recessive microcephaly (Gul et al., 2006) and the microcephalic primordial dwarfism disorder Seckel syndrome (Al-Dosari et al., 2010; McIntyre et al., 2012). Deimination has not bee described before in this protein to our knowledge.

is involved in the repair of double-strand DNA breaks in response to DNA damage (Silverman et al., 2004; Escribano-Díaz et al., 2013; Drané et al., 2017) and protects teleomeres (Fontana et al., 2018). It was identified here as deiminated in EVs only. RIF1 has been associated with cancer via activation of human telomerase reverse transcriptase expression (Liu et al., 2018). It is also required for immunoglobulin class-switch recombination during the germinal centre reaction for humoral antibody immunity (Di Virgilio et al., 2013). RIF1 is involved in telomere homeostasis and was recently found to be post-translationally S-acylated, identifying a novel posttranslational modification regulating DNA repair (Fontana et al., 2019). Post-translational deimination has hitherto not been described.

was here identified as a deimination protein candidate in both llama whole serum and EVs. It acts in early endocytic membrane fusion and membrane trafficking of recycling endosomes (Naslavsky et al., 2004; Rai et al., 2017). It plays a role in the lysosomal trafficking of CTSD/cathepsin D from the Golgi to lysosomes and also promotes the recycling of transferrin directly from early endosomes to the plasma membrane (Navaroli et al., 2012). Rabenosyn-5-like protein also binds phospholipid vesicles and plays roles in regulating protein sorting and recycling to the plasma membrane (Nielsen et al., 2000; de Renzis et al., 2002; Naslavsky et al., 2004). Rabenosyn-5 was recently identified as an upregulated urinary biomarker associated with malignant upper gastrointestinal cancer (Husi et al., 2019) and found to be increased in diabetic kidney disease (Dumont et al., 2017). Rabenosyn-5 has been shown to be phosphorylated, regulating its recruitment to membranes (Macé et al., 2005). Post-translational deimination has hitherto not been described..

was here identified as deiminated in whole llama serum only. It is expressed in testis and to a lesser extent in spleen, thymus, and prostate (Graziotto et al., 1999). SPATA2 has been identified as a bridging factor that regulates TNF-alphA–Induced necroptosis and is instrumental for TNF-induced cell death (Elliott et al., 2016; Kupka et al., 2016; Schlicher et al., 2016; Wagner et al., 2016; Schlicher et al., 2017). SPATA2 ensures normal secretory function of Sertoli cells (Zhao et al., 2017) and has recently been identified as a novel predictor in ovarian cancer outcome (Wieser et al., 2019). While SPATA2 has been linked to ubiquination (Schlicher et al., 2016), post-translational deimination has not been invesigated.

D- belongs to the albumin gene family, together with human serum albumin and alpha-fetoprotein. It is a multifaceted protein mainly produced in the liver, where its regulation is influenced by estrogen, glucocorticoids and inflammatory cytokines (Bikle and Schwartz, 2019). It is secreted into the blood circulation and is able to bind the various forms of vitamin D including ergocalciferol (vitamin D2) and cholecalciferol (vitamin D3), the 25-hydroxylated forms (calcifediol), and the active hormonal product 1,25-dihydroxyvitamin D (calcitriol) (Verboven et al., 2002; Norman, 2008). It transports vitamin D metabolites between skin, liver and kidney, and then on to various target tissues (Norman, 2008). VDBP is a macrophage activating factor (MAF) and has been tested as an anti-cancer agent via activation of macrophages against cancer cells (Yamamoto et al., 2008). Some association has also been made between polymorphisms of VDBP and the risk of coronary artery disease (Tarighi et al., 2017). Post-translational modifications (which still remain to be identified) of VDBP have been associated with multiple sclerosis (MS) (Perga et al., 2015), while protein deimination is well known to be associated with MS (Moscarello et al., 2013). VDBP has previously been identified to be glycosylated (Kilpatrick and Phinney, 2017) but was here identified as a deimination candidate in whole llama serum only. Post-translational deimination may contribute to various functions of VDBP in physiological as well as pathophysiological processes.

(PEAK1) was here identified as deiminated in EVs only. It is a tyrosine kinase and scaffold protein that transmits integrin-mediated extracellular matrix (ECM) signals to facilitate cell movement and growth. While aberrant expression of PEAK1 has been linked to cancer progression, its normal physiological role in vertebrate biology is still relatively unknown. Deletion of the PEAK1 gene in zebrafish, mice and human endothelial cells (ECs) induced severe defects in new blood vessel formation due to deficiencies in EC proliferation, survival and migration. PEAK1 seems therefore to play roles in modulation of cell adhesion and growth factor cues from the extracellular environment necessary for new vessel formation during vertebrate development and cancer (Wang et al., 2018). PEAK1 has not been described as deiminated before.

was here identified to be deiminated in llama EVs only. Desmoplakin is a unique and critical component of desmosomal cell-cell junctions and involved in integrity of the cytoskeletal intermediate filament network (Bendrick et al., 2019). It has been shown to be required for epidermal integrity and morphogenesis in the Xenopus laevis embryo (Bharathan and Dickinson, 2019) and novel roles in coordination of cell migration were recently established (Bendrick et al., 2019). Mutations in desmoplakin have been linked to multiple allergies, severe dermatitis and metabolic wasting (SAM) syndrome (Liang et al., 2019). It is also linked to Carvajal syndrome, involving altered skin and hair abnormalities, and heart diseases (Yermakovich et al., 2018), including cardiomyopathies (Reichl et al., 2018; Chen et al., 2019a, 2019b). Desmosomal proteins have been shown to have both tumour-promoting and tumour-suppressive functions, depending of cancer types and can regulate cell proliferation, differentiation, migration, apoptosis, and impact treatment sensitivity in different types of cancers (Zhou et al., 2017). As the roles of desmosomal proteins in cancer and metastasis are not fully understood, the identification of deiminated desmoplakin in llama EVs here, and not described before, may be of some interest and add to understanding of diverse functional ability via such post-translational modification.

, also called glucocorticoid-induced leucine zipper protein (GILZ), was here identified as a deimination candidate in llama EVs only. It is a glucocorticoid-responsive molecule involved in immune regulation and glucocorticoid actions. Its interactions with signal transduction pathways, many of which are operative in RA and other inflammatory diseases, suggest that it is a key endogenous regulator of the immune response including a key glucocorticoid-induced regulator of inflammation in rheumatoid arthritis (RA) (Beaulieu and Morand, 2011). GILTZ is a small, 135-amino acid protein with anti-inflammatory properties and has been shown to inhibit NF-κB and MAPK pathways (Bereshchenko et al., 2019; Ricci et al., 2019). It has also been shown to be induced in response to hypoxia by a HIF1α-dependent mechanism and in response to cholesterol starvation, leading to downstream shedding of procoagulant EVs in ovarian cancer (Koizume et al., 2019). Post-translational deimination of GILZ has not been described before but may indeed affect its multifaceted functions, including in inflammation and cancer.

In the current study we report deiminated proteins in llama serum and serum-derived EVs. Due to the fact that the llama genome is not yet fully annotated, the hits identified here may underestimate the amount of deiminated proteins present in llama serum and EVs. Therefore a wider protein-hit analysis was carried out including other members of the camelid family, using known sequences from Camelus ferus, Camelus dromedaries, Lama guanicoe and Vicugna pacos. Deimination of key immune factors of innate and adaptive immunity and key metabolic proteins is identified here for the first time in a camelid species, highlighting putative protein moonlighting functions via post-translational deimination. It must be noted that in relation to previously observed increases of deiminated proteins with age (Ding et al., 2017), the llama used in this study would be considered relatively old at 21 years of age, as typical llama lifespans are 15–25 years, with some individuals surviving 30 years or more.

Our findings presented here furthermore complement expanding research in the comparative EV research field. Previous studies on the camel urinary proteome revealed enriched proteins from EVs and relevance to stress and immune responses as well as antimicrobial activities (Alhaider et al., 2012). Research on EVs is a relatively new field in comparative immunology and to our knowledge; this is the first description of EVs in serum of a camelid species. Previous studies on EVs in camelids focussed on EVs in camel milk, which were shown to have anti-cancerous properties (Badawy et al., 2018). As PADs have been identified to play major roles in the regulation of EV release (Kholia et al., 2015; Kosgodage et al., 2017 and 2018), including in host-pathogen interactions (Gavinho et al., 2019), such EV-mediated communication may be of great relevance also for addressing diverse zoonotic diseases identified in camels (El-Alfy et al., 2019; Zhu et al., 2019).

In continuation of the current pilot study, the assessment of changes in deiminated proteins in camelid serum, and their lateral transfer via EVs, will be of great interest to assess animal health in response to infection and environmental stress. Our findings will further current understanding of the roles for EVs, PADs and posttranslational deimination throughout phylogeny and in relation to adaption to a range of, including extreme, environments. Furthermore, novel identification of post-translational deimination in key proteins of metabolism and immunity may reveal hitherto unrecognized moonlighting function of these proteins throughout phylogeny, in relation to physiological and pathological processes, as well as being translational to and informing inflammatory and metabolic diseases. PAD-mediated contribution to protein moonlighting and in EV-mediated communication in response to physiological and pathophysiological changes remains therefore a field of further studies.

Conclusion

This is the first report of deiminated proteins in serum and serum-EVs of a camelid species, using the llama as a model animal. Our findings highlight a hitherto unrecognized post-translational modification in key immune and metabolic proteins in camelids, which may be translatable to and inform a range of human metabolic and inflammatory pathologies.

CRediT authorship contribution statement

Michael F. Criscitiello: Resources, Funding acquisition, Validation, Writing - review & editing. Igor Kraev: Formal analysis, Resources, Validation, Visualization. Sigrun Lange: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing - original draft, Writing - review & editing.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.