Anisakis Sensitization in the Croatian fish processing workers: Behavioral instead of occupational risk factors?

Abstract

We undertook the first study systematically evaluating the risk of Anisakis-sensitization in Croatian fish-processing workers and potential genetic susceptibility to anisakiasis. Anti-Anisakis IgE seroprevalence and risk factors for 600 employees of Croatian fish processing facilities and 466 blood donor controls, were assessed by indirect ELISA targeted with: recombinant Ani s 1 and Ani s 7 allergens, an Anisakis crude extract, the commercial ImmunoCAP kit, and questionnaires. Genetic susceptibility to anisakiasis was evaluated by genotypisation of human leukocytes alleles (HLA). Anti-Anisakis seropositive and a fraction of negative subjects were also assessed by ELISA and Western Blot (WB) for IgG seroprevalence to Trichinella spp. Overall, the observed anti-Anisakis seroprevalence inferred by indirect ELISA was significantly higher in fish processing workers (1.8%, 95% CI 0.9–3.3%) compared to the controls (0%, 0–0.8%). Seven out of 11 Ani s 1 and Ani s 7-positives and none of selected 65 negative sera, tested positive on whole-Anisakis extract (ImmunoCAP), whereas Anisakis crude extract ELISA detected 3.9% (2.4–6.0%) seropositives in fish processing workers, three (14%) of which showed IgE reactivity to milk proteins. The highest risk associated with Anisakis-sensitization among workers was fishing in the free time, rather than any of attributes related to the occupational exposure. Although no association was observed between anti-Anisakis seropositivity and wearing gloves or protective goggles, the majority of workers (92%) wore protective gloves, minimizing the risk for Anisakis sensitization via skin contact. Six HLA alleles within DRB1 gene were significantly associated with seropositivity under dominant, allelic or recessive models. All sera confirmed negative for anti-Trichinella spp. IgG. The study exhaustively covered almost all marine fish processing workers in Croatia, reflecting real-time Anisakis sensitization status within the industry, already under the influence of wide array of allergens.

About the Authors:

Ana Jerončić

Contributed equally to this work with: Ana Jerončić, Diana Nonković

Roles Conceptualization, Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing

Affiliation: University of Split, School of Medicine, Split, Croatia

ORCID logo http://orcid.org/0000-0003-1621-1956

Diana Nonković

Contributed equally to this work with: Ana Jerončić, Diana Nonković

Roles Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Teaching Institute of Public Health, County of Dalmatia, Split, Croatia

Anamarija Vrbatović

Roles Data curation, Formal analysis, Investigation, Methodology

Affiliation: Institute of Oceanography and Fisheries, Split, Croatia

Jerko Hrabar

Roles Data curation, Formal analysis, Methodology

Affiliation: Institute of Oceanography and Fisheries, Split, Croatia

Ivana Bušelić

Roles Data curation, Formal analysis, Methodology, Visualization

Affiliation: Institute of Oceanography and Fisheries, Split, Croatia

Victoria Martínez-Sernández

Roles Data curation, Formal analysis, Investigation, Methodology

Affiliation: Faculty of Pharmacy, University of Santiago de Compostela, Santiago de Compostela, Spain

Santiago A. Lojo Rocamonde

Roles Formal analysis, Methodology

Affiliation: Servicio de Análisis Clínicos, Complejo Hospitalario Universitario de Santiago de Compostela (CHUS), Santiago de Compostela, Spain

ORCID logo http://orcid.org/0000-0002-5102-0621

Florencio M. Ubeira

Roles Conceptualization, Data curation, Formal analysis, Investigation, Writing – original draft, Writing – review & editing

Affiliation: Faculty of Pharmacy, University of Santiago de Compostela, Santiago de Compostela, Spain

Sonja Jaman

Roles Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Clinical Hospital Center Split, Split, Croatia

Esma Čečuk Jeličić

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Clinical Hospital Center Split, Split, Croatia

ORCID logo http://orcid.org/0000-0002-2744-2963

Marco Amati

Roles Data curation, Formal analysis, Methodology

Affiliation: Istituto Superiore di Sanità, Roma, Italy

Maria Angeles Gomez Morales

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Istituto Superiore di Sanità, Roma, Italy

Boris Lukšić

Roles Conceptualization, Validation, Writing – original draft

Affiliation: Clinical Hospital Center Split, Split, Croatia

Ivona Mladineo

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Institute of Oceanography and Fisheries, Split, Croatia

ORCID logo http://orcid.org/0000-0003-1364-3593

Introduction

Anisakiasis is a human disease contracted by consumption of inadequately thermally processed seafood harboring infective third-stage larvae of the nematode genus Anisakis. It is an important emerging foodborne disease, recently ranked as fifth in the European risk ranking and the second of 24 foodborne parasitoses with the highest "increasing illness potential" [1]. Four specific forms are recognized; gastric, intestinal, ectopic and gastro-allergic [2], but the frequent misdiagnoses and underreporting of cases leads to the biased data interpretation and a speculative incidence. Likewise, one of the latest sources estimated 20 cases per year at European level [3], while a recently developed quantitative risk model suggested only in Spain between 7700 and 8320 of cases annually [4]. In light of these findings, Moneo et al. [5] suggested a fifth anisakiasis form important from public health prospective that encompasses sensitized asymptomatic patients, otherwise healthy individuals with a high anti-Anisakis IgE titer that potentially underwent a subclinical or undiagnosed gastric anisakiasis without allergic symptoms. Due to the lack of clinical symptoms, this form of anisakiasis, which may account for more than 12% prevalence in some populations even using component-resolved diagnosis [6,7], is mainly detected in the course of epidemiological studies [8,9]. Also, since the presence of anti-Anisakis IgE antibodies in sensitized patients is detectable in serum during many years, which in some cases could be boosted by continuous ingestion of Anisakis allergens with food [10], it is not possible to know the time of primo-infection in such individuals.

It was reported that allergic reactions to Anisakis allergens can also occur without a previous infection either by inhalation or direct contact in the domestic or occupational environment [11]. In particular, contact and air borne exposures were mainly reported in occupational settings of fishery and aquaculture workers, cooks, fishmongers and anglers [12,13]. In such cases, patients exhibit symptoms of rhino-conjunctivitis, asthma, anaphylaxis and dermatitis [14,15], which frequently lead to increased incapacity and absenteeism from work [16]. Albeit it is less frequently reported compared to orally-exposed Anisakis allergy, non-oral exposure could have a significant impact within public health domain, as according to FAO almost 60 million people are engaged in primary sector of fisheries and aquaculture [17].

Croatian fish processing industry is a valuable sector for both national and international market, with an annual revenue of 112 MU$, encompassing circa 1150 employees in freshwater, marine fish, crustacean and bivalve processing sector [18]. It is important to highlight that these data encompass both the freshwater and marine fish, crustacean and bivalve processing sectors, since no official statistics in Croatia differentiate between specific sector types. However, our estimation is that at least two-thirds of total production is effectuated by marine fish processing facilities.

In a previous study we estimated that the overall Anisakis seroprevalence in Croatian apparently healthy population is 2% [9]. Now we aimed to investigate: i) the IgE seroprevalence to Anisakis allergens (Ani s 1 and Ani s 7 allergens, Anisakis crude extract allergens, allergens included in the ImmunoCAP kit) in workers that are professionally exposed to contact, air and food-borne in the fish processing industry, and ii) the potential risk factors for Anisakis-sensitization, including genetic susceptibility to anisakiasis by genotyping of human leukocytes alleles (HLA). Complementarily, we assessed i) the risk of suffering cutaneous and/ or respiratory allergic symptoms in fish-processing workers, and ii) the IgG sensitization to another foodborne helminth Trichinella spp., in a subsample of subjects including those testing positive for anisakiasis.

Materials and methods

Ethics approval and consent to participate

This research was approved by the Ethics Committee of the Teaching Institute of Public Health, County of Split-Dalmatia, class: 500-01/14-01/1, no: 2181-103-01-14. Written informed consent was obtained from all individual adult participants included in the study.

Participants and study design

This was a multicenter, cross-sectional study encompassing all marine fish processing facilities in Croatia, conducted from December 2015 till December 2016 (Fig 1). All employed workers were invited to join the study, and 600 (78%) were sampled. Also were sampled 446 apparently healthy Croatian working-age adults not associated with the industry who were sampled from the same underlying population (controls). A STROBE checklist is provided in S1 Checklist. More details on centers and facilities are presented in S1 Table. Detailed information on the study design and sampling protocol is given in S1 Data. A sample questionnaire is provided in S2 Data, while the distribution of responses to selected questionnaire' items is shown in S1 Fig.

[Figure omitted. See PDF.]

Fig 1. Geographic map of Croatia where facilities included in the seroepidemiological study of anti-Anisakis sensitization are shown as circles with a code.

1A: Centre 1 (Brac 43°19′N 16°38′E), facility A Sardina; 2B: Centre 2 (Sinj, 43°42′N 16°38′E), facility B Conex; 2C: Centre 2, facility C Felicita; 2D: Centre 2, facility D Trenton; 3E: Centre 3 (Zadar, 44°6′51''N 15°13′40''E), facility E Sali Mardesic; 3F: Centre 3, facility F Omega Benkovac; 3G: Centre 3, facility G Mislov; 3H: Centre 3, facility H Ostrea; 3I: Centre 3, facility I Noclerius; 4J: Centre 4 (Rovinj, 45°04′59''N 13°38′402''E), facility J Mirna. The map was adapted from freely available figure on https://d-maps.com/pays.php?num_pay=319&lang=en.

https://doi.org/10.1371/journal.pntd.0008038.g001

Briefly, among the subjects who finished filling the questionnaire (n = 563; 94%), two groups were established: Group A (n = 371; 67%) that comprised workers not self-reporting allergic diseases, and Group B (n = 192; 34%) that self-reported allergic diseases. According to the tasks carried out in the facilities, non-allergic and allergic workers were then subdivided in three subgroups each, respectively (see Table 1): subgroups A1 and B1 (n = 268 and 147, respectively) included subjects working in direct contact with fish; A2 and B2 (n = 39 and 22, respectively) those working in indirect contact with fish; A3 and B3 (n = 64 and 23, respectively) subjects working in tasks with no evident contact with fish. In addition, workers included in allergic group (Group B) were further subclassified depending on the type of allergic disease they reported (Table 1) into: subgroup that included workers reporting cutaneous allergy symptoms (BA1; n = 47), workers reporting manifestations of rhino-conjunctivitis (BA2; n = 56), those reporting respiratory symptoms (BA3; n = 47), and those that included workers reporting symptoms encompassed in two or more of the above categories (BA4; n = 42).

[Figure omitted. See PDF.]

Table 1. Number of the fish-processing workers grouped by allergy symptoms, contact with the fish according to work tasks and Trisakis-170 results.

https://doi.org/10.1371/journal.pntd.0008038.t001

IgE determinations to Ani s 1 and Ani s 7 allergens

A general screening for the presence of specific anti-Anisakis IgE antibodies in the overall sampled population (n = 600) working at marine fish processing facilities was done using the indirect ELISA Trisakis-170 kit (Universidad de Santiago de Compostela, Spain). This test uses two recombinant major allergens (Ani s 1 and Ani s 7) as target to assess IgE sensitization to Anisakis sp. [9]. Interestingly, this test proved to be adequate for detection of true Anisakis infections [19] and was able to detect 100% of gastroallergic anisakiasis and 95% of cases of Anisakis sensitization associated with chronic urticaria [20]. Moreover, this method can be considered as a gold standard for specificity as the Ani s 7 recombinant sequence used in the kit has no sequence identity with any known human allergen in line with the FAO/WHO criteria (https://fermi.utmb.edu/). Sensitivity defers between recombinant allergens; Ani s 1 is less sensitive (61%; 95% confidence interval, CI 54–68%) compared to Ani s 7 (93%; 90–98%) [21].

Briefly, wells in columns 1, 4, 7, and 10 of the ELISA plates provided contain rAni s 1, wells in columns 2, 5, 8, and 11 contain a truncated sequence of rAni s 7 and the remaining columns (3, 6, 9 and 12) contain the blocking solution alone (controls). To test the sera, 100 μl of undiluted serum was added to each well and incubated for 30 min at RT with shaking (750 rpm). Then, both the washing steps and the detection of IgE antibodies were done using the reagents provided with the kit and following the kit instructions. SigmaFast OPD (Sigma-Aldrich, Madrid) was used as substrate for peroxidase. Optical densities (ODs) at 492 nm were calculated by subtracting the OD value produced by the same serum in the absence of antigen.

IgE determinations to a whole Anisakis crude extract (CE)

Anisakis CE was prepared by Dr. Cuéllar (Departamento de Parasitología, Universidad Complutense, Madrid, Spain) using A. simplex s. l. (sensu lato) L3 larvae as previously described [21]. Briefly, the wells of Greiner Bio-One ELISA plates (Greiner Bio-One España, S.A.U., Madrid, Spain) were coated with Anisakis CE at 5 μg/mL in 0.05 M carbonate-bicarbonate buffer pH 9.6 at 4ºC overnight. Afterwards, the plates containing the antigen were washed three times with TBS-T buffer (50 mM Tris-HCl, pH 7.4; 150 mM NaCl containing 0.05% Tween-20) and blocked with 1% non-fat skimmed milk in TBS (TBS-SM) during 1 h at RT. In parallel, several empty ELISA plates were also blocked with TBS-SM to serve as control. Then, the plates (with or without antigen) were aspirated, the wells filled in duplicate with 100 μL of each serum, previously diluted 1:1 in TBS-SM containing 0.05% Tween-20, and incubated for 30 min at RT with agitation at 750 rpm. Finally, the plates were washed three times with TBS-T and the IgE anti-CE revealed using the same secondary reagents and procedures included with the Trisakis-170 kit. Mean OD values for each serum were obtained subtracting the OD values obtained with the corresponding wells without antigen from the wells containing antigen.

IgE determinations by ImmunoCAP

In addition to Trisakis-170 and Anisakis CE ELISA determinations, a subsample of sera (n = 76) was also tested for the presence of anti-Anisakis IgE antibodies using a commercial test (ImmunoCAP, ThermoFisher Scientific). ImmunoCAP determinations were done automatically at the Laboratory of Immunology (Complejo Hospitalario Universitario de Santiago-CHUS, Santiago de Compostela, Spain). Chosen serum samples for ImmunoCAP determinations included all sera tested positive by Trisakis-170 (n = 11, see Results section) and randomized sera from subgroups BA1 (n = 20), subgroup BA2 (n = 18), subgroup BA3 (n = 12) and subgroup BA4 (n = 15). ImmunoCAP test shows high sensitivity (98%; 95% CI 97–100%), but poor specificity (78%; 95% CI 73–83%) [22].

IgG determinations to Trichinella antigens

In order to explore sensitization to another food-borne helminth endemic in some Croatian areas, we tested anti-Anisakis seropositive and seronegative sera for the presence of IgG antibodies to Trichinella spp. A standard protocol of indirect ELISA using Trichinella excretory and secretory antigens (ESA) was employed and results were confirmed by Western Blotting (WB) [23].

Human leukocyte allergen (HLA) genotyping of susceptibility to Anisakis

In total 208 sera; 159 (27%) of fish processing workers and 49 (11%) of controls; were HLA typed, including all 11 anti-Anisakis IgE seropositive or suspicious participants. Genomic DNA was extracted from EDTA-blood samples using High Pure PCR Template Preparation Kit (Roche Diagnostics GmbH, Germany). HLA-DRB1, HLA-DQA1 and HLA-DQB1 alleles were detected by PCR-sequence specific oligonucleotide probing method, using commercially available Immucor Lifecodes HLA-SSO typing kit (Immucor Transplant Diagnostics, Inc, Stamford, USA) and standard polymerase chain reaction sequence specific priming protocol for Olerup SSP typing kits (Olerup GmbH, Vienna, Austria).

HLA-allele and genotype frequencies were calculated by direct counting. Departure from Hardy-Weinberg equilibrium (HWE) was assessed by standard frequencies tests: chi-square, exact, and permutation tests; within the R software version 3.5 (http://www.r-project.org/). The loci HLA-DQA1*01:02, HLA-DQA1*05:05 and HLA-DQB1*03:01 deviated from HWE and were excluded from analysis. The final set included 55 HLA-alleles: 25 (46%) common, 20 (36%) low frequency and 10 (18%) rare alleles.

Haplotypes were estimated based on HLA-DRB1, HLA-DQA1, and HLA-DQB1 4-digit typing data using Expectation-Maximization algorithm within the Arlequin v3.5.2.2 software (http://cmpg.unibe.ch/software/arlequin35).

Statistical calculations

When adequate, statistic comparisons between prevalences of anti-Anisakis IgE antibodies obtained either by Anisakis CE ELISA, Trisakis-170 ELISA or by ImmunoCAP (dependent variable) and its Fisher’s confidence intervals were calculated by WinPepi 11.65 (http://www.brixtonhealth.com/pepi4windows.html). Significance and strength of potential associations between the dependent, and independent variables (gender, age, occupational contact and safety, behavioral risk factors, diet, health status) were inferred with univariate logistic regression analysis using the function logistf (logistf package, the R software version 3.5). The method employs Bayesian approach to fit a logistic regression model using Firth's bias reduction method, standardly used to reduce bias due to low event rate in odds ratio estimates [24].

Odds ratios (OR) were considered statistically significant if their 95% confidence interval (CI) did not include 1. The statistical significance was further corroborated by p-values obtained from regression model’s log-likelihood-ratio statistics, and Fisher exact test. The usage of Firth's method and a sufficient number of events per parameter in univariate logistic regression model assured reliable estimates of regression coefficients and statistical significance [25].

As genetic mechanisms underlying potential associations between HLA-alleles and seropositivity are not known, we tested these associations using allelic, dominant, and recessive models while applying statistical methods described above. The additive model was not built as it includes more than two genotype categories, which increases the number of parameters in a logistic model affecting the reliability of regression coefficient estimates. Haplotype associations with seropositivity were tested within the Arlequin v3.5.2.2 software.

Given that anti-Anisakis seropositivity is a rare event in sampled fish-processing workers and that large sample size replication studies on this population are not likely, we applied exploratory association analysis using a significance level of 0.1.

Because not all Ani s 1 and Ani s 7-sensitized workers showed to be positive to the whole-nemetode extract, statistical calculations were performed only for sera testing positive with the Trisakis-170 kit.

Results

The mean ± SD age of fish processing workers, of whom the majority were women (73%, p<0.001), was 46 ± 10 years. Women also predominated in controls (63%, p<0.001) who were, overall, younger than the processing workers (39 ± 11 years, mean difference of 6 years, 95% CI 5.0–7.6). The median age of the seropositive fish processing workers, of whom five (45%) were women, was 56 years (range 30–63) (S2 Fig).

IgE response to Ani s 1 and Ani s 7 allergens

Including the borderline serum testing positive in a first run and borderline in a second run with the Trisakis-170 kit (Tables 1 and 2), 11 out of 600 sera from fish processing workers tested positive for Ani s 7 IgE antibodies (prevalence 1.8%, Fisher’s exact 95% CI 0.9–3.3%), of which four were both Ani s 1 and Ani s 7 positive (0.7%, Fisher’s exact 95% CI 0.2–1.7%).

[Figure omitted. See PDF.]

Table 2. Selected attributes of Anisakis spp. seropositive individuals working in Croatian fish processing industry.

https://doi.org/10.1371/journal.pntd.0008038.t002

Ten of 11 positive sera were from workers reporting tasks that necessitate a direct contact with fish (subgroups B1 = 7 sera, and A1 = 3 sera). However, these differences did not reach statistical significance (OR = 4.5; 95% CI 0.6–35.1, p = 0.186). Interestingly, most of positive sera (7/11) were also from workers reporting allergies in the questionnaire: cutaneous symptoms (subgroup BA1 = 3 sera), rhino-conjunctivitis (subgroup BA2 = 2 sera) or combined symptoms (subgroup BA4 = 2 sera). Comparing the proportions of seropositive sera in the population of allergic versus non-allergic workers we obtained a significant relationship when we excluded the borderline serum from the analysis (OR = 4.6; 95% CI 1.2–18.2, p = 0.036). In case that the said serum was included, the association was significant at the 0.1 level (OR = 3.5; 95% CI 1.0–12.0, p = 0.052).

In contrast to the population of fish processing workers, all 446 control sera tested anti-Anisakis seronegative with Ani s 1 and Ani s 7 (0%, Fisher’s exact 95% CI 0–0.8%) demonstrating significantly lower prevalence than in fish processing workers (Fisher’s exact p = 0.006).

IgE response to Anisakis CE antigens

The IgE response to Anisakis CE antigens in fishing processing workers with completed questionnaires (n = 563; 37 excluded; see Table 1), was positive in 21 tested sera. Among them, 10 sera tested positive in the group B (subjects reporting allergies in the questionnaire; n = 192; 5.2%) and 11 sera were positive in group A (subjects not reporting allergic symptoms; n = 371; 3%). No significant differences (p = 0.240; Fisher’s exact test) were obtained between both groups. Also, was observed that all positive sera to CE antigens in group B corresponded to individuals working in direct contact with fish (subgroup B1; n = 146, 6.9%), while in the group A (not reporting allergies in the questionnaire) 9 subjects worked in direct contact with fish (subgroup A1; n = 268, 3.3%) and 2 subjects worked in other tasks (subgroup A3; n = 64, 3.1%). Again, no statistical differences were observed between these subgroups (p = 0.140; Fisher’s exact test). Regarding the positive sera in subgroup B1 (n = 10), 4 sera were from subjects reporting cutaneous allergic symptoms (subgroup BA1; n = 47), 3 subjects reported rhinoconjuctivitis (subgroup BA2; n = 56), 1 subject indicated respiratory symptoms (subgroup BA3; n = 47) and 2 subjects reported symptoms from several categories (subgroup BA4; n = 42). On the other hand, excluding the sera testing positive to Ani s 7 (i.e., those suffering a previous infections by Anisakis spp.) only 5/192 (2.6%, 95% CI 0.9–6.0%) sera tested positive to CE in the allergic subpopulation (subgroups BA1 = 2; BA2 = 1; BA3 = 1 and BA4 = 1). These figures were similar to those obtained for the non-allergic subpopulation for which 7/371 (1.9%, 95% CI 0.8–3.9) were positive to CE antigens (subgroups A1 = 4; A2 = 1 and A3 = 1).

As all positive sera to Anisakis CE in group B (n = 10) and 9/11 positive sera in group A fell within subgroups B1 and A1, respectively (both corresponding to workers in direct contact with fish), we analyzed whether seroprevalences in subjects working in direct contact with fish (either reporting allergic or non-allergic symptoms; subgroups B1 and A1, respectively) have a significantly higher number of sensitizations to Anisakis CE than the remaining workers (subgroups B2 + B3 + A2 + A3; n = 149). Again, the Fisher’s exact test revealed no statistically significant association (OR = 0.31; 95% CI 0.072–1.37, p = 0.120).

Anti-Anisakis IgE response measured by ImmunoCAP

In total seven out of 11 Trisakis-170 positive sera also tested positive to ImmunoCAP (sensitivity 63%, 95% CI 31–88%), whereas the results of 65 Trisakis-170 negative sera were fully concordant with ImmunoCAP (specificity 100%, 95–100%). By categories (see Table 1), the ImmunoCAP positive sera were from subgroups BA1/B1 (2 sera), BA2/B1 (2 sera), BA4/B1 (1 serum), and A1 (2 sera). ImmunoCAP values (kU/L) are included in Table 2.

Anti-Trichinella IgG reactivity

Individual anti-Anisakis positive sera by Trisakis-170, including the borderline negative serum (n = 11), and a random sample of anti-Anisakis negative sera (n = 11) were also tested for IgG antibodies to Trichinella spp. While the ELISA method detected eight anti-Trichinella positive sera, WB was negative for all of them (Table 2).

Risk factors in fish processing workers

The results of univariate analyses of associations between anti-Anisakis seropositivity and sociodemographic, behavioral, or occupational safety attributes of workers are presented in Table 3.

[Figure omitted. See PDF.]

Table 3. Sociodemographic, behavioral, and occupational safety attributes of Croatian fish processing workers (n = 600) associated with the anti-Anisakis seropositivity at the significance level of 0.1.

https://doi.org/10.1371/journal.pntd.0008038.t003

The strongest risk factor for acquiring anti-Anisakis seropositivity is fishing in the free time, followed by the traditional fish processing techniques used at home (Table 3). Contrary, no association was observed for direct contact with fish on a production line in a factory, wearing gloves or protective goggles. It should be stressed, however, that the majority of workers (92%) wore protective gloves even if not in direct contact with fish, minimizing the risk for Anisakis sensitization via skin contact.

Except for the grilling and the cooking, other fish processing techniques were rarely used by participants at home (S2 Table). Nevertheless, marinating or raw fish preparation showed strong association with the seropositivity (Table 3). No seropositive workers salted or smoked fish as it requires particular equipment, therefore the estimation of their risk was unfeasible.

Dietary habits showed that meat consumption and fish consumption less than once weekly were strong protective factors against Anisakis sensitization (Table 3). Contrary, fish consumption per se or consumption of thermally processed fish weren’t identified as significant factors. Interestingly, the type of fish or the place where fish is prepared and consumed were not associated with seropositivity, even though more than 63% participants declared consumption of the marine small oily (pelagic) fish, which traditionally is the only fish type prepared without gutting.

Male gender and older age were risk factors; men were 3.4 times more likely to become seropositive than women, and older age increased this risk by 10% annually (Table 3). A significant interaction between sex and age was shown as seropositive men were much older than seropositive women (Mann-Whitney test, p = 0.002, S2 Fig). Given the fact that men are more inclined towards fishing in spare time, after removal of those who fish in their free time, there were 1% of seropositive participants both in men (1 out of 111) and women (4 out of 391). However, as the observation is derived from just five events in total of 502 participants which is less than minimum 10, no formal testing was conducted.

In total 563 fish processing workers answered the question on chronic allergies; of whom 32 (6%) self-reported chronic allergies. The chronic allergy issue was explained to workers as the presence of any type of allergy symptoms they had before their employment (i.e. childhood and teenage), and acute allergy as an entity they have acquired since being employed in fish processing facility. Seronegative workers self-reported chronic allergies less frequently than seropositive workers (6% vs 27%; OR 0.2, 95% CI 0.04–0.6) or the controls (6% vs 14%; OR 0.4, 95% CI 0.2–0.6). When we compared the prevalence of exclusively acute allergic symptoms in workers with that of allergy symptoms in controls we found that urticaria and rhinitis were frequent in seropositive workers, less so in seronegative workers, and were present in just 2% of healthy controls (Table 4 and S3 Fig).

[Figure omitted. See PDF.]

Table 4. Prevalence of acute allergy related symptoms in seropositive and seronegative fish processing workers, and of allergy symptoms in healthy controls in Croatia.

https://doi.org/10.1371/journal.pntd.0008038.t004

Human leukocyte allergen (HLA) and Anisakis sensitization

To investigate if the Anisakis-induced allergic response might be mediated by HLA gene polymorphisms, we assessed the association between seropositivity and HLA alleles coding for different protein types within DRB1, DQA1, and DQB1 loci. We observed a high level of HLA allelic diversity (homozygotes prevalence 0–4%). The majority of alleles were not associated with the seropositivity under any genetic model tested (logistic regression, p≥0.110). However, six HLA alleles within the DRB1 gene increased the risk of seropositivity under: dominant, allelic and/or recessive models (Table 5). All of these associations were identified as strong (OR≥3.5) and were additionally supported by the results of Fisher exact tests, which were significant at the 0.1 level (Table 5), with the precision of estimates being somewhat lower. For other genes, we found similar, but fewer significant associations with only one allele associated with increased risk for acquiring seropositivity in the DQA1 gene, and three alleles identified in the DQB1 gene. Of the latter, only the allele HLA-DQB1*05:02 was allied with no IgE hypersensitivity.

[Figure omitted. See PDF.]

Table 5. Associations to anti-Anisakis seropositivity in Croatia among 10 HLA alleles from 11 seropositive (including the borderline case), 148 seronegative fish processing workers and 49 controls, at the significance level of 0.1.

https://doi.org/10.1371/journal.pntd.0008038.t005

In respect to the haplotype association, we demonstrated significant associations of anti-Anisakis seropositivity with five haplotypes, containing majority of HLA alleles listed in Table 6. The exception were HLA-DRB1*03:01, HLA-DRB1*11:01, and HLA-DQB1*05:02 loci.

[Figure omitted. See PDF.]

Table 6. Associations to anti-Anisakis seropositivity in Croatia among HLA haplotypes from 11 seropositivea, 148 seronegative fish processing workers and 49 controls.

https://doi.org/10.1371/journal.pntd.0008038.t006

Discussion

This is the first study systematically evaluating the risk of Anisakis allergen sensitization in a representative sample of 600 employees of fish processing facilities and 466 working-age blood donor controls, using highly specific and sensitive serodiagnostic test. Also tested are associations of Anisakis-sensitization in fish processing workers to a number of sociodemographic, behavioral, health-related, genetic, and professional attributes. Given that the knowledge on risk factors for Anisakis sensitization is incomplete, often contradictory and vague, and since the replication of the study is unlikely as it would have required repeated efforts to achieve formal cooperation of entire industry, we performed exploratory association analysis. Exploratory nature together with a low number of observed events might have affected the precision of some estimates and conditioned the emergence of some false-positive interactions. Nevertheless, association results presented here are mainly corroborated with the strength of associations, their consistency, biological plausibility, and/or coherence with previous knowledge, as described in detail below. More importantly, we argue that the study exhaustively covered almost all marine fish processing workers in Croatia, precisely reflecting real-time Anisakis sensitization status within the industry. This sets a baseline for the future tracking of sensitization in the sector already under the influence of wide array of allergens.

Risk factors for Anisakis sensitization in fish processing workers

Higher anti-Anisakis seroprevalences in workers (1.8%, 95% CI 0.9–3.3%) compared to controls (0%, 0–0.8%) suggests an increased risk of acquiring Anisakis seropositivity in fish processing industry. This was also postulated in France, Italy, South Africa and Spain [14,15,26–29], but the body of evidence was limited to case reports, studies with small sample size and/or non-representative sample or cross-sectional studies without control groups of individuals outside of the fish processing industry. However, in our study the highest risk for Anisakis sensitization among the workers was fishing in the free time (OR 11.4), suggesting that the professional exposure in studied environment does not represent a major risk factor for fish-processing workers when appropriate safety measures such as gloves, masks and goggles are in place. Rather, the association between IgE sensitization and fishing in free time, which is found in this subpopulation and accompanied with 5 times higher frequency of eating raw fish in amateur fisherman than in those who do not, implicate this behaviour as risky (Table 2). The fact that all seropositive workers tested by Trisakis-170 were positive to Ani s 7, which reveals true Anisakis infections [19], confirms that anti-Anisakis IgE antibodies in our study were induced by an unsuspected previous Anisakis infection rather than to an environmental exposition of workers to Anisakis allergens. Although fish workers may be more exposed to Anisakis aeroallergens than normal population, this seems not to be relevant to sensitize against major Ani s 1 and Ani s 7 secretory allergens.

While eating fish per se (OR 0.6, 95% CI 0.1–82.3, p = 0.768) or eating thermally unprocessed fish (OR 2.3, 95% CI 0.7–8.1, p = 0.162) were not risk factors in our study, eating meat and the frequency of fish consumption were. The reason for this contradiction may be that workers who didn’t eat meat actually ate fish more frequently (no vegetarians or vegans were observed, reflecting dietary habits of a typical Croatian).

An alternative explanation for the lack of association between Anisakis IgE sensitization and exposure to Anisakis allergens in the population of fish processing workers is that Ani s 1 and Ani s 7 allergen present in the Trisakis-170 kit are not adequate to measure sensitization by cutaneous and/or respiratory routes. To cover this possibility, the sera of the entire population working in the fish industry, and reporting allergy status in the questionnaire (n = 563), were tested for IgE antibodies to Anisakis CE antigen in indirect ELISA. In this case, an increased number of fish-processing employees has been found positive to Anisakis CE (3.9%, 2.4–6.0%) in comparison with Trisakis-170 (1.8%, 95% CI 0.9–3.3%). At a first glance, such difference may suggest: i) that the working population might be sensitized to other relevant Anisakis allergens present in CE and not present in the Trisakis-170 kit (Ani s 1 and Ani s 7, only), or alternatively, ii) that some allergens present in the Anisakis CE cross-reacted with other allergens to which some workers are sensitized. To distinguish between these two possibilities, the frequency of IgE sensitization to Anisakis CE was compared within the subpopulations of workers reporting versus those non-reporting allergic symptoms, as well as between subjects working in direct contact with fish with respect to the remaining subjects working in other tasks. We found no significant differences in the proportion of Anisakis IgE sensitization for any of the above comparisons, which suggests that the greater amount of positive patients to Anisakis CE, compared to Trisakis-170, is more probably due to cross-reactivity with other allergens than to sensitizations to Anisakis by skin and/or respiratory routes. Similarly, but with the lower level of evidence, Mazzucco et al. [14] showed that seroprevalence among fishermen/ sailors compared to fish industry workers was 6.7-fold greater, whereas the skin contact frequency was same among two groups. Authors unexpectedly found an inverse association between dietary fish intake and seroprevalence, and interpreted the anomaly by potentially reduced fish consumption in workers suffering allergy. Nevertheless, since a great majority of subjects in that study [14] worked in direct contact with fish (417 versus 148) a bias could be introduced and this conclusion should be taken with some caution.

Another reason for dominant behavioral instead of occupational risk factors could be that this subpopulation shows a more frequent habit of eating raw small pelagic fish, which unlike the processed anchovy and sardine, does not undergo veterinary inspection. The protective clothing (goggles, gloves, mask), being a prerequisite for the application of HACCP (Hazard Analysis of Critical Control Points) principles in Croatian fish processing industry [30], could also play a role. Workers wearing a face mask seemed to be protected (p<0.1), but because all Anisakis-sensitized workers wore gloves and only 20% of all workers wore masks, we were unable to precisely determine the risk of not wearing protective clothing in Anisakis sensitization. The effect of protective clothing and occupational allergies in general is not fully resolved. It can be insignificant in adults with atopic dermatitis [31] or represent an occupational allergen per se (natural rubber latex, [32]). Therefore, the odds of protection clothing in Anisakis allergy should be carefully considered in future in a differently designed study.

Unlike for other acute allergic symptoms in workers, the occurrence of urticaria and rhinitis (Table 4) was significantly higher in seropositive than in seronegative individuals (ORs from 4.9 to 5.0) or controls (ORs from 14.5 to 37.0). This is not surprising because urticaria is typically associated to helminths and blood-sucking arthropods with a cutaneous phase triggering Th2-type response [33]. We also found implication of other allergens in our study by differentiating between chronic and acute allergy, as seronegative workers reported more frequently all of the listed acute allergic symptoms than the controls. Similar to our finding, Jeebhay et al. [34] found that asthma and ocular-nasal symptoms are frequent in marine fish processing facilities. The implication of multiple-allergens sensitization and clinical manifestation of the allergy in case of Anisakis is still unresolved. Daschner et al. [35] observed that Anisakis-sensitized patients suffering chronic urticaria had also high prevalence of pollen, mold and dander sensitization. House dust mites sensitization was less prevalent in those patients, whom then exhibited allergic rhino-conjunctivitis and bronchial asthma, even in a region with a low burden of mites allergens.

No difference in Anisakis sensitization was observed between smoking and non-smoking workers, but ex-smokers were at significantly higher risk of sensitization than smokers, probably due to changes in the nasal mucociliary epithelium typical for ex-smokers. Namely, while the expelling of cigarette irritants through mucosal ciliature in smokers is inefficient, cigarettes induce coughing that enables the clearing process. In ex-smokers, it can take a year for the cilia to recover after cessation, during which ex-smokers can be susceptible to chest infections [36]. Whether the same mechanism has a role in Anisakis sensitization still remains to be tested.

Older male workers were at higher risk to Anisakis sensitization, contrasting a retrospective survey in French hospitals on anisakiasis incidence which demonstrated female dominance and younger average age of cases [37]. However, the latter encompassed broader range of species (genera Pseudoterranova and Anisakis) and surveyed different disease types (“esophageal, gastroduodenal, allergic, eosinophilic granuloma or another form”) in contrast to our study that focuses on sensitization. The predominance of Anisakis-sensitized women was previously found in Japan and South Korea [38,39] related to their increased consumption of raw fish, but not in USA or Croatia [9,40]. This suggests that age and gender susceptibility are strongly influenced by country-specific diversity of sociodemographic, behavioral and occupational habits of the population under study.

Seroprevalence of anti-Anisakis IgE and anti-Trichinella IgG

The Ani s 7 allergen [41] was present in all seropositive individuals, which is in agreement with previous studies that recognize it as the most relevant Anisakis allergen, including patients with Anisakis-induced chronic urticaria [20]. By contrast, Ani s 1 tested positive only in four individuals, but the inclusion of this allergen in the Trisakis-170 is of interest since IgE antibodies to Ani s 1 may be present in a few number of cases not detected by Ani s 7 [20]. The measurement of the IgE response to this allergen is also of interest since in the absence of an additional exposure to infection, IgE levels to Anis s 1 seem to persist, or even increase after sustained oral exposition to Anisakis allergens present in fishery products [10].

Interestingly, when the sensitization to the nematode whole extract presented in ImmunoCAP has been measured in a sera subsample (n = 76; fish-processing employees that were Ani s 1 and Ani s 7 seropositive or seronegative, but suffered one or more allergy symptoms), positive reaction was observed only in seven Trisakis-170-seropositive sera. This is surprising at the first glance since ImmunoCAP may contain Anisakis allergens recognized by cross-reactive antibodies [23], therefore giving false positive results. A possible explanation for these results is that Ani s 1 and/or Ani s 7 major allergens might not be well represented in the current pool of antigens used in ImmunoCAP, which may provoke false negative results when the response to these allergens is dominant. The relatively low IgE response observed with ImmunoCAP for some sera with respect to Trisakis-170 (e.g., #321 in Table 2) is in agreement with this hypothesis. Moreover, the fact that all ImmunoCAP positive sera fell within the group of Trisakis-170 positive workers reinforces our conclusion that the fish-processing environment per se is not associated with risk of Anisakis-sensitization. Otherwise, it would be expected that a significant number of workers in contact with Anisakis antigens (e.g. through ingestion of processed fish or aerosols in the working environment), tested positive by ImmunoCAP and negative by Trisakis-170.

In our previous study [9] we observed that 60% (6/10 cases) of Trisakis-170-positive sera were also positive to Ascaris sp., Toxocara sp. or both (tested by a commercial ELISA for detection of IgG, encompassing a pool of antigens with 95% specificity), suggesting that the output was merely a cross-reactivity, rather than the real status of seroprevalence. This time we tested Trichinella-sensitization for two reasons; firstly, because of the historical incidence of acute trichinellosis in the coastal area, although with low incidence (1.7–4.8/ 100,000 inhabitants in Croatia [42]); and secondly, because the migratory larvae of Trichinella spp. resembles the human infection model by L3 Anisakis. As Anisakis seropositive sera tested negative for anti-Trichinella IgG, no association was investigated. Cross-reactivity is a major problem in the serological diagnosis of parasitic infections (especially those caused by nematodes), particularly when crude parasite extracts are used, and to a lesser extent the excretory/secretory antigens (ESA). Such cross-reactivity has previously been described among A. simplex and other nematodes [43,44], and among Trichinella spp. and other zoonotic parasites [45,46]. However, studies demonstrated that WBs based on homologous ESA could distinguish between Anisakidae and Trichinella seropositivity [47] and that a WB based on T. spiralis ESA allows the definition of a unique diagnostic pattern for Trichinella-infected humans [24].

Difference in seroprevalence of healthy population estimated in Mladineo et al. [9] (2%) and in this study (0%) might have arisen from different structure of two groups. While the former recruited persons ongoing professional and systematic medical examination using stratified sampling, the latter recruited blood donors at sites close to the factories who, by definition, are frequently scrutinized for health and better reflect working population in terms of age. Compared to the sample from 2014 with the mean age of 58 years, participants from this sample were on average 19 years younger (95% CI 17–21).

HLA and Anisakis sensitization

HLA class I and class II proteins play a pivotal role in the adaptive human immune system, the latter mediating specific immunization to the antigen. In this study we identified several HLA class II loci strongly associated with Anisakis seropositivity, majority increasing the risk of acquiring seropositivity. This suggests that identified HLA-alleles might be implicated in functional improvement of MHC class II antigen-presenting role by i.e. sustained antigen presentation, in line with described hypersensitivity to Anisakis [20]. Similar to our study, Sánchez-Velasco et al. [48] reported that HLA-DRB1*04:04 increased the risk of Anisakis allergy by 4.9 times, while DQB1*02:02 was exclusive for Anisakis-allergic patients. Moreover, the epitope-binding motif of HLA molecules coded by HLA DRB1*04:04 has been evidenced at least once in the homologous segments of all the characterized allergens of A. simplex [49]. In addition to HLA-alleles already detected by Sánchez-Velasco et al. [48], we identified seven other risk HLA-alleles (Table 5). While some of these alleles might represent spurious association or be in linkage disequilibrium, several findings support their identification. Specifically, DRB1*11:01 allele is involved in presentation of group 2 allergens of Dermatophagoides spp. [50] to which all the Anisakis allergic patients previously showed sensitization [51]. This is not surprising as both Ani s 1 and 7 are encompassed within major allergens group, being products of the nematode excretory and secretory activity [5]. DQB1*05:02 that we identified as a strong protective factor against Anisakis sensitization is also a protective factor against house dust mite–sensitive allergic rhinitis [52]. The list of HLA-alleles and haplotypes we found important for Anisakis sensitization, associated with other disorders, is presented in S3 Table. It suggests that HLA association with Anisakis sensitization is primarily the result of general susceptibility of the seropositive individual to the impairment of immune system, including allergies to major allergens, rather than the nematode per se.

Conclusion

This is the largest seroepidemiological study that evaluated the risk of Anisakis sp. allergens exposure in employees of fish processing facilities, testing the relation of a number of sociodemographic, behavioural, health-related and professional traits of the scrutinised group, and the nematode seroprevalence. While we evidenced significantly higher seroprevalence in fish processing employees compared to the healthy control population, the major risk for the former group turned to be fishing in the free time. The main limitation of our study is the low number of seropositive individuals within fish processing employees, and the disequilibrium between the number of subjects working in direct contact with fish with respect to other workers, but we argue that data-analysis methodology secured reliable estimates and optimal cost/benefit ratio of the research, and that the study exhaustively covered almost all fish processing employees in Croatia, therefore reflecting accurately the contemporary Anisakis sensitisation status in the industry.

Supporting information

S1 Checklist. STROBE checklist.

https://doi.org/10.1371/journal.pntd.0008038.s001

(DOCX)

S1 Table. Demographic structure of facilities included in the study and attributes of centres encompassing target facilities.

https://doi.org/10.1371/journal.pntd.0008038.s002

(DOCX)

S2 Table. Frequency distribution of fish processing techniques that the sample of fish processing Croatian workers employs at home.

https://doi.org/10.1371/journal.pntd.0008038.s003

(DOCX)

S3 Table. The list of HLA alleles and haplotypes we found associated with Anisakis sensitisation, which have been associated with other disorders.

https://doi.org/10.1371/journal.pntd.0008038.s004

(DOCX)

S1 Data. Study design and sampling protocol of Anisakis seroprevalence study in fish processing workers.

https://doi.org/10.1371/journal.pntd.0008038.s005

(DOCX)

S2 Data. A sample questionnaire distributed to the employees in fish processing industry in Croatia enrolled in this study.

https://doi.org/10.1371/journal.pntd.0008038.s006

(DOCX)

S1 Fig. Pie charts showing responses of the employees in Croatian fish processing industry to the distributed questionnaire.

https://doi.org/10.1371/journal.pntd.0008038.s007

(TIFF)

S2 Fig. Age distribution according to seropositivity status in Croatia (Positive, Negative) and gender (M-men, W-women).

https://doi.org/10.1371/journal.pntd.0008038.s008

(DOCX)

S3 Fig. Prevalence of different acute allergy symptoms in anti-Anisakis seropositive and seronegative Croatian workers, and in controls.

https://doi.org/10.1371/journal.pntd.0008038.s009

(PDF)

Citation: Jerončić A, Nonković D, Vrbatović A, Hrabar J, Bušelić I, Martínez-Sernández V, et al. (2020) Anisakis Sensitization in the Croatian fish processing workers: Behavioral instead of occupational risk factors? PLoS Negl Trop Dis 14(1): e0008038. https://doi.org/10.1371/journal.pntd.0008038

References

1. Bouwknegt M, Devleesschauwe B, Graham H, Robertson LJ, van der Giessen J. Prioritization of foodborne parasites in Europe. Euro Surveill. 2018;23: 17–00161. https://doi.org/10.2807/1560-7917.ES.2018.23.9.17–00161https://doi.org/10.2807/1560-7917.ES.2018.23.9.17-00161

2. Nieuwenhuizen NE. Anisakis–immunology of a foodborne parasitosis. Parasite Immunol. 2016;38: 548–557. pmid:27428817

3. D’amico P, Malandra R, Costanzo F, Castigliego L, Guidi A, Gianfaldoni D, et al. Evolution of the Anisakis risk management in the European and Italian context. Food Res Int. Elsevier Ltd; 2014;64: 348–362. pmid:30011661

4. Bao M, Pierce GJ, Pascual S, González-Muñoz M, Mattiucci S, Mladineo I, et al. Assessing the risk of an emerging zoonosis of worldwide concern: anisakiasis. Sci Rep. 2017;7: 43699. pmid:28287609

5. Moneo I, Carballeda-Sangiao N, González-Muñoz M. New perspectives on the diagnosis of allergy to Anisakis spp. Curr Allergy Asthma Rep. 2017;17: 27. pmid:28429304

6. Puente P, Anadón AM, Rodero M, Romarís F, Ubeira FM, Cuéllar C. Anisakis simplex: the high prevalence in Madrid (Spain) and its relation with fish consumption. Exp Parasitol. 2008;118: 271–274. pmid:17825294

7. Ubeira FM, Anadón AM, Salgado A, Carvajal A, Ortega S, Aguirre C, López-Goikoetxea MJ, Ibanez L, Figueiras A. Synergism between prior Anisakis simplex infections and intake of NSAIDs, on the risk of upper digestive bleeding: a case-control study. PLoS Negl Trop Dis. 2011;5: e1214. pmid:21738810

8. Valiñas B, Lorenzo S, Eiras A, Figueiras A, Sanmartín ML, Ubeira FM. Prevalence of and risk factors for IgE sensitization to Anisakis simplex in a Spanish population. Allergy. 2001;56: 667–671. pmid:11421926

9. Mladineo I, Poljak V, Martínez-Sernández V, Ubeira FM. Anti-Anisakis IgE seroprevalence in the healthy Croatian coastal population and associated risk factors. PLoS Negl Trop Dis. 2014;8: e2673. pmid:24516680

10. Carballeda-Sangiao N, Rodríguez-Mahillo AI, Careche M, Navas A, Moneo I, González-Muñoz M. Changes over time in IgE sensitization to allergens of the fish parasite Anisakis spp. PLoS Negl Trop Dis. 2016;10: e0004864. pmid:27448190

11. Nieuwenhuizen NE, Lopata AL. Allergic reactions to Anisakis found in fish. Curr Allergy Asthma Rep. 2014;14: 455. pmid:25039016

12. Caffarelli C, Garrubba M, Greco C, Mastrorilli C, Povesi Dascola C. Asthma and food allergy in children: is there a connection or interaction? Front Pediatr. 2016;4: 34. pmid:27092299

13. Patelis A, Janson C, Borres MP, Nordvall L, Alving K, Malinovschi A. Aeroallergen and food IgE sensitization and local and systemic inflammation in asthma. Allergy. 2014;69: 380–387. pmid:24397423

14. Mazzucco W, Lacca G, Cusimano R, Provenzani A, Costa A, Di Noto AM, et al. Prevalence of sensitization to Anisakis simplex among professionally exposed populations in Sicily. Arch Environ Occup Heal. 2012;67: 91–97. pmid:22524649

15. Scala E, Giani M, Pirrotta L, Guerra EC, Cadoni S, Girardelli CR, et al. Occupational generalised urticaria and allergic airborne asthma due to Anisakis simplex. Eur J Dermatology. 2001;11: 249–50.

16. Jeebhay MF, Robins TG, Lopata AL. World at work: Fish processing workers. Occup Environ Med. 2004;61: 471–474. pmid:15090672

17. Food and Agriculture Organisation of the United Nations (FAO). The state of world fisheries and aquaculture. 2016. [Internet]. Available from: http://www.fao.org/3/a-i5555e.pdf.

18. Croatian Bureau of Statistics (CBS). Average monthly gross earning per person in paid employment in legal entities, according to NKD 2017. Activity sections—period, results from JOPPD forms. (in Croatian and English). 2017. Requested from https://www.dz.

19. Anadón AM, Romarís F, Escalante M, Rodríguez E, Gárate T, Cuéllar C, Ubeira FM. The Anisakis simplex Ani s 7 major allergen as an indicator of true Anisakis infections. Clin Exp Immunol. 2009;156: 471–478. pmid:19438600

20. Cuéllar C, Daschner A, Valls A, De Frutos C, Fernández-Fígares V, Anadón AM, et al. Ani s 1 and Ani s 7 recombinant allergens are able to differentiate distinct Anisakis simplex-associated allergic clinical disorders. Arch Dermatol Res. 2012;304: 283–288. pmid:22249742

21. González-Fernández J, Daschner A, Nieuwenhuizen NE, Lopata AL, Frutos CD, Valls A, Cuéllar C. Haemoglobin, a new major allergen of Anisakis simplex. Int J Parasitol. 2015;45: 399–407. Epub 2015 Feb 12. pmid:25683373

22. Anadón AM, Rodríguez E, Gárate MT, Cuéllar C, Romarís F, Chivato T, et al. Diagnosing human anisakiasis: Recombinant Ani s 1 and Ani s 7 allergens versus the UniCAP 100 fluorescence enzyme immunoassay. Clin Vaccine Immunol. 2010;17: 496–502. pmid:20107002

23. Gomez-Morales MA, Ludovisi A, Amati M, Cherchi S, Pezzotti P, Pozio E. Validation of an enzyme-linked immunosorbent assay for diagnosis of human trichinellosis. Clin Vaccine Immunol. 2008;15: 1723–1729. pmid:18827188

24. Heinze G, Schemper M. A solution to the problem of separation in logistic regression. Stat Med. 2002;21: 2409–2419. pmid:12210625

25. Vittinghoff E, McCulloch CE. Relaxing the rule of ten events per variable in logistic and cox regression. Am J Epidemiol. 2007;165: 710–718. pmid:17182981

26. Armentia A, Lombardero M, Callejo A, Martín Santos JM, Martín Gil FJ, Vega J, et al. Occupational asthma by Anisakis simplex. J Allergy Clin Immunol. 1998;102: 831–4. pmid:9819301

27. Barbuzza O, Guarneri F, Galtieri G, Gangemi S, Vaccaro M. Protein contact dermatitis and allergic asthma caused by Anisakis simplex. Contact Dermatitis. 2009;60: 239–240. pmid:19338604

28. Chord-Auger S, Miegeville M, Le Pape P. [Anisakiasis in the Nantes area. From fishmongers’ stalls to medical offices]. [In French]. Parasite. 1995;2: 395–400. pmid:8745739

29. Purello-D’Ambrosio F, Pastorello E, Gangemi S, Lombardo G, Ricciardi L, Fogliani O, et al. Incidence of sensitivity to Anisakis simplex in a risk population of fishermen/fishmongers. Ann Allergy, Asthma Immunol. 2000;84: 439–444.

30. Croatian Chamber of Commerce (HGK). Vodič za dobru higijensku praksu u industriji proizvoda ribarstva. [In Croatian] (A guide for good hygiene practices in fish processing industry). 2015. [Internet]. Available: https://www.hgk.hr/documents/vodiczadobruhigijenskupraksuuindustrijiproizvodaribarstvav2595785ef38a6ed4.pdf

31. Wei J, Gerlich J, Vogelberg C, von Mutius E, Windstetter D, Genuneit J, et al. Do young adults with atopic dermatitis avoid harmful workplace exposure at their first job? A prospective cohort study. Int Arch Occup Environ Health. 2016;89: 397–406. pmid:26188734

32. Vandenplas O, Raulf M. Occupational latex allergy: the current state of affairs. Curr Allergy Asthma Rep. 2017;17: 14. pmid:28251426

33. Kolkhir P, Balakirski G, Merk HF, Olisova O, Maurer M. Chronic spontaneous urticaria and internal parasites—a systematic review. Allergy. 2016;71: 308–322. pmid:26648083

34. Jeebhay M, Robins T, Miller M, Bateman E, Smuts M, Baatjies R, et al. Occupational allergy and asthma among salt water fish processing workers. Am J Ind Med. 2008;910: 899–910. pmid:18726880

35. Daschner A, Rodero M, De Frutos C, Valls A, Cuéllar C. Chronic urticaria is associated with a differential helminth-arthropod-related atopy phenotype. J Dermatol. 2010;37: 780–785. pmid:20883361

36. Utiyama D, Yoshida C, Goto D, Carvalho T, Santos U, Koczulla A, et al. The effects of smoking and smoking cessation on nasal mucociliary clearance, mucus properties and inflammation. Clinics. 2016;71: 344–350. pmid:27438569

37. Yera H, Fréalle É, Dutoit E, Dupouy-Camet J. A national retrospective survey of anisakidosis in France (2010–2014): decreasing incidence, female predominance, and emerging allergic potential. Parasite. 2018;25: 23. pmid:29637891

38. Sohn WM, Na BK, Kim TH, Park TJ. Anisakiasis: Report of 15 gastric cases caused by Anisakis type I larvae and a brief review of Korean anisakiasis cases. Korean J Parasitol. 2015;53: 465–470. pmid:26323845

39. De Silva D, Yamao M. A yen for sushi: an analysis of demographic and behavioural patterns of sushi consumption in Japan. J Foodserv. 2006;17: 63–76.

40. Shiferaw B, Verrill L, Booth H, Zansky SM, Norton DM, Crim S, et al. Sex-based differences in food consumption: Foodborne diseases active surveillance network (FoodNet) population survey, 2006–2007. Clin Infect Dis. 2012;54: S453–S457. pmid:22572669

41. Rodríguez E, Anadón AM, García-Bodas E, Romarís F, Iglesias R, Gárate T, et al. Novel sequences and epitopes of diagnostic value derived from the Anisakis simplex Ani s 7 major allergen. Allergy Eur J Allergy Clin Immunol. 2008;63: 219–225. pmid:18186812

42. Murrell KD, Pozio E. Worldwide occurrence and impact of human trichinellosis, 1986–2009. Emerg Infect Dis. 2011;17: 2194–2202. pmid:22172230

43. Perteguer MJ, Cuéllar C, Guillén JL, Águila C, Fenoy S, Chivato T, et al. Cross-reactivity between Anisakis simplex sensitization and visceral larva migrans by Toxocara canis. Acta Trop. 2003;89: 85–89. pmid:14636986

44. Lozano M, Martin H, Diaz S, Manas A, Valero L, Campos B. Cross-reactivity between antigens of Anisakis simplex s. l. and other ascarid nematodes. Parasite. 2004;11: 219–223. Available: pmid:15224584

45. Intapan PM, Maleewong W, Sukeepaisarnjaroen W, Morakote N. Potential use of Trichinella spiralis antigen for serodiagnosis of human capillariasis philippinensis by immunoblot analysis. Parasitol Res. 2006;98: 227–231. pmid:16333664

46. Gómez-Morales MA, Ludovisi A, Amati M, Blaga R, Zivojinovic M, Ribicich M, et al. A distinctive Western blot pattern to recognize Trichinella infections in humans and pigs. Int J Parasitol. 2012;42: 1017–1023. pmid:22995743

47. Møller LN, Krause TG, Koch A, Melbye M, Kapel CMO, Petersen E. Human antibody recognition of Anisakidae and Trichinella spp. in Greenland. Clin Microbiol Infect. 2007;13: 702–708. pmid:17484764

48. Sánchez-Velasco P, Mendizábal L, Antón EM, Ocejo-Vinyals G, Jerez J, Leyva-Cobián F. Association of hypersensitivity to the nematode Anisakis simplex with HLA class II DRB1*1502-DQB1*0601 haplotype. Hum Immunol. 2000;61: 314–319. pmid:10689122

49. Guarneri F, Guarneri C, Benvenga S. Cross-reactivity of Anisakis simplex: Possible role of Ani s 2 and Ani s 3. Int J Dermatol. 2007;46: 146–150. pmid:17269964

50. Verhoef A, Higgins JA, Thorpe CJ, Marsh SGE, Hayball JD, Lamb JR, et al. Clonal analysis of the atopic immune response to the group 2 allergen of Dermatophagoides spp.: Identification of HLA-DR and -DQ restricted T cell epitopes. Int Immunol. 1993;5: 1589–1597. pmid:7508743

51. Fæste CK, Jonscher KR, Dooper MMWB, Egge-Jacobsen W, Moen A, Daschner A, et al. Characterisation of potential novel allergens in the fish parasite Anisakis simplex. EuPA Open Proteomics.); 2014;4: 140–155. pmid:27110489

52. Zhao Y, Zhao Y, Li J, Zhang Y, Zhang L. HLA-DRB1*08:03:02 and HLA-DQB1*06:01:01 are associated with house dust mite–sensitive allergic rhinitis in Chinese subjects. Int Forum Allergy Rhinol. 2016;6: 854–861. pmid:27013183