Next Article in Journal
Potential Immunogenic Activity of Computationally Designed mRNA- and Peptide-Based Prophylactic Vaccines against MERS, SARS-CoV, and SARS-CoV-2: A Reverse Vaccinology Approach
Previous Article in Journal
Structural Properties of Phenylalanine-Based Dimers Revealed Using IR Action Spectroscopy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 27
Issue 7
10.3390/molecules27072371
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Functional Protein Composition in Femoral Glands of Sand Lizards (Lacerta agilis)

by
Alejandro Ibáñez
1,2,*,
Bozena Skupien-Rabian
3,
Urszula Jankowska
3,
Sylwia Kędracka-Krok
4,
Bartłomiej Zając
5 and
Maciej Pabijan
1
1
Department of Comparative Anatomy, Institute of Zoology and Biomedical Research, Faculty of Biology, Jagiellonian University, 30-387 Kraków, Poland
2
Department of Ecology and Vertebrate Zoology, Faculty of Biology and Environmental Protection, University of Łódź, 90-237 Łódź, Poland
3
Proteomics and Mass Spectrometry Core Facility, Malopolska Centre of Biotechnology, Jagiellonian University, 30-387 Kraków, Poland
4
Department of Physical Biochemistry, Faculty of Biochemistry, Biophysics and Biotechnology, Jagiellonian University, 30-387 Kraków, Poland
5
Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, 30-387 Kraków, Poland
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(7), 2371; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27072371
Submission received: 16 February 2022 / Revised: 29 March 2022 / Accepted: 31 March 2022 / Published: 6 April 2022
(This article belongs to the Section Chemical Biology)
Browse Figures
Versions Notes

Abstract

:
Proteins are ubiquitous macromolecules that display a vast repertoire of chemical and enzymatic functions, making them suitable candidates for chemosignals, used in intraspecific communication. Proteins are present in the skin gland secretions of vertebrates but their identity, and especially, their functions, remain largely unknown. Many lizard species possess femoral glands, i.e., epidermal organs primarily involved in the production and secretion of chemosignals, playing a pivotal role in mate choice and intrasexual communication. The lipophilic fraction of femoral glands has been well studied in lizards. In contrast, proteins have been the focus of only a handful of investigations. Here, we identify and describe inter-individual expression patterns and the functionality of proteins present in femoral glands of male sand lizards (Lacerta agilis) by applying mass spectrometry-based proteomics. Our results show that the total number of proteins varied substantially among individuals. None of the identified femoral gland proteins could be directly linked to chemical communication in lizards, although this result hinges on protein annotation in databases in which squamate semiochemicals are poorly represented. In contrast to our expectations, the proteins consistently expressed across individuals were related to the immune system, antioxidant activity and lipid metabolism as their main functions, showing that proteins in reptilian epidermal glands may have other functions besides chemical communication. Interestingly, we found expression of the Major Histocompatibility Complex (MHC) among the multiple and diverse biological processes enriched in FGs, tentatively supporting a previous hypothesis that MHC was coopted for semiochemical function in sand lizards, specifically in mate recognition. Our study shows that mass spectrometry-based proteomics are a powerful tool for characterizing and deciphering the role of proteins secreted by skin glands in non-model vertebrates.

1. Introduction

Proteins are structurally complex macromolecules with a wide variety of functions. The chemical versatility of proteins makes them ideal candidates for intraspecific chemical communication, including their putative involvement in the social and sexual interactions of animals [1,2]. For instance, darcin—a male-specific protein—mediates female attraction towards male urine scent marks in mice [3]. In Urodela—e.g., salamanders and newts—proteinaceous pheromones are expressed in cutaneous and cloacal glands [4,5,6]. Indeed, many vertebrates possess integumentary glands that secrete a wide variety of chemical substances. Yet, the vast majority of molecules, especially proteins, remain uncharacterized for non-model vertebrates.
Lizards (order Squamata) have both an olfactory epithelium and Jacobson’s (vomeronasal) organ that discriminate a range of compounds and/or scents [7,8,9]. In addition, some species of lizards possess femoral glands (FGs), epidermal structures located on the ventral epidermis of the hind limbs that secrete a rich mixture of chemicals [10,11,12]. Femoral glands are usually dimorphic and, in general, more developed in males than in females [10,12,13]. Recent comparative analyses suggested that FGs and other epidermal glands may have played a role in the evolution of sociality in this vertebrate group [14]. Indeed, it is well known that FG secretions are involved in many social interactions (reviewed in [12]). For instance, they may be used to mark territories and establish dominance relationships in territorial species [15,16,17]. Additionally, accumulating evidence indicates that femoral secretions may provide information about the signaler’s condition, facilitating the assessment of potential partners and thus ultimately influencing mate choice [12,18,19].
In the desert iguana, Dipsosaurus dorsalis, up to four fifths of the compounds in FG secretions are proteins while just around one fifth are lipids [20]. Although the lipophilic fraction is less abundant than the proteinaceous one, the former has been characterized and relatively well studied in many lizard species [12,21,22]. Studies focusing on the proteinaceous fraction, especially concerning the role of proteins, are much more scarce [20,23,24,25,26,27]. Recent research has shown that femoral gland proteins may play an important role in individual recognition in lacertid lizards, and may therefore be involved in intraspecific communication [24,26]. In addition, the proteinaceous fraction of FGs has been studied in a few lizard species, showing species-specific protein patterns, giving further support to the idea that proteins may be involved in chemical signaling [27]. Nevertheless, FG secretions, and in particular proteins, might have other unknown functions as well. A recent study did not identify any proteins directly connected to chemical communication in the Galápagos marine iguana [28]. Unexpectedly, proteins from FG secretions were linked to the innate immune system response, indicating a potential defensive role against pathogenic agents or an inhibitory role towards microbial agents that could degrade chemosignals [28].
The sand lizard Lacerta agilis Linnaeus (1758) is widely distributed across Europe, inhabiting semi-open to open areas, such as dunes, heathland, meadows, farmland and forest margins, to heavily modified areas including roadsides, railway embankments, abandoned mineral extraction sites and vegetated urban areas [29]. This species often lives at high densities, with many lizards utilizing overlapping home ranges with numerous opportunities for encounters [29]. Confrontations may lead to visual displays and fights among individuals. The FGs of male sand lizards, as in many other lacertids [12], are most prominent during and immediately after the spring mating season. Matthey (1929) [13] demonstrated that FGs are sexual secondary organs regulated by testicular hormones that boost gland development in male sand lizards. Laboratory experiments showed that female sand lizards choose male odors based on their underlying MHC class I genetic background, aiding in communication among or within the sexes of this species [30]. Although major lipophilic classes of compounds are present in FG secretions of L. agilis [31], the protein has not been studied.
Here, we characterized the protein component of FG secretions in male sand lizards using a mass spectrometry-based proteomic approach. Our main goals were to identify proteins expressed in FG secretions and explore patterns of variation and stability of protein profiles among individuals. Then, we examined the potential functions of proteins identified in lizard FGs. As a general prediction, we expected a high diversity of proteins in FG secretions, including some that are involved in intraspecific communication. We deciphered the potential functions of proteins by using bioinformatics tools and compared our results with previous research on FG proteins.

2. Results and Discussion

2.1. Protein Identification and Inter-Individual Variation

Femoral pores and yellow secretions produced by epidermal glands are easily discernible to the naked eye (see Figure 1A). Mass spectrometry-based proteomics allowed us to quantitatively profile the proteins present in the FG secretions of 15 male sand lizards. The average (±SD) number of identified protein groups—i.e., proteins grouped together according to their shared set or subset of peptides [32]—was 815 (±195) per FG sample. The variation among individuals was relatively high, ranging from 422 to 1097 identified protein groups (see Figure 1B).
Based on all identified proteins, there were 667 protein groups having valid (i.e., non “NaN”) intensity-Based Absolute Quantification (iBAQ [33]) values in at least 10 out of 15 samples (Supplementary File S1) that were selected for further inspection. These 667 protein groups have, in general, a relatively low molecular weight with a distribution similar to that of the marine iguana proteome (see Figure 2 and Tellkamp et al. [28] for a comparison). A relatively high correlation was observed across most of the samples according to iBAQ intensities for all 667 protein groups (Figure 2). These results show that, despite the variable number of proteins expressed in FGs, there is a subset of proteins more commonly expressed across individual lizards, indicative of functional stability (see Table 1 and below for proteins grouped by function). In the common wall lizard (Podarcis muralis), individuals from the same population or from the same color morphs possess similar electrophoretic banding patterns, probably due to genetic relatedness, but remarkable inter-individual variation in protein bands was also observed [24]. Although the large variability across samples in the total number of proteins in the sand lizard is in line with previous results (see Figure 1B), our findings also highlight the presence of a group of proteins that are more commonly expressed across individuals, whose functions are explored in more detail below.

2.2. Protein Functionality and Comparison with the Galápagos Marine Iguana Proteome

The 667 protein groups were ranked according to the mean value of iBAQ intensity (Figure 3A). The 30 proteins with the highest iBAQ values (Top30) were grouped according to their functions in Table 1 (see Figure 3B for a visual representation of the proportion of Top30 proteins per functional category). We compared our results with the only other FG secretion proteome available, that of the distantly related Galápagos marine iguana [28].
The most abundant protein was cathepsin D (Figure 3C), which is a soluble lysosomal aspartic endopeptidase involved in numerous physiological processes, including degradation of intracellular proteins, activation and degradation of polypeptide hormones and growth factors, activation of enzymatic precursors, and brain antigen processing [35]. Another protein from the cathepsin family, cathepsin A, a serine protease [36], was also found among the Top30 proteins. Cathepsins play a significant role in immune responses through their ability to process antigens [37,38]. Cathepsin D was also one of the most abundant proteins found in FG secretions of the marine iguana [28].
The second most abundant protein in FG secretions of male sand lizards was carbonic anhydrase. This protein catalyzes conversion of carbon dioxide to bicarbonate, which might be responsible, e.g., for pH regulation [39].
Interestingly, the group of Top30 proteins contained two key antioxidant enzymes, superoxide dismutase (SOD1) and peroxiredoxin (PRDX4). Other antioxidant enzymes were repeatedly identified in our samples, including catalase with intensity rank (IR) of 112 and glutathione peroxidase (GPX6 and GPX1, IR of 180 and 252, respectively). These proteins may potentially have a protective role against the oxidation of the other components of FG secretion; e.g., antioxidants could protect lipids from peroxidation due to environmental agents. For instance, it has been hypothesized that differences in the lipophilic components of different populations of Zootoca vivipara could be driven by environmental humidity that may lead to the production of compounds preventing lipid oxidation in FG secretions [40]. Our results provide provisional evidence that proteins secreted in FGs may shield lipid-based semiochemicals from oxidative damage, and hence increase their resistance to environmental fade-out.
The Top30 proteins also included chaperones, cytoskeletal proteins, histones and others (see Table 1 and Figure 3B). As mentioned in Tellkamp et al. [28], the presence of intracellular proteins in FG secretion is expected due to the holocrine nature of the FGs. Moreover, proteins indicating the epidermal origin of FGs were detected in our experiment, as in Tellkamp et al. [28]. Several keratins (also in the Top30 group), proteins from the S100A family (IR: 23, 133 and 289), 14-3-3 sigma protein (IR: 84) as well as Tgm1 (IR: 373) were detected. Consistent with the hydrophobic nature of FG secretions, several lipid-related proteins were recognized as some of the most abundant (see Figure 3A,B).
The dataset of repeatedly identified proteins (667 protein groups) was also analyzed in terms of function using STRING (Search Tool for the Retrieval of Interacting Genes/Proteins; see Material and Methods). This analysis showed an enrichment in several biological processes, including immune system process (GO:0002376) and lipid metabolic process (GO:0006629) (see Supplementary File S1), which is in line with the study of the Galápagos marine iguana. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment also showed numerous ribosomal and proteasomal proteins in the dataset (see Supplementary File S1). Considering protein intensity within the samples, proteasomal proteins were less abundant than the other mentioned protein groups (Figure 3A). Immune- and lipid-related proteins were distributed across the whole range of intensity, while ribosomal and epidermal proteins attained higher positions in the intensity ranking (see Figure 3A). Analysis with the DAVID tool confirmed the STRING findings and revealed strong over-representation of extracellular exosome proteins within the analyzed FG protein set (proteins of the GO TERM “extracellular exosome” are shown in Supplementary File S2; see sheet “GO SLIMMER_ LIST OF TERMS”). Exosomes are membranous vesicles secreted by cells involved in intercellular signaling processes [41]. Femoral secretions are made up of components and organelles derived from cell degradation, with small vesicles also present [42,43]. Therefore, the observed over-representation of extracellular exosome proteins in FGs is a clear indication of secretory activity.
Our approach to the annotation of the sand lizard proteome relies on the use of high-quality but taxonomically restricted, annotated databases for human proteins. Our analysis therefore assumes similar gene functions in lacertids as in humans, which is a clear simplification. However, since only about 1% of protein sequences in the UniProt database have experimentally verified functions, the transfer of protein annotations between different species is a widespread approach [44,45,46,47]. Nonetheless, our results should be taken as a first approximation of the potential functions of FG proteins in sand lizards, and are in need of validation in future studies.

2.3. Are Sand Lizard FG Proteins Involved in Intraspecific Chemical Communication?

Contrary to our expectations, no proteins obviously linked to pheromones were found in this study. Several hypotheses could explain this result: (1) proteins related to chemical signaling remain in the non-identified proteinaceous fraction; (2) proteins might have multiple functions including a possible role in pheromone communication that cannot yet be discerned or annotated due to limitations in the databases; (3) proteins in sand lizard FGs are not directly involved in the semiochemical ecology of this species. Proteins related to lipid metabolism were well represented across the spectrum of functional groups (see Figure 3A). Among Top30 proteins, we found phospholipase A2 (EC 3.1.1.4) (IR: 24) and an aldo_ket_red domain-containing protein (IR: 28). These proteins might be connected to metabolism of the lipophilic compounds present in FGs. Phospholipases are enzymes that hydrolyze phospholipids into fatty acids and other metabolites. Aldo-keto reductase family members play a pivotal role in the biosynthesis and metabolism of steroids [48]. In the case of sand lizards, steroids and fatty acids account for roughly 50% of the lipophilic fraction of femoral secretions [31]. Therefore, phospholipases and aldo-keto reductases could be involved in the metabolism or production of specific lipids in FGs, e.g., steroids and fatty acids, that could be further used in chemical signaling [11,22]. In addition, a fatty acid binding protein (FABP5; IR: 41) was found among the 667 highly expressed proteins. Fatty acid binding proteins play an important role in the transport of some lipids. Although the connection of lipid metabolic proteins with chemical communication is unknown, these, or similar proteins, could be directly or indirectly involved in the biosynthesis and/or transport of chemosignals.
Antigen processing and presentation of exogenous peptide antigen via MHC (Major Histocompatibility Complex) class II (GO:0019886), as well as antigen processing and presentation of peptide antigen via MHC class I (GO:0002474), were detected among the multiple and diverse biological processes enriched in FGs (see Supplementary File S1). Previous research has demonstrated that the MHC class I genotype is correlated with mate acquisition and reproductive success in sand lizards [49]. Moreover, females perceive olfactory substances released by males, and prefer males with a more distinct MHC class I genotype [30]. The MHC genotype is associated with sexual secondary traits and chemical signaling in vertebrates [50]. Therefore, a tentative hypothesis is that MHC proteins present in FGs are recognized and used by females to choose potential partners (e.g., to choose males with a more distinct MHC background [30]). However, we note that caution should be used when interpreting GO enrichment analysis relying on human databases, for distantly related vertebrates. Whether MHC proteins might have a pheromonal function remains unknown, but our results provide a framework to further explore the potential connection between the FG proteome, MHC and mate choice.
To determine whether proteins are directly involved in intraspecific communication, large-scale studies on FG proteomics under a macroecological framework are needed. For instance, a good approach could be to first explore the relationship between chemical signal complexity (e.g., richness of proteins) and traits that are under sexual selection, such as body size, using a large number of lizard species [51]. In addition, other factors such as climatic or environmental conditions, as well as diet, could also influence protein composition in lizards. For example, in Schreiber’s green lizard (Lacerta schreiberi), it has been shown that experimental manipulation of the natural diet may affect the lipophilic composition of FG secretions, and could ultimately influence mate choice [52]. Similar experimental studies could be used to test whether diet may modulate protein expression in FGs. Additionally, further research integrating mass-spectrometry proteomics with evolutionary comparative analyses, as well as experimental studies, are necessary to shed light on protein functionality in the context of chemical signaling.

3. Materials and Methods

3.1. Lizard Capture and Sampling

Most lizards were collected within the city of Kraków in the surroundings of an old quarry (50.034° N, 19.865° E); only male sand lizards were sampled. Two animals were collected 30 km east of this area, at the edge of a forest (50.004° N, 20.274° E). Lizards were sampled in spring, coinciding with the mating season of this species. Femoral gland secretions were obtained by gently squeezing the femoral gland pores. The exudates were stored in 100 µL of lysis buffer (5% SDS, 0.1 M Tris, pH 7.5, 1 mM PMSF). The secretions were ground in tubes with a pestle, incubated at approximately 100 °C for 5 min and refrigerated before further processing. Proteome analysis was conducted in samples of femoral secretions originating from 15 male lizards.
The samples were sonicated for 15 min at 320 W (intensity setting—high) with time interval 30 s/30 s ON/OFF using a Bioruptor UCD-200 sonicator (Diagenode, Liege, Belgium). Then, samples were incubated at 95 °C for 5 min and centrifuged (20,000× g, 10 min, RT). Supernatants were saved. Protein digestion was performed using S-Trap™ micro spin columns according to the manufacturer’s protocol. Briefly, samples were solubilized in 5% SDS, 50 mM triethylammonium bicarbonate (TEAB; pH 7.55), reduced with 50 mM DTT solution, and alkylated with the addition of iodoacetamide to a final concentration of 40 mM, then acidified with phosphoric acid. Afterwards, 6 volumes of Strap binding buffer (90% aqueous methanol, 100 mM TEAB, pH 7.1) were added and the mixture was placed on the S-Trap by centrifugation at 4000× g for 10 s. Samples were purified by washing with S-Trap binding buffer and centrifugation. Proteins were digested overnight with trypsin (Promega, Madison, WI, USA) at 20:1 protein to enzyme wt:wt ratio. Peptides were eluted with 50 mM TEAB, followed by 0.2% aqueous formic acid and 50% acetonitrile, containing 0.2% formic acid. Lastly, peptides were vacuum dried.

3.2. LC-MS/MS Analysis

The LC-MS/MS analysis was conducted with the use of a nanoHPLC system (UltiMate 3000 RSLCnano System, Thermo Fisher Scientific, Waltham, MA, USA) coupled with a Q Exactive mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). Peptides were suspended in 2% acetonitrile (ACN) with 0.05% trifluoroacetic acid (TFA) and loaded onto a C18 trap column (AcclaimPepMap100 C18, Thermo Scientific; ID 75 μm, length 20 mm, particle size 3 μm, pore size 100 Å). Then, they were separated on a C18 analytical column (AcclaimPepMapRLSC C18, Thermo Scientific; ID 75 μm, length 500 mm, particle size 2 μm, pore size 100 Å) in a 4 h gradient of acetonitrile (2–40%) in the presence of 0.05% formic acid. Eluting peptides were ionized in a PicoView nanospray source (DPV-550, New Objective, Woburn, MA, USA) and measured with the mass spectrometer in a data-dependent mode using the Top12 method. The resolution of full MS and MS/MS scans was 70,000 and 17,500, respectively. The performance of the LC-MS/MS platform was monitored during data acquisition using the QCloud quality control system [53].

3.3. LC-MS/MS Data Processing

Data were processed using MaxQuant (version 1.6.7.0) [54] and searched by the integrated search engine, Andromeda [55], against Lacertidae SwissTrEMBL database (45,072 sequences, downloaded on 23 September 2020), supplemented with common protein contaminant sequences. iBAQ quantification was enabled and, except for this parameter, standard software settings were used, which included a false discovery rate (FDR) below 1% for peptide and protein identification. The resulting table with identified protein groups was further processed using Perseus [56]. Protein groups from the reverse database, common protein contaminants, as well as proteins only identified by site, were filtered out (with 1487 protein groups left). Data were log transformed and the matrix was filtered for protein groups that had valid iBAQ values in at least 10 out of 15 samples (667 protein groups). Means of the iBAQ values were calculated for each protein group to assess their abundance within femoral gland secretions. A heat map of the Pearson correlation coefficients, based on iBAQ values, was constructed with Perseus software to visually explore inter-individual variation in protein expression.
Protein groups were ranked according to their mean iBAQ values (from highest to lowest). Thirty proteins with the highest values (Top30) were manually grouped by function with the help of the information available in the UniProt database [34].
Functional enrichment analysis for the whole dataset (667 protein groups) was conducted using STRING [57,58]. Gene names were retrieved from the UniProt database for first major protein ID for each protein group. They were available for 586 proteins and all of them were uploaded to STRING. Human taxonomy was chosen to ensure a high degree of protein annotation. 510 of the genes were recognized by the software. Additionally, the DAVID bioinformatics resource [59,60] was applied, similarly using human taxonomy (506 recognized genes). The transfer of functional annotations between species is commonplace in proteomic and genomic studies [44,45,46,47]. Although the use of human reference databases (or other well-annotated databases) may increase the number of annotations, extrapolation of protein function to other organisms should be cautiously interpreted. Considering this limitation, we use a similar approach to provide tentative biological functions of FG proteins in sand lizards.

4. Conclusions

The use of mass spectrometry-based proteomics as a tool in chemical ecology may help identify and functionally annotate proteins found in the skin glands of vertebrates. We find that FG proteomes in sand lizards have a broad and complex composition, and reveal substantial variation in the total number of proteins among individuals. Nonetheless, some proteins were abundantly expressed across multiple individuals. At a functional level, FG proteins expressed at high intensity are related to the immunological response, lipid metabolism and antioxidant activity, among others, showing a striking similarity to protein groups identified in Galápagos marine iguanas. Although it is currently unknown whether the expression of immune-regulatory proteins occurs in FG secretions of other lizard species, our dataset demonstrates that these proteins are present in at least one other family of squamates (Lacertidae) beyond marine iguanas (family Iguanidae). Since marine iguanas are phylogenetically distant from lacertid lizards, it would be interesting to investigate whether other species secrete similar functional protein groups, and whether the expression of proteins related to immune response is a widespread phenomenon in squamate epidermal glands.
Proteins directly involved in intraspecific chemical communication were not identified in sand lizard FGs. However, some of the described proteins could have multiple functions, and thus further research is needed to validate the FG proteins found in sand lizards. Here, we provide a first assessment on the number, identity and tentative functions of FG proteins that should be explored in greater detail to understand other potential roles they may have in the context of intraspecific communication.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/molecules27072371/s1, File S1. Information on the number of identified proteins, protein groups that had valid iBAQ values for at least 10 out of 15 samples, proteins retrieved from the UniProt database together with the intensity rank, STRING analysis and values for Figure 3A; File S2. Results of the analysis with the DAVID tool.

Author Contributions

A.I. and M.P. designed the study; U.J., B.S.-R. and S.K.-K. performed the proteomics and bioinformatics analyses; B.Z. and M.P. took part in field work and logistics; A.I and M.P. prepared the manuscript with input from all the authors. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by a National Science Centre of Poland (Narodowe Centrum Nauki, NCN) OPUS grant awarded to AI (UMO-2017/25/B/NZ8/01498).

Institutional Review Board Statement

Lizards were captured under authorization issued by the Regional Directorate for Environmental Protection in Kraków (no: OP.6401.99.2020.PKw). All lizards were released at the same site of collection immediately after sampling.

Informed Consent Statement

Not applicable.

Data Availability Statement

The mass spectrometry data were deposited to the ProteomeXchange via the MassIVE repository with the dataset identifier PXD031566.

Acknowledgments

We want to thank the National Science Centre of Poland (Narodowe Centrum Nauki, NCN) for financial support (UMO-2017/25/B/NZ8/01498) and the Regional Directorate for Environmental Protection in Kraków for authorization to capture the lizards used in this study.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the author.

References

  1. Touhara, K. Sexual communication via peptide and protein pheromones. Curr. Opin. Pharmacol. 2008, 8, 759–764. [Google Scholar] [CrossRef] [PubMed]
  2. Wyatt, T.D. Proteins and peptides as pheromone signals and chemical signatures. Anim. Behav. 2014, 97, 273–280. [Google Scholar] [CrossRef]
  3. Roberts, S.A.; Simpson, D.M.; Armstrong, S.D.; Davidson, A.J.; Robertson, D.H.; McLean, L.; Beynon, R.J.; Hurst, J.L. Darcin: A male pheromone that stimulates female memory and sexual attraction to an individual male’s odour. BMC Biol. 2010, 8, 75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Houck, L.D. Pheromone communication in amphibians and reptiles. Annu. Rev. Physiol. 2009, 71, 161–176. [Google Scholar] [CrossRef] [Green Version]
  5. Maex, M.; Van Bocxlaer, I.; Mortier, A.; Proost, P.; Bossuyt, F. Courtship pheromone use in a model urodele, the Mexican axolotl (Ambystoma mexicanum). Sci. Rep. 2016, 6, 20184. [Google Scholar] [CrossRef] [Green Version]
  6. Van Bocxlaer, I.; Maex, M.; Treer, D.; Janssenswillen, S.; Janssens, R.; Vandebergh, W.; Proost, P.; Bossuyt, F. Beyond sodefrin: Evidence for a multi-component pheromone system in the model newt Cynops pyrrhogaster (Salamandridae). Sci. Rep. 2016, 6, 21880. [Google Scholar] [CrossRef]
  7. Baeckens, S.; Herrel, A.; Broeckhoven, C.; Vasilopoulou-Kampitsi, M.; Huyghe, K.; Goyens, J.; Van Damme, R. Evolutionary morphology of the lizard chemosensory system. Sci. Rep. 2017, 7, 10141. [Google Scholar] [CrossRef] [Green Version]
  8. Cooper, W.E. Correlated evolution of prey chemical discrimination with foraging, lingual morphology and vomeronasal chemoreceptor abundance in lizards. Behav. Ecol. Sociobiol. 1997, 41, 257–265. [Google Scholar] [CrossRef]
  9. Schwenk, K. Of tongues and noses: Chemoreception in lizards and snakes. Trends Ecol. Evol. 1995, 10, 7–12. [Google Scholar] [CrossRef]
  10. Cole, C.J. Femoral glands in lizards: A review. Herpetologica 1966, 22, 199–206. [Google Scholar]
  11. Martín, J.; Lopez, P. Pheromones and chemical communication in lizards. In Reproductive Biology and Phylogeny of Lizards and Tuatara; Rheubert, J.L., Siegel, D.S., Trauth, S.E., Eds.; CRC Press: Boca Raton, FL, USA, 2014; Volume 10, pp. 43–77. [Google Scholar]
  12. Mayerl, C.; Baeckens, S.; Van Damme, R. Evolution and role of the follicular epidermal gland system in non-ophidian squamates. Amphibia-Reptilia 2015, 36, 185–206. [Google Scholar] [CrossRef] [Green Version]
  13. Matthey, R. Caractères sexuels secondaires du lézard mâle (‘Lacerta agilis’). Bull. Soc. Vaud. Sci. Nat. 1929, 57, 71–81. [Google Scholar]
  14. Baeckens, S.; Whiting, M.J. Investment in chemical signalling glands facilitates the evolution of sociality in lizards. Proc. R. Soc. B 2021, 288, 20202438. [Google Scholar] [CrossRef] [PubMed]
  15. Alberts, A.C.; Pratt, N.C.; Phillips, J.A. Seasonal productivity of lizard femoral glands: Relationship to social dominance and androgen levels. Physiol. Behav. 1992, 51, 729–733. [Google Scholar] [CrossRef]
  16. Carazo, P.; Font, E.; Desfilis, E. Chemosensory assessment of rival competitive ability and scent-mark function in a lizard, Podarcis hispanica. Anim. Behav. 2007, 74, 895–902. [Google Scholar] [CrossRef]
  17. Khannoon, E.; Breithaupt, T.; El-Gendy, A.; Hardege, J.D. Sexual differences in behavioural response to femoral gland pheromones of Acanthodactylus boskianus. Herpetol. J. 2010, 20, 225–229. [Google Scholar]
  18. Martín, J.; López, P. Chemoreception, symmetry and mate choice in lizards. Proc. R. Soc. B 2000, 267, 1265–1269. [Google Scholar] [CrossRef]
  19. Martín, J.; López, P. Condition-dependent pheromone signaling by male rock lizards: More oily scents are more attractive. Chem. Senses 2010, 35, 253–262. [Google Scholar] [CrossRef] [Green Version]
  20. Alberts, A.C. Chemical properties of femoral gland secretions in the desert iguana, Dipsosaurus dorsalis. J. Chem. Ecol. 1990, 16, 13–25. [Google Scholar] [CrossRef]
  21. Baeckens, S. Evolution of animal chemical communication: Insights from non-model species and phylogenetic comparative methods. Belg. J. Zool. 2019, 149, 63–93. [Google Scholar] [CrossRef]
  22. Weldon, P.J.; Flachsbarth, B.; Schulz, S. Natural products from the integument of nonavian reptiles. Nat. Prod. Rep. 2008, 25, 738–756. [Google Scholar] [CrossRef]
  23. Alberts, A.C. Phylogenetic and adaptive variation in lizard femoral gland secretions. Copeia 1991, 1991, 69–79. [Google Scholar] [CrossRef]
  24. Mangiacotti, M.; Fumagalli, M.; Cagnone, M.; Viglio, S.; Bardoni, A.M.; Scali, S.; Sacchi, R. Morph-specific protein patterns in the femoral gland secretions of a colour polymorphic lizard. Sci. Rep. 2019, 9, 8412. [Google Scholar] [CrossRef]
  25. Mangiacotti, M.; Fumagalli, M.; Scali, S.; Zuffi, M.A.; Cagnone, M.; Salvini, R.; Sacchi, R. Inter-and intra-population variability of the protein content of femoral gland secretions from a lacertid lizard. Curr. Zool. 2017, 63, 657–665. [Google Scholar] [CrossRef] [Green Version]
  26. Mangiacotti, M.; Gaggiani, S.; Coladonato, A.J.; Scali, S.; Zuffi, M.A.L.; Sacchi, R. First experimental evidence that proteins from femoral glands convey identity-related information in a lizard. Acta Ethol. 2019, 22, 57–65. [Google Scholar] [CrossRef]
  27. Mangiacotti, M.; Baeckens, S.; Scali, S.; Martín, J.; Van Damme, R.; Sacchi, R. Evolutionary and biogeographical support for species-specific proteins in lizard chemical signals. Biol. J. Linn. Soc. 2021, 134, 912–928. [Google Scholar] [CrossRef]
  28. Tellkamp, F.; Lang, F.; Ibáñez, A.; Abraham, L.; Quezada, G.; Günther, S.; Looso, M.; Tann, F.J.; Müller, D.; Cemic, F.; et al. Proteomics of Galápagos marine iguanas links function of femoral gland proteins to the immune system. Mol. Cell. Proteom. 2020, 19, 1523–1532. [Google Scholar] [CrossRef]
  29. Blanke, I.; Fernley, H. The Sand Lizard—Between Light and Shadow; Laurenti-Verlag: Bielefeld, Germany, 2015. [Google Scholar]
  30. Olsson, M.; Madsen, T.; Nordby, J.; Wapstra, E.; Ujvari, B.; Wittsell, H. Major histocompatibility complex and mate choice in sand lizards. Proc. R. Soc. B 2003, 270, S254–S256. [Google Scholar] [CrossRef] [Green Version]
  31. Baeckens, S.; Martín, J.; García-Roa, R.; Pafilis, P.; Huyghe, K.; Van Damme, R. Environmental conditions shape the chemical signal design of lizards. Funct. Ecol. 2018, 32, 566–580. [Google Scholar] [CrossRef]
  32. Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized ppb-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 2008, 26, 1367–1372. [Google Scholar] [CrossRef]
  33. Schwanhäusser, B.; Busse, D.; Li, N.; Dittmar, G.; Schuchhardt, J.; Wolf, J.; Chen, W.; Selbach, M. Global quantification of mammalian gene expression control. Nature 2011, 473, 337–342. [Google Scholar] [CrossRef] [Green Version]
  34. The UniProt Consortium. UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res. 2021, 49, D480–D489. [Google Scholar] [CrossRef]
  35. Benes, P.; Vetvicka, V.; Fusek, M. Cathepsin D—Many functions of one aspartic protease. Crit. Rev. Oncol. Hematol. 2008, 68, 12–28. [Google Scholar] [CrossRef] [Green Version]
  36. Kolli, N.; Garman, S.C. Proteolytic activation of human cathepsin A. J. Biol. Chem. 2014, 289, 11592–11600. [Google Scholar] [CrossRef] [Green Version]
  37. Conus, S.; Simon, H.U. Cathepsins and their involvement in immune responses. Swiss Med. Wkly. 2010, 140, 4–11. [Google Scholar] [CrossRef]
  38. Reich, M.; Spindler, K.-D.; Burret, M.; Kalbacher, H.; Boehm, B.O.; Burster, T. Cathepsin A is expressed in primary human antigen-presenting cells. Immunol. Lett. 2010, 128, 143–147. [Google Scholar] [CrossRef]
  39. Hassan, M.I.; Shajee, B.; Waheed, A.; Ahmad, F.; Sly, W.S. Structure, function and applications of carbonic anhydrase isozymes. Bioorg. Med. Chem. 2013, 21, 1570–1582. [Google Scholar] [CrossRef]
  40. Gabirot, M.; López, P.; Martín, J.; De Fraipont, M.; Heulin, B.; Sinervo, B.; Clobert, J. Chemical composition of femoral secretions of oviparous and viviparous types of male common lizards Lacerta vivipara. Biochem. Syst. Ecol. 2008, 36, 539–544. [Google Scholar] [CrossRef]
  41. Mathivanan, S.; Ji, H.; Simpson, R.J. Exosomes: Extracellular organelles important in intercellular communication. J. Proteom. 2010, 73, 1907–1920. [Google Scholar] [CrossRef]
  42. Imparato, B.A.; Antoniazzi, M.M.; Rodrigues, M.T.; Jared, C. Morphology of the femoral glands in the lizard Ameiva ameiva (Teiidae) and their possible role in semiochemical dispersion. J. Morphol. 2007, 268, 636–648. [Google Scholar] [CrossRef]
  43. Khannoon, E.R.; Dollahon, N.R.; Bauer, A.M. Comparative study of the pheromone-manufacturing femoral glands in two sympatric species of lacertid lizards (Acanthodactylus). Zool. Sci. 2013, 30, 110–117. [Google Scholar] [CrossRef]
  44. Dietrich, M.A.; Irnazarow, I.; Adamek, M.; Jurecka, P.; Teich, L.; Rakus, K.; Kodzik, N.; Chadzińska, M.; Steinhagen, D.; Ciereszko, A. 2D-DIGE proteomic analysis of blood plasma reveals changes in immune-and stress-associated proteins following hormonal stimulation of carp males. Fish Shellfish Immunol. 2021, 118, 354–368. [Google Scholar] [CrossRef]
  45. Nynca, J.; Żarski, D.; Fröhlich, T.; Köster, M.; Bobe, J.; Ciereszko, A. Comparative proteomic analysis reveals the importance of the protective role of ovarian fluid over eggs during the reproduction of pikeperch. Aquaculture 2022, 548, 737656. [Google Scholar] [CrossRef]
  46. Rast, J.P.; D’Alessio, S.; Kraev, I.; Lange, S. Post-translational protein deimination signatures in sea lamprey (Petromyzon marinus) plasma and plasma-extracellular vesicles. Dev. Comp. Immunol. 2021, 125, 104225. [Google Scholar] [CrossRef]
  47. Tietgen, L.; Hagen, I.J.; Kleven, O.; Bernardi, C.D.; Kvalnes, T.; Norén, K.; Hasselgren, M.; Wallén, J.F.; Angerbjörn, A.; Landa, A.; et al. Fur colour in the Arctic fox: Genetic architecture and consequences for fitness. Proc. R. Soc. B 2021, 288, 20211452. [Google Scholar] [CrossRef]
  48. Penning, T.M.; Wangtrakuldee, P.; Auchus, R.J. Structural and functional biology of aldo-keto reductase steroid-transforming enzymes. Endocr. Rev. 2019, 40, 447–475. [Google Scholar] [CrossRef]
  49. Olsson, M.; Madsen, T.; Wapstra, E.; Silverin, B.; Ujvari, B.; Wittzell, H. MHC, health, color, and reproductive success in sand lizards. Behav. Ecol. Sociobiol. 2005, 58, 289–294. [Google Scholar] [CrossRef]
  50. Ruff, J.S.; Nelson, A.C.; Kubinak, J.L.; Potts, W.K. MHC signaling during social communication. In Self and Nonself; Springer: New York, NY, USA, 2012; pp. 290–313. [Google Scholar]
  51. Baeckens, S.; Martín, J.; García-Roa, R.; Van Damme, R. Sexual selection and the chemical signal design of lacertid lizards. Zool. J. Linn. Soc. 2018, 183, 445–457. [Google Scholar] [CrossRef] [Green Version]
  52. Kopena, R.; López, P.; Martín, J. Relative contribution of dietary carotenoids and vitamin E to visual and chemical sexual signals of male Iberian green lizards: An experimental test. Behav. Ecol. Sociobiol. 2014, 68, 571–581. [Google Scholar] [CrossRef]
  53. Chiva, C.; Olivella, R.; Borras, E.; Espadas, G.; Pastor, O.; Sole, A.; Sabido, E. QCloud: A cloud-based quality control system for mass spectrometry-based proteomics laboratories. PLoS ONE 2018, 13, e0189209. [Google Scholar] [CrossRef] [Green Version]
  54. Tyanova, S.; Temu, T.; Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 2016, 11, 2301–2319. [Google Scholar] [CrossRef]
  55. Cox, J.; Neuhauser, N.; Michalski, A.; Scheltema, R.A.; Olsen, J.V.; Mann, M. Andromeda: A peptide search engine integrated into the MaxQuant environment. J. Proteome Res. 2011, 10, 1794–1805. [Google Scholar] [CrossRef]
  56. Tyanova, S.; Temu, T.; Sinitcyn, P.; Carlson, A.; Hein, M.Y.; Geiger, T.; Mann, M.; Cox, J. The Perseus computational platform for comprehensive analysis of (prote) omics data. Nat. Methods 2016, 13, 731–740. [Google Scholar]
  57. STRING. Functional Protein Association Networks. Available online: https://string-db.org/ (accessed on 11 March 2021).
  58. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P. STRING v11: Protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [Green Version]
  59. DAVID. Bioinformatics Resources (2021 Update). Available online: https://david.ncifcrf.gov/home.jsp (accessed on 18 March 2022).
  60. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar]
Molecules 27 02371 g001 550
Figure 1. (A) Macroscopic aspect of femoral glands (highlighted with red circles) in a male sand lizard (Lacerta agilis). Photo by A. Ibáñez. (B) Total number of proteins identified in each femoral gland sample. Numbers in the abscissa are sample ID numbers (i.e., sample identifiers).
Figure 1. (A) Macroscopic aspect of femoral glands (highlighted with red circles) in a male sand lizard (Lacerta agilis). Photo by A. Ibáñez. (B) Total number of proteins identified in each femoral gland sample. Numbers in the abscissa are sample ID numbers (i.e., sample identifiers).
Molecules 27 02371 g001
Molecules 27 02371 g002 550
Figure 2. Top left: Density plot of the molecular weight (kDa) for protein groups with a valid iBAQ in at least 10 out of 15 samples (i.e., a total of 667 proteins selected). Heat map showing Pearson correlation coefficients for the 15 samples based on the iBAQ values of the 667 selected proteins. Scale shows color correspondence for the Pearson coefficient values.
Figure 2. Top left: Density plot of the molecular weight (kDa) for protein groups with a valid iBAQ in at least 10 out of 15 samples (i.e., a total of 667 proteins selected). Heat map showing Pearson correlation coefficients for the 15 samples based on the iBAQ values of the 667 selected proteins. Scale shows color correspondence for the Pearson coefficient values.
Molecules 27 02371 g002
Molecules 27 02371 g003 550
Figure 3. (A) Protein intensity rank (x-axis) based on the mean log2 iBAQ values (y-axis). Proteins of selected functional groups are color-coded (top: immune related and epidermal identity; middle: lipid related; bottom: ribosome and proteasome). (B) Detailed functions for Top30 proteins. The pie chart shows the number of proteins per category (see Table 1). (C) MS/MS spectrum of the peptide SYALTPQQYALK used for the identification of Cathepsin D (D3JHM6_9SAUR). The score value for this spectrum was 207.29.
Figure 3. (A) Protein intensity rank (x-axis) based on the mean log2 iBAQ values (y-axis). Proteins of selected functional groups are color-coded (top: immune related and epidermal identity; middle: lipid related; bottom: ribosome and proteasome). (B) Detailed functions for Top30 proteins. The pie chart shows the number of proteins per category (see Table 1). (C) MS/MS spectrum of the peptide SYALTPQQYALK used for the identification of Cathepsin D (D3JHM6_9SAUR). The score value for this spectrum was 207.29.
Molecules 27 02371 g003
Table
Table 1. Thirty most abundant proteins (Top30) detected in the femoral gland secretions of sand lizards. Top30 proteins were manually grouped by function with the help of the information available in the UniProt database [34]. IR: Intensity Rank.
Table 1. Thirty most abundant proteins (Top30) detected in the femoral gland secretions of sand lizards. Top30 proteins were manually grouped by function with the help of the information available in the UniProt database [34]. IR: Intensity Rank.
Protein IDProtein NameGene NameLengthOrganismIR *IR in Tellkamp et al. **Comment
Cathepsins
D3JHM6_9SAURCathepsin D (Fragment)CTSD24Lacerta schreiberi110
A0A670IA46_PODMUPeptidase A1 domain-containing proteinCTSD399Podarcis muralis310
A0A670IH97_PODMUCarboxypeptidase (EC 3.4.16.-)CTSA530Podarcis muralis14112
Chaperones
A0A670J0X0_PODMUUncharacterized proteinPDIA6449Podarcis muralis4not FG-specific ***
A0A670KPF3_PODMU78 kDa glucose-regulated protein (Binding-immunoglobulin protein) (Heat shock protein 70 family protein 5) (Heat shock protein family A member 5) (Immunoglobulin heavy chain-binding protein)HSPA5654Podarcis muralis517
A0A670I8P5_PODMUSHSP domain-containing protein 169Podarcis muralis8not FG-specific ***belongs to the small heat shock protein (HSP20) family
A0A670K8T4_PODMUUncharacterized proteinHSPA8646Podarcis muralis11not FG-specific ***
Lipid related
A0A670JAW3_PODMUGP-PDE domain-containing proteinGDPD3318Podarcis muralis6not FG-specific ***
A0A670JD81_PODMUBeta-hexosaminidase (EC 3.2.1.52)HEXA522Podarcis muralis10105
A0A670HXK7_PODMUPhospholipase A2 (EC 3.1.1.4) 731Podarcis muralis24158
A0A670KKI2_PODMUAldo_ket_red domain-containing protein 278Podarcis muralis28not FG-specific ***members of aldo-keto reductase family are involved in steroid metabolic processes
A0A670IG25_PODMUInositol-1-monophosphatase (EC 3.1.3.25)IMPA1347Podarcis muralis29not FG-specific ***
Cytoskeletal, Keratins
A0A1D9CFN7_9SAURBeta-actin 375Eremias argus7261
A0A670KC37_PODMUIF rod domain-containing proteinLOC114593000611Podarcis muralis94, 7, 11, 12, 25, 46, 109, 764keratin
A0A670JPK3_PODMUIF rod domain-containing proteinLOC114581896472Podarcis muralis274, 7, 11, 12, 25, 46, 109, 764keratin
A0A670JN74_PODMUIF rod domain-containing proteinLOC114582343462Podarcis muralis304, 7, 11, 12, 25, 46, 109, 764keratin
Histones
A0A670K9E9_PODMUHistone H2ALOC114582712129Podarcis muralis15not FG-specific ***
A0A670JD32_PODMUHistone H3LOC114581771136Podarcis muralis16not FG-specific ***
A0A670K7U8_PODMUHistone H4 101Podarcis muralis21not FG-specific ***
A0A670K3F4_PODMUHistone H2BLOC114581794126Podarcis muralis22not FG-specific ***
Antioxidant defense
A0A670I2N3_PODMUSuperoxide dismutase [Cu-Zn] (EC 1.15.1.1)SOD1155Podarcis muralis17not FG-specific ***
A0A670I4X7_PODMUThioredoxin domain-containing proteinPRDX4276Podarcis muralis25113
Others
A0A670JGA4_PODMUCarbonic anhydrase (EC 4.2.1.1)LOC114583714316Podarcis muralis22pH regulation
A0A670JY74_PODMUGlyceraldehyde-3-phosphate dehydrogenase (EC 1.2.1.12)GAPDH333Podarcis muralis12not FG-specific ***
A0A670HW54_PODMU40S ribosomal protein S27a (Ubiquitin carboxyl extension protein 80) (Ubiquitin-40S ribosomal protein S27a)RPS27A156Podarcis muralis13204
A0A670ICE5_PODMUElongation factor 1-alphaEEF1A1462Podarcis muralis18398
A0A670HUL4_PODMUVoltage-dependent anion-selective channel protein 2VDAC2283Podarcis muralis19not FG-specific ***
A0A670HU92_PODMUMannosyl-glycoprotein endo-beta-N-acetylglucosaminidase (EC 3.2.1.96)ENGASE811Podarcis muralis2094glycoprotein modification
A0A670K0C8_PODMUS_100 domain-containing proteinS100A14102Podarcis muralis2347cell survival and apoptosis/epidermal identity
A0A670JM37_PODMUSERPIN domain-containing proteinSERPINE2386Podarcis muralis26416protease inhibitor
* iBAQ Intensity Rank—based on the mean value from at least 10 out of 15 samples. ** within femoral gland (FG) specific proteins in Galápagos marine iguanas [28]. *** identified in FG sample but not classified as FG-specific in Galápagos marine iguanas [28].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ibáñez, A.; Skupien-Rabian, B.; Jankowska, U.; Kędracka-Krok, S.; Zając, B.; Pabijan, M. Functional Protein Composition in Femoral Glands of Sand Lizards (Lacerta agilis). Molecules 2022, 27, 2371. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27072371

AMA Style

Ibáñez A, Skupien-Rabian B, Jankowska U, Kędracka-Krok S, Zając B, Pabijan M. Functional Protein Composition in Femoral Glands of Sand Lizards (Lacerta agilis). Molecules. 2022; 27(7):2371. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27072371

Chicago/Turabian Style

Ibáñez, Alejandro, Bozena Skupien-Rabian, Urszula Jankowska, Sylwia Kędracka-Krok, Bartłomiej Zając, and Maciej Pabijan. 2022. "Functional Protein Composition in Femoral Glands of Sand Lizards (Lacerta agilis)" Molecules 27, no. 7: 2371. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27072371

Article Metrics

Back to TopTop