# Within and among population differences in cuticular hydrocarbons in the seabird tick Ixodes uriae

Marlène Dupraz, Chloé Leroy, Thorkell Lindberg Thórarinsson, Patrizia

D'ettorre, Karen Mccoy

► To cite this version:

Marlène Dupraz, Chloé Leroy, Thorkell Lindberg Thórarinsson, Patrizia D'ettorre, Karen Mccoy. Within and among population differences in cuticular hydrocarbons in the seabird tick Ixodes uriae. Peer Community Journal, 2022, 2, pp.e51. 10.24072/pcjournal.164. hal-03864261

HAL Id: hal-03864261

https://hal.science/hal-03864261

Submitted on 20 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

Published
2022-08-30

**Cite as**
Marlène Dupraz, Chloé Leroy, Thorkell Lindberg Thórarinsson, Patrizia d'Ettorre and Karen D. McCoy (2022) *Within and among population differences in cuticular hydrocarbons in the seabird tick Ixodes uriae*, Peer Community Journal, 2: e51.

**Correspondence**
marlenedupraz@gmail.com

**Peer-review**
Peer reviewed and
recommended by
PCI Zoology,
https://doi.org/10.24072/pci.zool.100014

(cc) BY

This article is licensed
under the Creative Commons
Attribution 4.0 License.

# Within and among population differences in cuticular hydrocarbons in the seabird tick *Ixodes uriae*

Marlène Dupraz<sup>1,1</sup>, Chloé Leroy<sup>1,2</sup>, Thorkell Lindberg Thórarinsson<sup>3</sup>, Patrizia d'Ettorre<sup>1,2</sup>, and Karen D. McCoy<sup>1,1</sup>

Volume 2 (2022), article e51

https://doi.org/10.24072/pcjournal.164

## Abstract

The hydrophobic layer of the arthropod cuticle acts to maintain water balance, but can also serve to transmit chemical signals via cuticular hydrocarbons (CHC), essential mediators of arthropod behavior. CHC signatures typically vary qualitatively among species, but also quantitatively among populations within a species, and have been used as taxonomic tools to differentiate species or populations in a variety of taxa. Most work in this area to date has focused on insects, with little known for other arthropod groups such as ticks. The worldwide distribution and extensive host-range of the seabird tick *Ixodes uriae* make it a good model to study the factors influencing CHC composition. Genetically differentiated host-races of *I. uriae* have evolved across the distribution of this species but the factors promoting sympatric population divergence are still unknown. To test for a potential role of host-associated CHC in population isolation, we collected *I. uriae* specimens from two of its seabird hosts, the Atlantic puffin (*Fratercula arctica*) and the common guillemot (*Uria aalge*) in different colonies in Iceland. Using gas-chromatography and mass-spectrometry, we detected a complex cuticular mixture of 22 hydrocarbons, including n-alkanes, methyl-alkanes and alkenes ranging from 17 to 33 carbons in length. We found that each population had a distinct CHC profile. The host group explained the greatest amount of population divergence, with long-chain hydrocarbons being more abundant in puffin tick populations compared to guillemot tick populations. Future work will now be required to test whether the different CHC signals reinforce assortative mating, thereby playing a role in generating *I. uriae* population divergence patterns, and to evaluate diverse hypotheses on the origin of distinct population signatures.

<footnote n="1"><sup>1</sup>MIVEGEC, IRD, CNRS, Univ. Montpellier, Domaine La Valette - 900, rue Jean François BRETON 34090 Montpellier, France, <sup>2</sup>Laboratoire d'Ethologie Expérimentale et Comparée UR 4443 (LEEC), Université Sorbonne Paris Nord, F-93430 Villetaneuse, France, <sup>3</sup>Northeast Iceland Nature Research Centre, 640 Húsavík, Iceland</footnote>
Introduction

The arthropod cuticle acts as both an exoskeleton and a barrier from the external environment (Andersen,
1979). The outermost layer, the epicuticle, is covered by a lipid layer (Lockey, 1988) made up of esters,
carboxylic acids, alcohols, carbonyls and long-chain hydrocarbons (Andersen, 1979) that protects the organism
from desiccation (Filshie, 1982). However, cuticular hydrocarbons (CHC) are also involved in chemical
communication, serving as sex pheromones, kairomones and/or signature mixtures allowing recognition of
social identity (van Zweden and d'Ettorre, 2010; Wyatt, 2010). Using analytic techniques such as gas
chromatography-mass spectrometry (GC-MS) and MALDI-TOF mass spectrometry, researchers have described
hydrocarbons of up to 70 carbons in chain length in insects, principally *n*-alkanes, methyl-branched alkanes
and alkenes (Blomquist and Bagnères, 2010). The array of CHC on the cuticle constitute a species-specific
chemical signature, varying qualitatively between species and quantitatively within species (Lockey, 1988).
CHC patterns are genetically controlled, but the relative abundance of particular components can be linked to
environmental conditions (Estrada-Peña et al., 1993; Gibbs et al., 1991). CHC patterns have been used as
taxonomic tools to characterize hundreds of arthropod species (Howard and Blomquist, 2005), and to
discriminate closely-related populations (Bagnères et al., 1991; Bartelt et al., 1986; Jallon and David, 1987;
Kruger et al., 1991; Simmons et al., 2014).

Ticks are hematophagous arthropods widely distributed across the globe and parasitize a diverse array of vertebrate species (McCoy and Boulanger, 2015). However, little work has been performed to describe CHC profiles and their variation among species and populations. The only studies to date have used CHC profiles in an attempt to differentiate closely-related taxa (Estrada-Peña et al., 1992, 1994, 1996; Hunt, 1986; Estrada-Peña and Dusbabek, 1993). However, tick CHCs may play essential roles in several aspects of tick life histories. First, the tick epicuticle is perforated with numerous channels providing a large surface for exchange with the external environment. As most tick species spend the major part of their life cycle in the off-host environment, maintaining water balance across this surface under different environmental conditions will directly dictate survival (Randolph and Storey, 1999); the presence of CHCs likely plays an important role in this. Second, as obligate parasites, access to the vertebrate host for the bloodmeal is a key aspect of the tick life cycle and, as such, ticks have adapted key traits to locate and successfully exploit their host (McCoy and Boulanger, 2015), as for example Haller's organ, a unique chemosensory structure located on the 1st pairs of legs (Leonovich, 2004). Shimshoni et al. (2013) found *Rhipicephalus* tick species had different cuticular fatty acid compositions in relation to host use. However, whether these differences are linked to adaptive survival or are by-product variation due to the host resource is unknown as yet. Finally, ticks aggregate both on hosts and in the off-host environment (Randolph, 1998). This behavior is thought to facilitate blood feeding on the host and increase survival in the off-host environment. It may also enable ticks to find appropriate mates for reproduction. Some ticks, such as the tropical bont tick *Amblyomma variegatum*, produce a multicomponent pheromone to provoke aggregation (Schöni et al., 1984) but the potential role of these pheromones in assortative mating has never been examined.

Seabird ticks are generally nidicolous, exploiting different local host species that use diverse micro-habitats within the colony (Dietrich et al., 2011). Due to this diversity, these ticks may experience diverse selective pressures coming from both the hosts and from the temperature and humidity conditions of the nest micro-habitat. In particular, *Ixodes uriae*, a tick associated with seabird colonies in the polar regions of both hemispheres, is known to form host-specific races that show genetic (McCoy et al., 2005, 2003, 2001) and morphological (Dietrich et al., 2013) differences in relation to host use. Differential performance on alternative hosts has also been experimentally demonstrated (Dietrich et al., 2014). Nevertheless, the factors driving divergence within this species have yet to be specifically identified. A potential role for isolating mechanisms that lead to assortative mating based on host use have been suggested (McCoy et al., 2013).

Here, we examine the degree of cuticular hydrocarbon diversity in *I. uriae* and the factors influencing these
chemical signatures. Based on current knowledge, we predicted that cuticular hydrocarbon patterns could vary
in ticks 1) exploiting different host species and 2) sampled in different geographic locations. If host exploitation
modifies CHC diversity and abundance, we expected that signatures in ticks from the same seabird host in
different locations should be more similar than signatures in ticks from different host species in the same
geographic location. To test these predictions, we collected *I. uriae* specimens in the nesting area of two
seabird host species, the common guillemot (CG) *Uria aalge*, and the Atlantic puffin (PF) *Fratercula arctica*, at
three locations in Iceland. We then extracted cuticular hydrocarbons and analyzed them using gas chromatography-mass spectrometry (GC-MS).

## Methods

### Population samples

Flat adult female ticks were collected off-host in three sites in Iceland in June 2016: common guillemot ticks (CG) were collected under rocks in the middle of the colony in Langanes (L) (66°22'07.2"N 14°38'33.0"W) and on Grimsey Island (GRI) (66°32'57.3"N 17°59'31.1"W); puffin ticks (PF) were collected in burrows on Lundy Island (LI) (66°06'53.2"N 17°22'13.1"W) and on Grimsey Island (GRI) (66°32'39.1"N 18°01'12.9"W) (Fig.1). Ticks were kept alive in plastic tubes before hydrocarbon extraction. Four replicates of 10 flat ticks were extracted for each host species at each site.

### Hydrocarbon extraction

Several extraction attempts were made with ticks kept alive or stored in alcohol and with only one individual or in batches of five, but the quantities of extracted cuticular compounds were too low to use in analyses. Cuticular compounds were therefore extracted by immersing 10 living flat female ticks in 200 µl of pentane (HPLC grade, Sigma-Aldrich) in glass vials. The vials were agitated for 1 minute, set down to rest for 5 minutes and then re-agitated for 1 minute. Ticks were then removed from the vials and preserved in 70% ethanol. The pentane was then evaporated off over 30 minutes to obtain dried eluent and the vials were closed and kept at 4°C until analysis. A negative control containing only 200 µl of pentane was added to each extraction session to control for potential contamination. We did not use an internal standard because the aim of our analysis was not to quantify the absolute amount of CHCs, but to compare the relative proportions of CHCs among samples.

### Chemical analyses

Cuticular hydrocarbons are not very volatile, but to reduce variation during the evaporation process, all samples were treated in the same way, placed under the same fume hood in a temperature controlled room. Samples were re-diluted in 40µl of pentane and 3µl were injected into an Agilent Technologies 7890A gas chromatograph (capillary column: Agilent HP-5MS, 30 m × 0.25 mm × 0.25 µm; split-splitless injector; carrying helium gas at 1 mL/min) coupled with an Agilent 5975C mass spectrometer with 70 eV electron impact ionization. The oven temperature was programmed at 70°C for 1min, and was increased at 30°C/min to 200°C, then to 320°C at 5°C/min and held for 5 min. Compounds were then identified (raw data available on Zenodo at https://doi.org/10.5281/zenodo.6497483) and sorted on the basis of their mass spectra and retention time by comparison with standards (alkane standard solutions, Sigma-Aldrich, Saint Louis, MO, USA) and published spectra (NIST Library). The abundance of each detected cuticular hydrocarbon was used to build a box plot using the ggplot 2 package with R software (version 4.1.1) (script is available on Zenodo at https://doi.org/10.5281/zenodo.7018260). The areas under the peaks were extracted using the Agilent MSD ChemStation software (E.02.01.1177) and were used to calculate the relative amount of each hydrocarbon to normalize data coming from different origins and of different qualities. The mean peak area of the four population replicates was calculated for each hydrocarbon to obtain a colony average.

### Multivariate analyses

For analytical purposes, if the relative abundance of a CHC was 0, the value was replaced by 0.00001 which is several times lower than the lowest quantity found for a CHC in our dataset. The data were then transformed using a centered log ratio procedure (data available in supplementary materials at https://doi.org/10.5281/zenodo.5889077) and analyzed by Partial Least Squares coupled with a Discriminant analysis (PLS-DA) using the R software (v 3.4.3) and the "RVAidememoire" package (Hervé, 2014) (https://cran.r-project.org/web/packages/RVAideMemoire/index.html) (script is available on Zenodo at https://doi.org/10.5281/zenodo.7018260). PLS-DA is an appropriate analysis when the dataset contains fewer groups than explanatory variables, as in the case of the present quantitative dataset. PLS-DA is a supervised technique, meaning that class memberships of the CHCs need to be predefined. Here, we only used the first eight axes produced by the PLS to perform two PLS-DA tests: the first analysis considered the four population samples as four different classes (LCG, GCG, LiPF, GPF, outlined on Fig. 1).
<figure><img src="image_2.png" /><figcaption>**Figure 1:** Map showing the sampling locations for *Ixodes uriae* ticks and hosts in Iceland in June 2016. PF refers to the Atlantic puffin (*Fratercula arctica*) and CG to the common guillemot (*Uria aalge*). The acronyms of the sampling locations are as follows: Grimsey Island - G; Lundey Island - Li; Llanganes - L.</figcaption></figure>
The second analysis was based on two classes only, representing the two host types (PF, CG). The number of significant PLS components was determined by cross model validation. We also calculated a numerical value representing the importance of each CHC variable in the projection (VIP), i.e. VIP values larger than 1 are most influential (Hervé, 2014).

## Results

### Hydrocarbon profiles

We detected a complex pattern of 22 hydrocarbons and two non-identified components (I1 and I2) in the cuticle of *I. uriae* ticks. The hydrocarbon mixture was composed of: 14 *n*-alkanes ranging from 17 to 33 carbon atoms in length, 4 monomethyl-alkanes with 17 to 24 carbon atoms and 4 alkenes with 21 to 27 carbon atoms (Table 1).

All hydrocarbons were shared by the four tick populations, except 2meC<sub>24</sub> and C<sub>33</sub> which were only found in puffin ticks (Fig.2).

A high degree of variation among samples in mean abundance of cuticular components was obvious (Fig.2). Hydrocarbons C<sub>22:1</sub>, *n*-C<sub>29</sub> and *n*-C<sub>31</sub> were the most predominant, whereas many hydrocarbons were detected in low quantities, as for example *n*-C<sub>20</sub>, 2meC<sub>22</sub> and *n*-C<sub>26</sub>.

<table><caption>**Table 1:** List of detected compounds in the cuticle of *I. uriae* adult female ticks, including name and type.</caption><thead><tr><td>Cuticular component</td><td>Name</td><td>Type</td></tr></thead><tbody><tr><td>C<sub>17</sub></td><td>n-heptadecane</td><td>alkane</td></tr><tr><td>meC<sub>17</sub></td><td>methylheptadecane</td><td>methyl-branched alkane</td></tr><tr><td>C<sub>18</sub></td><td>n-octadecane</td><td>alkane</td></tr><tr><td>C<sub>20</sub></td><td>n-eicosane</td><td>alkane</td></tr><tr><td>C<sub>21:1</sub></td><td>heneicosane</td><td>alkene</td></tr><tr><td>C<sub>21</sub></td><td>n-heneicosane</td><td>alkane</td></tr><tr><td>C<sub>22:1</sub></td><td>docosene</td><td>alkene</td></tr><tr><td>C<sub>22</sub></td><td>n-docosane</td><td>alkane</td></tr><tr><td>2meC<sub>22</sub></td><td>2methyldocosane</td><td>methyl-branched alkane</td></tr><tr><td>C<sub>23</sub></td><td>n-tricosane</td><td>alkane</td></tr><tr><td>C<sub>24</sub></td><td>n-tetracosane</td><td>alkane</td></tr><tr><td>9meC<sub>24</sub></td><td>9methyltetracosane</td><td>methyl-branched alkane</td></tr><tr><td>2meC<sub>24</sub></td><td>2methyltetracosane</td><td>methyl-branched alkane</td></tr><tr><td>C<sub>25:1</sub></td><td>pentacosene</td><td>alkene</td></tr><tr><td>C<sub>25</sub></td><td>n-pentacosane</td><td>alkane</td></tr><tr><td>C<sub>26</sub></td><td>n-hexacosane</td><td>alkane</td></tr><tr><td>C<sub>27:1</sub></td><td>heptacosene</td><td>alkene</td></tr><tr><td>C<sub>27</sub></td><td>n-heptacosane</td><td>alkane</td></tr><tr><td>C<sub>28</sub></td><td>n-octacosane</td><td>alkane</td></tr><tr><td>C<sub>29</sub></td><td>n-nonacosane</td><td>alkane</td></tr><tr><td>C<sub>31</sub></td><td>n-hentriacontane</td><td>alkane</td></tr><tr><td>C<sub>33</sub></td><td>n-tritriacontane</td><td>alkane</td></tr><tr><td>I1</td><td>Non-identified</td><td>-</td></tr><tr><td>I2</td><td>Non-identified</td><td>-</td></tr></tbody></table>
The abundance pattern of cuticular hydrocarbons varied between ticks of the two host species (Fig.3), but also among the same host in different colony sites. Nevertheless, results of all pairwise profile comparisons were non-significant. No CHC was specific to CG samples, although *n*-C₁₇ and C₂₂:₁ were highly abundant in the GRI colony (Fig.2 and 3). In contrast, long chain cuticular hydrocarbons tended to be present in both PF samples: *n*-C₂₇, *n*-C₂₉ and *n*-C₃₁. CHC components tended to be in lower overall abundance in samples from Langanes (Langanes CG), whereas samples from Lundey Island (Lundey PF) had the highest overall abundance (Fig. 2).

<figure><img src="image_3.png" /><figcaption>**Figure 2:** Boxplot representing the mean abundance of each detected cuticular hydrocarbon from *I. uriae* ticks. Compounds ranged from 17 to 33 carbon atoms in length. Mean abundance was calculated based on values from four replicate tick pools. First and third quartiles, median and data points are shown. See Table 1 for CHC abbreviations. PF refers to the Atlantic puffin (*Fratercula arctica*) and CG to the Common guillemot (*Uria aalge*).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 3:** Gas-chromatograms showing the comparison of cuticular hydrocarbon profiles for pools of 10 female ticks from Atlantic puffins (PF: black line) and Common guillemots (CG: orange line) on Grimsey Island, Iceland.</figcaption></figure>
## PLS-DA analyses

### Population samples

Cross model validation for population samples showed that 50.6% of samples were assigned to the population of origin. PLS-DA analysis showed a discriminating population effect (p=0.008) separating the two host species on first axis (Fig.4). The three first axes explained respectively 11.94, 7.99 and 3.99 % of the total variance among samples.

<figure><img src="image_5.png" /><figcaption>**Figure 4:** Graphical representation of the four population samples on the two first axes of the PLS-DA analysis. The first and second axes explained respectively 11.94 and 7.99% of the total variance among samples (p=0.008). The acronyms of the sampling locations are as follows: Grimsey - G; Langanes - L; Lundey Island - Li. PF refers to Atlantic puffin (*Fratercula arctica*) and CG to Common guillemot (*Uria aalge*).</figcaption></figure>

Hydrocarbons *n*-C₃₃, 2meC₂₄ and *n*-C₂₉ appeared as the most influential components separating the four population samples (Fig.5).

### Host types

The cross-validation test for host type revealed that 75% of the tick pools were assigned to the host of origin (PF, CG). However, the PLS-DA analysis revealed no significant difference in CHC profiles between the host types, although a tendency was evident (PLS-DA: p=0.064).
<figure><img src="image_6.png" /><figcaption>**Figure 5:** VIP classification showing the influence of CHC in discriminating the four population samples. VIP values larger than 1 (dashed red line) are most influential. See Table 1 for CHC abbreviations.</figcaption></figure>

# Discussion

## Cuticle composition and function

Cuticular hydrocarbons, including linear and methyl-branched alkanes and alkenes, are biosynthesized from an acetyl-CoA molecule in specialized secretory cells, the oenocytes which are mainly located in the epidermis of insects (Howard and Blomquist, 2005). They are shuttled through the hemolymph to the epicuticular surface via specialized pore canals penetrating the cuticular layers (Holze et al, 2021). Once fixed on the arthropod cuticle within a complex mixture of alcohols, esters, aldehydes, fatty acids, etc, hydrocarbons help prevent desiccation and serve in chemical communication, constituting essential mediators of insect behavior (Blomquist and Bagnères, 2010). Using GC-MS techniques to analyze extracts from female *I. uriae* ticks, we detected a complex mixture of cuticular hydrocarbons containing linear and monomethyl-alkanes and alkenes ranging from 17 to 33 carbon atoms. The qualitative composition is similar to that reported for many other arthropods, comprising predominantly linear alkanes (n-C23, n-C25, n-C27, n-C29 and n-C31) (Howard and Blomquist, 2005; Lockey, 1988). A majority of the detected cuticular hydrocarbons were already reported in other hard ticks: *I. persulcatus* (Tkachev et al., 2000), *Amblyomma variegatum* (Estrada-Peña et al., 1994) and *Rhipicephalus* spp. (Estrada-Peña et al., 1992). Nevertheless, alkenes were only detected once in low quantities in *Rhipicephalus* spp. (Estrada-Peña et al., 1992). Here,
we report the presence of alkenes as heneicosene, pentacosene, heptacosene and a high quantity of docosene, particularly in one of the CG samples. Alkenes were demonstrated to act as sex pheromones in the Alfalfa leaf-cutter bee *Megachile rotundata* and the rove beetle *Aleochara curtula* (Paulmier et al., 1999; Peschke and Metzler, 1987). Heneicosene is also described as an aggregation pheromone in *Drosophila* species (Bartelt et al., 1988; Bartelt and Jackson, 1984). Pentacosene is implicated in the mating process of different fly species as a stimulant pheromone (Uebel et al., 1978). The presence of these alkenes in *I. uriae* could correspond to the biological state of the female ticks when they were collected. However, as mating in this tick species frequently takes place prior to feeding on the host, and in particular when nymphal ticks emerge as flat females (McCoy and Tirard, 2002), we cannot conclude on the mating status of collected females. Indeed, some of the profile variation in our data may be due to the inclusion of females in different reproductive states. Experimental studies will now be necessary to determine the potential role of the observed CHC compounds in mating behavior.

We found also large amount of long-chain hydrocarbons in PF samples (n-C29, n-C31 and n-C33). Mixed with other compounds, n-C31 was demonstrated to induce copulation in males of the stable fly *Stomoxys calcitrans* (Uebel et al., 1975). Other long chain hydrocarbons, including n-C23 to n-C31, produced in large quantities and acting in combination, were found to serve in colony recognition by bumblebees *Bombus terrestris* (Rottler et al., 2013). Studies have also highlighted the importance of the environment in the acquisition of new hydrocarbons (d'Ettorre et al., 2006; d'Ettorre et al., 2002), acting for example as a recognition signal of their own nest for the social wasp *Polistes metricus* Say (1981) (Singer and Espelie, 1996; Espelie et al., 1990). Moreover, 2meC24, which we only detected in PF CHC profiles, has been shown to serve as a contact pheromone in peach twig borers *Anarsia lineatella* (Schlamp, 2005). As puffin burrows are often deep, and densely distributed (Harris and Wanless, 2011), the production of these cuticular hydrocarbons in *I. uriae* may enable ticks to find conspecifics in the host nesting environment. The use of this particular blend of CHCs may also help ticks to find individuals that smell similar, favoring assortative mating (van Zweden and d'Ettorre, 2010). The production of a specific CHC blend may not be necessary for CG ticks because guillemots breed in extremely dense numbers on cliff ledges with no constructed nest; ticks tend to aggregate under and around rocks on the cliff ledge such that finding a mate from the same host type may be easier compared to ticks exploiting puffin hosts. More generally, the presence of complex compounds in the CHC pattern of *I. uriae* suggests that chemical communication may be important in this tick species, potentially enabling ticks to find a suitable mating partner that would enable them to maintain host-adapted gene complexes (Sonenshine and Roe, 2014).

## Alternative hypotheses to explain the origin of CHC patterns

As expected, chemical analyses revealed that each tick population had a distinct CHC profile, but that specific CHCs were also associated with different host types. In particular, n-C33, 2meC24 and n-C29 were most frequently or exclusively detected in PF populations. The quantitative variability among detected hydrocarbons could be related to different factors.

First, aging and development have been demonstrated to impact cuticular hydrocarbon patterns in different taxa (Desena et al., 1999; Ichinose and Lenoir, 2009). For example, aging favors the production of longer hydrocarbon chains and decreases sexual attractiveness in *Drosophila melanogaster* (Kuo et al., 2012). The high quantities of long-chain hydrocarbons (n-C27, n-C29, n-C31) observed in the puffin tick samples could be a consequence of age variability between host-associated tick samples. Samples used in the present study were collected in the field and the relative age of the specimens could not be determined. However, as tick activity is relatively synchronous with seabird breeding, we did not expect the overall timing in adult female activity to differ in a systematic way between ticks exploiting the different host species, nor among distinct colony locations within Iceland.

Second, climatic conditions can also shape CHC composition. Differences in cuticular hydrocarbon composition were observed in geographically distant tick populations of *Amblyomma variegatum* (Estrada-Peña et al., 1994) and some extreme climatic parameters were shown to be correlated with methylated-alkanes in *Rhipicephalus sanguineus* (Estrada-Peña, 1993). These results suggest a potential link between variation in CHC compounds and adaptation to environmental temperatures. The different types of off-host environments used by *I. uriae* (rock or burrow) may display significant variation in terms of temperature and relative humidity due to differential exposure to climatic factors such as sun, rain and/or snow (Buckley and Buckley, 1980). These selective factors can impact tick survival and could lead to the
quantitative variation in cuticular components observed in this study. Howard et al (1978) argued that the effectiveness of CHCs in preventing desiccation depends on the quantity of CHCs and that saturated CHCs are more effective components in protecting against water loss. By positioning temperature and humidity captors during one year in the off-host environment in a heterospecific seabird breeding colony in northern Norway (70°22' N, 31°10' E), we observed that the average temperature and relative humidity (HR) ranged respectively from -7.5 to 18°C and from 0 to 110% HR in breeding sites of Common guillemots, and from -7.5 to 10°C and from 80 to 105% HR in those occupied by puffins (data not shown). The micro-habitat used by CG ticks may therefore be more exposed to temperature and humidity variation than the more stable deep burrows used by PF ticks. However, the presence and abundance of saturated CHCs, such as 9meC<sub>24</sub> and 2meC<sub>24</sub>, were not more frequent in CG samples compared to PF samples. More detailed environmental data from each location and more sampled locations are therefore necessary to more fully evaluate this hypothesis.

Third, diet also appears to be an important factor shaping cuticular hydrocarbon profiles. In the Argentine ant *Linepithema humile*, colonies eating different prey items present particular hydrocarbon profiles that include components coming directly from the prey (Liang and Silverman, 2000). In the same way, Geiselhardt et al., (2012) showed that males of the phytophagous mustard leaf beetle *Phaedon cochleariae* preferred to mate with females reared on the same host plant compared to females from a different host plant, even though they originated from the same laboratory stock population. Studies have also shown that the ectoparasite *Varroa destructor* can mimic its host's cuticular hydrocarbons, enabling the parasite to escape the hygienic behavior of the host honeybee (Le Conte et al., 2015; Nazzi and Le Conte, 2016). In ticks, it was demonstrated that the cuticular fatty acid profiles of *Rhipicephalus* spp. presented significant differences in abundance according to host use (Shimshoni et al., 2013). In our system, avian erythrocyte membrane lipids might provide a direct pathway for the tick to synthesize cuticular hydrocarbons from pre-existing shorter-chain fatty acids and could result in the acquisition of specific hydrocarbon mixtures. This type of acquisition could explain both the variation among ticks from different host species, and variation among colony locations, if seabird diets shift among locations. Future studies could therefore be expanded to look for correlations between avian erythrocyte membrane lipid composition and tick cuticular signatures by comparing ticks and host blood samples. In addition, looking at patterns in tick cuticular lipids, such as fatty acids and steroids, which have been reported as pheromones in metastriate ticks (Sonenshine, 2004; Sonenshine et al, 1985), would also be of particular interest.

Fourth, CHCs can be altered by infection with different micro-organisms. For example, both the corn borer (*Ostrinia nubilalis*) and the common cockchafer (Melolontha melolontha) show altered CHC profiles when infected by fungi of the genus Beauveria (Lecuona et al, 1991). Moreover, Yoder and Domingus (2003) demonstrated that production of the CHCs *n*-C<sub>20</sub> and *n*-C<sub>24</sub> by *Dermacentor variabilis* ticks could protect against predatory behavior by ants. The abundances of *n*-C<sub>20</sub>, *n*-C<sub>21</sub>, *n*-C<sub>23</sub> and *n*-C<sub>24</sub> in the puffin tick pools of the present study may represent different infection statuses and/or that puffin burrows provide a less protective environment against predators, such as ants or spiders, than stones in common guillemot colonies.

In general, making inferences on environmental versus host effects will require the examination of ticks from additional sites and host types. However, regardless of the origin of host-associated differences in CHC profiles of *I. uriae*, it is important to now determine whether these differences reinforce assortative mating patterns, favoring the divergence of sympatric populations and the rapid formation of local host races. Future analyses should focus on the characterization and isolation of the main components of the cuticular mixture from ticks of each seabird host type to test for their biological activity and potential role in tick mating behavior.

## Acknowledgements

We would like to thank Yann Kolbeinsson from Northeast Iceland Nature Research Centre, Iceland for sampling help. We thank "Tiques et Maladies à Tiques" Working Group of the Société Française d'Ecologie et d'Evolution (SFE2) for stimulating discussions and support. A preprint of this article has been peer-reviewed and recommended by Peer Community In Zoology (https://doi.org/10.24072/pci.zool.100014).
## Data, scripts and codes availability

Data are available online: https://doi.org/10.5281/zenodo.6497483S. Scripts are also available online: https://doi.org/10.5281/zenodo.7018260

## Supplementary material

Supplementary material are available online: https://doi.org/10.5281/zenodo.5889077

## Conflict of interest disclosure

The authors of this preprint declare that they have no conflict of interest with the content of this article.

## Funding

Funding for this study was provided by the ANR grant ESPEVEC (ANR-13-BSV7-0018-01) to KDM. MD was supported by a fellowship from the French Ministry for National Education and Research at University of Montpellier.

## References

* Andersen, S.O., 1979. Biochemistry of Insect Cuticle. Annu. Rev. Entomol. 24, 29-59. https://doi.org/10.1146/annurev.en.24.010179.000333

* Bagnères, A.-G., Killian, A., Clement, J.-L., Lange, C., 1991. Interspecific recognition among termites of the genus *Reticulitermes*: Evidence for a role for the cuticular hydrocarbons. J. Chem. Ecol. 17, 2397-2420. https://doi.org/10.1007/BF00994590

* Bartelt, R.J., Armold, M.T., Schaner, A.M., Jackson, L.L., 1986. Comparative analysis of cuticular hydrocarbons in the *Drosophila virilis* species group. Comp. Biochem. Physiol. Part B Comp. Biochem. 83, 731-742. https://doi.org/10.1016/0305-0491(86)90138-0

* Bartelt, R.J., Jackson, L.L., 1984. Hydrocarbon Component of the *Drosophila virilis* (Diptera: Drosophilidae) Aggregation Pheromone: (Z)-10-Heneicosene. Ann. Entomol. Soc. Am. 77, 364-371. https://doi.org/10.1093/aesa/77.4.364

* Bartelt, R.J., Schaner, A.M., Jackson, L.L., 1988. Aggregation pheromones in *Drosophila borealis* and *Drosophila littoralis*. J. Chem. Ecol. 14, 1319-1327. https://doi.org/10.1007/BF01019355

* Blomquist, G.J., Bagnères, A.-G., 2010. Insect Hydrocarbons: Biology, Biochemistry, and Chemical Ecology. Cambridge University Press.

* Buckley, F.G., Buckley, P.A., 1980. Habitat Selection and Marine Birds, in: Behavior of Marine Animals. Springer, Boston, MA, pp. 69-112. https://doi.org/10.1007/978-1-4684-2988-6_3

* Desena, M.L., Edman, J.D., Clark, J.M., Symington, S.B., Scott, T.W., 1999. *Aedes aegypti* (Diptera: Culicidae) Age Determination by Cuticular Hydrocarbon Analysis of Female Legs. J. Med. Entomol. 36, 824-830. https://doi.org/10.1093/jmedent/36.6.824

* D'Ettorre, P., Mondy, N., Lenoir, A., Errard, C., 2002. Blending in with the crowd: social parasites integrate into their host colonies using a flexible chemical signature. Proc. R. Soc. Lond. B Biol. Sci. 269, 1911-1918. https://doi.org/10.1098/rspb.2002.2110

* D'Ettorre, P., Wenseleers, T., Dawson, J., Hutchinson, S., Boswell, T., Ratnieks, F.L.W., 2006. Wax combs mediate nestmate recognition by guard honeybees. Anim. Behav. 71, 773-779. https://doi.org/10.1016/j.anbehav.2005.05.014

* Dietrich, M., Beati, L., Elguero, E., Boulinier, T., McCoy, K.D., 2013. Body size and shape evolution in host races of the tick *Ixodes uriae*. Biol. J. Linn. Soc. 108, 323-334. https://doi.org/10.1111/j.1095-8312.2012.02021.x
Dietrich, M., Gómez-Díaz, E., McCoy, K.D., 2011. Worldwide Distribution and Diversity of Seabird Ticks: Implications for the Ecology and Epidemiology of Tick-Borne Pathogens. Vector-Borne Zoonotic Dis. 11, 453-470. https://doi.org/10.1089/vbz.2010.0009

Dietrich, M., Lobato, E., Boulinier, T., McCoy, K.D., 2014. An experimental test of host specialization in a ubiquitous polar ectoparasite: a role for adaptation? J. Anim. Ecol. 83, 576-587. https://doi.org/10.1111/1365-2656.12170

Espelie, K.E., Wenzel, J.W., Chang, G., 1990. Surface lipids of social wasp *Polistes melricus* say and its nest and nest pedicel and their relation to nestmate recognition. J. Chem. Ecol. 16, 2229-2241. https://doi.org/10.1007/BF01026933

Estrada-Peña, A., 1993. Climate and cuticular hydrocarbon variation in *Rhipicephalus sanguineus* ticks (Acari: Ixodidae). Parasitol. Res. 79, 512-516. https://doi.org/10.1007/BF00931594

Estrada-Peña, A., Castella, J., Morel, P.C., 1994. Cuticular Hydrocarbon Composition, Phenotypic Variability, and Geographic Relationships in Allopatric Populations of Amblyomma variegatum (Acari: Ixodidae) from Africa and the Caribbean. J. Med. Entomol. 31, 534-544. https://doi.org/10.1093/jmedent/31.4.534

Estrada-Peña, A., Estrada-Peña, R., Peiró, J.M., 1992. Differentiation of *Rhipicephalus* Ticks (Acari: Ixodidae) by Gas Chromatography of Cuticular Hydrocarbons. J. Parasitol. 78, 982. https://doi.org/10.2307/3283217

Estrada-Peña, A., Gray, J.S., Kahl, O., 1996. Variability in cuticular hydrocarbons and phenotypic discrimination of Ixodes ricinus populations (Acarina: Ixodidae) from Europe. Exp. Appl. Acarol. 20, 457-466. https://doi.org/10.1007/BF00053309

Estrada-Peña, A., Guglielmone, A.A., Mangold, A.J., Castellá, J., 1993. Patterns of cuticular hydrocarbon variation and genetic similarity between natural populations of *Amblyomma cajennense* (Acari: Ixodidae). Acta Trop. 55, 61-78. https://doi.org/10.1016/0001-706X(93)90049-H

Filshie, B.K., 1982. Fine Structure of the Cuticle of Insects and Other Arthropods, in: Insect Ultrastructure. Springer, Boston, MA, pp. 281-312. https://doi.org/10.1007/978-1-4615-7266-4_10

Geiselhardt, S., Otte, T., Hilker, M., 2012. Looking for a similar partner: host plants shape mating preferences of herbivorous insects by altering their contact pheromones. Ecol. Lett. 15, 971-977. https://doi.org/10.1111/j.1461-0248.2012.01816.x

Gibbs, A., Mousseau, T.A., Crowe, J.H., 1991. Genetic and acclimatory variation in biophysical properties of insect cuticle lipids. Proc. Natl. Acad. Sci. 88, 7257-7260. https://doi.org/10.1073/pnas.88.16.7257

Harris, M.P., Wanless, S., 2011. The Puffin. Bloomsbury Publishing.

Hervé, M., 2014. Aide-mémoire de statistique appliquée à la biologie. Constr Son Étude Anal Résultats À Aide Logiciel R Version 5.

Holze, H., Schrader, L., Buellesbach, J., 2021. Advances in deciphering the genetic basis of insect cuticular hydrocarbon biosynthesis and variation. Heredity 126, 219-234. https://doi.org/10.1038/s41437-020-00380-y

Howard, R.W., Blomquist, G.J., 2005. Ecological, Behavioral, and Biochemical Aspects of Insect Hydrocarbons. Annu. Rev. Entomol. 50, 371-393. https://doi.org/10.1146/annurev.ento.50.071803.130359

Howard, R.W., McDaniel, C.A., Blomquist, G.J., 1978. Cuticular hydrocarbons of the eastern subterranean termite, Reticulitermes flavipes (Kollar) (Isoptera: Rhinotermitidae). J. Chem. Ecol. 4, 233-245. https://doi.org/10.1007/BF00988058

Hunt, L.M., 1986. Differentiation between three species of Amblyomma ticks (Acari: Ixodidae) by analysis of cuticular hydrocarbons. Ann. Trop. Med. Parasitol. 80, 245-249. https://doi.org/10.1080/00034983.1986.11812008

Ichinose, K., Lenoir, A., 2009. Ontogeny of hydrocarbon profiles in the ant Aphaenogaster senilis and effects of social isolation. C. R. Biol. 332, 697-703. https://doi.org/10.1016/j.crvi.2009.04.002

Jallon, J.-M., David, J.R., 1987. Variation in Cuticular Hydrocarbons Among the Eight Species of the *Drosophila melanogaster* Subgroup. Evolution 41, 294-302. https://doi.org/10.2307/2409139

Kruger, E.L., Pappas, C.D., Howard, R.W., 1991. Cuticular Hydrocarbon Geographic Variation Among Seven North American Populations of *Aedes albopictus* (Diptera: Culicidae). J. Med. Entomol. 28, 859-864. https://doi.org/10.1093/jmedent/28.6.859
Kuo, T.-H., Yew, J.Y., Fedina, T.Y., Dreisewerd, K., Dierick, H.A., Pletcher, S.D., 2012. Aging modulates cuticular hydrocarbons and sexual attractiveness in *Drosophila melanogaster*. *J. Exp. Biol.* 215, 814-821. https://doi.org/10.1242/jeb.064980

Le Conte, Y, Huang, Y. Z., Roux, M., Zeng, Z.-J., Christidès, J.-P., Bagnères, A.-G., 2015. *Varroa destructor* changes its cuticular hydrocarbons to mimic new hosts. *Biol. Lett.* 11, 20150233. https://doi.org/10.1098/rsbl.2015.0233

Lecuona, R., Riba, G., Cassier, P., Clement, J.L., 1991. Alterations of insect epicuticular hydrocarbons during infection with *Beauveria bassiana* or *B. brongniartii*. *Journal of Invertebrate Pathology* 58, 10-18. https://doi.org/10.1016/0022-2011(91)90156-K

Léger, E., Vourc'h, G., Vial, L., Chevillon, C., McCoy, K.D., 2013. Changing distributions of ticks: causes and consequences. *Exp. Appl. Acarol.* 59, 219-244. https://doi.org/10.1007/s10493-012-9615-0

Leonovich, S.A., 2004. Phenol and lactone receptors in the distal sensilla of the Haller's organ in *Ixodes ricinus* ticks and their possible role in host perception. *Exp. Appl. Acarol.* 32, 89-102. https://doi.org/10.1023/B:APPA.0000018200.24760.78

Liang, D., Silverman, J., 2000. "You are what you eat": Diet modifies cuticular hydrocarbons and nestmate recognition in the Argentine ant, Linepithema humile. *Naturwissenschaften* 87, 412-416. https://doi.org/10.1007/s001140050752

Lockey, K.H., 1988. Lipids of the insect cuticle: origin, composition and function. *Comp. Biochem. Physiol. Part B Comp. Biochem.* 89, 595-645. https://doi.org/10.1016/0305-0491(88)90305-7

McCoy, K.D., Boulanger, N. (Eds.), 2015. Tiques et maladies à tiques : biologie, écologie évolutive, épidémiologie, Didactiques. IRD, Marseille.

McCoy, K.D., Boulinier, T., Tirard, C., Michalakis, Y., 2001. Host specificity of a generalist parasite: genetic evidence of sympatric host races in the seabird tick *Ixodes uriae*. *J. Evol. Biol.* 14, 395-405. https://doi.org/10.1046/j.1420-9101.2001.00290.x

McCoy, K.D., Chapuis, E., Tirard, C., Boulinier, T., Michalakis, Y., Bohec, C.L., Maho, Y.L., Gauthier-Clerc, M., 2005. Recurrent evolution of host-specialized races in a globally distributed parasite. *Proc. R. Soc. B Biol. Sci.* 272, 2389-2395. https://doi.org/10.1098/rspb.2005.3230

McCoy, K.D., Léger, E., Dietrich, M., 2013. Host specialization in ticks and transmission of tick-borne diseases: a review. *Front. Cell. Infect. Microbiol.* 3, 57. https://doi.org/10.3389/fcimb.2013.00057

McCoy, K.D., Tirard, C., 2002. Reproductive Strategies of the Seabird Tick *Ixodes uriae* (Acari: Ixodidae). *J. Parasitol.* 88, 813-816. https://doi.org/10.1645/0022-3395(2002)088[0813:RSOTST]2.0.CO;2

McCoy, K.D., Tirard, C., Michalakis, Y., 2003. Spatial genetic structure of the ectoparasite *Ixodes uriae* within breeding cliffs of its colonial seabird host. *Heredity* 91, 422-429. https://doi.org/10.1038/sj.hdy.6800339

Nazzi, F., Le Conte, Y., 2016. Ecology of *Varroa destructor*, the Major Ectoparasite of the Western Honey Bee, *Apis mellifera*. *Annu. Rev. Entomol.* 61, 417-432. https://doi.org/10.1146/annurev-ento-010715-023731

Paulmier, I., Bagnères, A.-G., Afonso, C.M., Dusticier, G., Rivière, G., Clément, J.-L., 1999. Alkenes as a sexual pheromone in the alfalfa leaf-cutter bee *Megachile rotundata*. *J. Chem. Ecol.* 25, 471-490. https://doi.org/10.1023/A:1020993518226

Peschke, K., Metzler, M., 1987. Cuticular hydrocarbons and female sex pheromones of the rove beetle, *Aleochara curtula* (Goeze) (Coleoptera:Staphylinidae). *Insect Biochem.* 17, 167-178. https://doi.org/10.1016/0020-1790(87)90157-0

Randolph, S.E., 1998. Ticks are not Insects: Consequences of Contrasting Vector Biology for Transmission Potential. *Parasitol. Today* 14, 186-192. https://doi.org/10.1016/S0169-4758(98)01224-1

Randolph, S.E., Storey, K., 1999. Impact of Microclimate on Immature Tick-Rodent Host Interactions (Acari: Ixodidae): Implications for Parasite Transmission. *J Med Entomol* 36, 741-748. https://doi.org/10.1093/jmedent/36.6.741

Rottler, A.-M., Schulz, S., Ayasse, M., 2013. Wax Lipids Signal Nest Identity in Bumblebee Colonies. *J. Chem. Ecol.* 39, 67-75. https://doi.org/10.1007/s10886-012-0229-0

Schlamp, K.K., 2005. Contact pheromone components and diel periodicity of sexual communication in peach twig borers, *Anarsia lineatella* (Lepidoptera: Gelechiidae) (PhD Thesis). *Biological Sciences Department-Simon Fraser University*.
Schöni, R., Hess, E., Blum, W., Ramstein, K., 1984. The aggregation-attachment pheromone of the tropical bont tick *Amblyomma variegatum* Fabricius (Acari, Ixodidae): Isolation, identification and action of its components. J. Insect Physiol. 30, 613–618. https://doi.org/10.1016/0022-1910(84)90045-3

Shimshoni, J.A., Erster, O., Rot, A., Cuneah, O., Soback, S., Shkap, V., 2013. Cuticular fatty acid profile analysis of three *Rhipicephalus* tick species (Acari: Ixodidae). Exp. Appl. Acarol. 61, 481–489. https://doi.org/10.1007/s10493-013-9713-7

Simmons, L.W., Thomas, M.L., Gray, B., Zuk, M., 2014. Replicated evolutionary divergence in the cuticular hydrocarbon profile of male crickets associated with the loss of song in the Hawaiian archipelago. J. Evol. Biol. 27, 2249–2257. https://doi.org/10.1111/jeb.12478

Singer, T.L., Espelie, K.E., 1996. Nest surface hydrocarbons facilitate nestmate recognition for the social wasp, *Polistes metricus* Say (Hymenoptera: Vespidae). J. Insect Behav. 9, 857–870. https://doi.org/10.1007/BF02208974

Sonenshine, D.E., 2004. Pheromones and other semiochemicals of ticks and their use in tick control. Parasitology 129, S405–S425. https://doi.org/10.1017/S003118200400486X

Sonenshine, D.E., Roe, R.M. (Eds.), 2014. Biology of ticks, 2nd ed. ed. Oxford University Press, New York.

Sonenshine, D.E., Silverstein, R.M., Brossut, R., Davis, E.E., Taylor, D., Carson, K.A., Homsher, P.J., Wang, V.B., 1985. Genital sex pheromones of ixodid ticks: 1. Evidence of occurrence in anterior reproductive tract of American dog tick, *Dermacentor variabilis* (Say) (Acari: Ixodidae). J Chem Ecol 11, 1669–1694. https://doi.org/10.1007/BF01012120

Tkachev, A.V., Dobrotvorsky, A.K., Vjalkov, A.I., Morozov, S.V., 2000. Chemical composition of lipophylic compounds from the body surface of unfed adult *Ixodes persulcatus* ticks (Acari: Ixodidae). Exp. Appl. Acarol. 24, 145–158. https://doi.org/10.1023/A:1006430323587

Uebel, E.C., Schwarz, M., Sonnet, P.E., Miller, R.W., Menzer, R.E., 1978. Evaluation of the mating stimulant pheromones of *Fannia canicularis*, *F. pusio*, and *F. femoralis* as attractants. Fla. Entomol. 139–143. https://doi.org/10.2307/3494228

Uebel, E.C., Sonnet, P.E., Bierl, B.A., Miller, R.W., 1975. Sex pheromone of the stable flyl: Isolation and preliminary identification of compounds that induce mating strike behavior. J. Chem. Ecol. 1, 377–385. https://doi.org/10.1007/BF00988839

van Zweden, J.S., d'Ettorre, P., 2010. Nestmate recognition in social insects and the role of hydrocarbons. In Insect Hydrocarbons: Biology, Biochemistry, and chemical Ecology, 222–243. http://doi.org/10.1017/CBO9780511711909.012

Wyatt, T.D., 2010. Pheromones and signature mixtures: defining species-wide signals and variable cues for identity in both invertebrates and vertebrates. J. Comp. Physiol. A 196, 685–700. https://doi.org/10.1007/s00359-010-0564-y

Yoder, J.A., Domingus, J.L., 2003. Identification of hydrocarbons that protect ticks (Acari: Ixodidae) against fire ants (Hymenoptera: Formicidae), but not lizards (Squamata: Polychrotidae), in an allomonal defense secretion. International Journal of Acarology 29, 87–91. https://doi.org/10.1080/01647950308684326