- 1VISAVET Health Surveillance Centre, Universidad Complutense Madrid, Madrid, Spain
- 2Division of Microbiology, Department of Biotechnology and Food Science, Universidad de Burgos, Burgos, Spain
- 3Department of Food Science and Technology, University of Córdoba, Córdoba, Spain
- 4Laboratorio de Referencia e Investigación en Hepatitis Víricas, Centro Nacional de Microbiología, Instituto de Salud Carlos III, Madrid, Spain
- 5Departamento de Microbiología e Inmunología, Universidad Nacional de Río Cuarto, Córdoba, Argentina
- 6Department of Microbiology, Hospital Universitario Rio Hortega, Valladolid, Spain
- 7Laboratory of Applied Virology, Universidade Federal de Santa Catarina, Florianópolis, Brazil
- 8Department of Animal Health, Faculty of Veterinary, Universidad Complutense Madrid, Madrid, Spain
Zoonotic hepatitis E, mainly caused by hepatitis E virus (HEV) genotype (gt) 3, is a foodborne disease that has emerged in Europe in recent decades. The main animal reservoir for genotype 3 is domestic pigs. Pig liver and liver derivates are considered the major risk products, and studies focused on the presence of HEV in pig muscles are scarce. The objective of the present study was to evaluate the presence of HEV in different organs and tissues of 45 apparently healthy pigs from nine Spanish slaughterhouses (50% national production) that could enter into the food supply chain. Anti-HEV antibodies were evaluated in serum by an ELISA test. Ten samples from each animal were analyzed for the presence of HEV RNA by reverse transcription real-time PCR (RT-qPCR). The overall seroprevalence obtained was 73.3% (33/45). From the 450 samples analyzed, a total of 26 RT-qPCR positive samples were identified in the liver (7/45), feces (6/45), kidney (5/45), heart (4/45), serum (3/45), and diaphragm (1/45). This is the first report on detection of HEV RNA in kidney and heart samples of naturally infected pigs. HEV RNA detection was negative for rib, bacon, lean ham, and loin samples. These findings indicate that pig meat could be considered as a low risk material for foodborne HEV infection.
Introduction
Hepatitis E virus (HEV) is a small non-enveloped positive sense single-stranded RNA virus classified in Hepeviridae family (Emerson and Purcell, 2003) and is the main cause of viral acute hepatitis in humans worldwide (EFSA, 2017). In developing countries, the virus is mainly transmitted through contaminated water, whereas in industrialized countries, sporadic cases are basically related to animals, and hepatitis E is currently considered an emerging zoonotic disease (EFSA, 2017). In Spain, anti-IgG-VHE prevalence ranges from 0.6 to 10% in the general population (Echevarria et al., 2015), whereas it could reach up to 19% in persons exposed to pigs (Galiana et al., 2008). Generally, hepatitis E is a self-limiting disease, but it can become chronic or cause a severe disease in immunocompromised patients or with previous liver or chronic diseases (Pavio et al., 2010).
Pigs are an important zoonotic source of HEV. The swine population is considered endemic for HEV-genotype 3 (gt3) in many European countries (Pavio et al., 2010). Human cases in which foodborne route was implicated have been increasing during the last decade (EFSA, 2017). Pigs are susceptible to infection, but they do not suffer clinical disease, so visual inspections are not valid to detect possible HEV infection in the necropsy or during slaughter (Meng et al., 1997; van der Poel et al., 2001). Thus, laboratory tests are necessary to determine the potential contamination of tissues, organs, muscles, and fluids, which could enter the food chain. The development of precise strategies for the prevention and control of HEV infection should be based on advances in the knowledge of source, epidemiology, and control methods. As pork products (including meat) are highly prevalent in European food markets, it is necessary to evaluate the potential risks they represent relatively to HEV infection of humans.
Pork meat is the most widely consumed meat type in the world, with 112. 472 thousand tons consumed in 2018 (United States Department of Agriculture Foreign Agricultural Service [USDA], 2019). In Europe, Spain is the second producer after Germany, with 19% of the total pork sector production (MAPAMA, 2018). Moreover, Spain exported 2,196,648 tons of pork products in 2018 (MAPAMA, 2018). HEV-gt3 is present in swine populations in different European countries and has been linked to cases of hepatitis E in several countries (Lapa et al., 2015). In Spain, HEV has been circulating in pig populations at least since the 1980s, reaching a farm seroprevalence up to 98% [95% confidence interval (CI) 96.1–99.9%] (de Deus et al., 2008; Seminati et al., 2008; Casas et al., 2009a; Jimenez de Oya et al., 2011). HEV prevalence in Spanish domestic pig serum samples was determined at 18.8% (64/341) (Jimenez de Oya et al., 2007).
The European Food Safety Authority (EFSA) has recommended integrated studies in the food chain to determine the potential risk of HEV in pork products (EFSA, 2017). Although some studies have evaluated HEV presence in pig liver, bile, feces, or serum in abattoirs, extensive data are lacking about the HEV presence in other organs or muscles that can enter the food chain. Consequently, the objective of the present work was to investigate the presence of HEV in pig products at the moment of slaughter, in an endemic country, to determine the potential risk of pork products, especially pork meat.
Materials and Methods
Sampling Strategy
A cross-sectional study on pigs being slaughtered between November and December 2017 was undertaken through a sampling strategy with a national coverage. Nine Spanish slaughterhouses were selected according to their slaughter capacity, which represents 50% of national pig production, and were located in different regions within the country. In each slaughterhouse, five animals (between 5 and 6 months old) were randomly selected. From each animal, 10 different samples were obtained: 10 ml of blood, 100–125 g feces, and approximately 25 g of heart, kidney, liver, ribs, bacon, diaphragm, lean ham (femoral biceps), and loin head. To avoid cross contamination, sterile scalpel blades and disposable material were used for each sample, and samples were taken with the appropriate hygienic precautions. Samples were refrigerated and sent to the laboratory in less than 8 h and frozen at −80°C until processing.
Antibody Detection by ELISA
Once in the laboratory, blood was conserved at refrigeration temperature until the next day, when serum was obtained and stored at −20°C. Serum samples were tested for the presence of antibodies using the ID Screen Hepatitis E multi species indirect ELISA (IDvet, Montpellier, France) validated for swine based on recombinant gt3 capsid antigens. This ELISA kit detects IgG anti-HEV. Test procedures and interpretation of results were performed according to the manufacturer’s instructions. ELISA tests were repeated three times when the sample tested negative.
Detection of HEV by Real-Time PCR
Sample Process Control Virus
A sample process control virus (SPCV) was added to each sample immediately before the start of the analysis. The SPCV was murine norovirus 1 (MNV-1) (Diez-Valcarce et al., 2011a), which had been propagated in RAW264.7 cells to a concentration of 107 TCID50/ml, and a spike containing approximately 3 × 103 TCID50 was added to each sample.
Virus Concentration and Nucleic Acid Extraction From Pork Organs and Cuttings
The meat samples (1 cm3 from three different locations) were collected and stored in a sterile plastic bag. The extraction procedure was based on a mechanical disruption of the tissues followed by a silica-membrane-based RNA extraction (Di Bartolo et al., 2012; Rodriguez-Lazaro et al., 2012). Briefly, each sample (approximately 1 g) was finely chopped using a sterile razor blade, and then 100 mg of homogenate was transferred into a Fast Prep tube containing 200 μl phosphate-buffered saline (PBS) and 2 g of sterile 1-mm zirconia beads (BioSpec Products, Inc., Bartlesville, OK, United States). Twenty microliters of MNV-1 (∼3 × 103 TCID50) was added to each tube. The tube was then placed into a mechanical disruptor (FastPrep-24, MP Biomedicals, Santa Ana, CA, United States) and subjected to two cycles at a speed of 4 m × s–1 for 40 s. Afterward, 1 ml QIAzol (Qiagen, Hilden, Germany) was added to the tubes and vortexed 30 s, and the mixture was dispensed into a new Fast Prep tube and incubated 5 min at room temperature. At that point, 200 μl chloroform:isoamyl alcohol (24:1 Sigma Aldrich) was added, the mixture was vortexed 15 s and incubated 5 min at room temperature and then centrifuged at 10,000 × g for 15 min at 4°C. The resulting supernatant was used for immediate nucleic acid extraction using RNeasy Lipid Tissue Mini kit (Qiagen, Hilden, Germany) following manufacturer instructions, and the final 100 μl RNA extract was assayed immediately or stored at −80°C.
Virus Concentration and Nucleic Acid Isolation From Pork Feces and Serum
Two hundred fifty milligrams of samples were transferred to a 15-ml centrifuge tube and suspended in 2.25 ml of PBS containing gentamycin (10 mg/ml), and 20 μl of the SPCV (∼3 × 103 TCID50) was added to the sample. The suspension was vortexed for 60–90 s and centrifuged at 3,000 × g for 15 min. For serum analysis, 20 μl of the SPCV (∼3 × 103 TCID50) was added to 1 ml of blood, and the blood sample was centrifuged at 2,500 × g for 10 min. The supernatants from fecal or serum samples were then immediately used for nucleic acid isolation or stored at −80°C. Nucleic acids were extracted using a QIAamp viral RNA mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The final elution was performed twice with 50 μl elution buffer, resulting in a 100-μl nucleic acid extract. The nucleic acid extract was assayed immediately or stored at −80°C until analysis.
Virus Detection by RT-qPCR
The presence of the target virus (HEV) and the SPCV (MNV-1) was evaluated using reverse transcription real-time PCR (RT-qPCR). All reaction mixes included an internal amplification control (IAC), which was constructed as described by Diez-Valcarce et al., 2011a, b.
One-step duplex RT-qPCRs were performed using the oligonucleotides, controls, and conditions previously described (Diez-Valcarce et al., 2011b, 2012; Martinez-Martinez et al., 2011; Di Bartolo et al., 2012; Rodriguez-Lazaro et al., 2015). The thermocycling conditions varied slightly: 15 min at 50°C, 2 min at 95°C, followed by 45 cycles of 15 s at 95°C and 1 min at 60°C. All RT-qPCRs were conducted in a duplex format, targeting the specific viruses (HEV or MNV-1) with a FAM-labeled probe and the chimerical IAC using a VIC-labeled probe. All tests also included negative controls for viruses and for IACs.
Reporting and Interpretation of Data
For a proper interpretation of the results, four different signals were considered: (i) the target virus; (ii) the SPCV virus; (iii) the target IAC; (iv) the SPCV IAC (D’Agostino et al., 2011). When a PCR assay showed a Cq (quantification cycle) value ≤ 40, independently of the corresponding IAC Cq value, the result was interpreted as positive. When an assay showed a Cq value ≥ 40 with the corresponding IAC Cq value ≤ 40, the result was interpreted as negative. When both the target and its corresponding IAC showed Cq values ≥ 40, the reaction was considered to have failed. When at least one of the replicate HEV assays was positive, the sample was considered to be positive. In the absence of signals for SPCV and its IAC, the pre-amplification process (virus concentration and extraction steps) was concluded to have failed (D’Agostino et al., 2011). When signals for SPCV and its IAC and target IAC were present, the absence of target virus signal was conclusively considered a test negative result.
Extraction Efficiency
The extraction efficiency was calculated by comparing the Cq value of the sample containing the control (SPCV) with the Cq value of the SPCV alone, just spiked in the reagents used for concentration and extraction of the sample but without any food matrix, using the following formula: 2 (CqTNPC –Cqsample) × 100 (Diez-Valcarce et al., 2012). Efficiency results were classified as insufficient (extraction efficiency <5%), acceptable (5–25%), good (25–50%), and very good (>50%). Extraction efficiencies lower than 5% were not acceptable, and the pre-amplification process (virus concentration and extraction) of the given sample was repeated.
HEV Genotyping
Positive samples for HEV were subjected to sequence analysis, partially amplifying and sequencing ORF2, as described previously (Munoz-Chimeno et al., 2016). Phylogenetic analyses were performed with the Mega 7.0 using the method neighbor-joining with 1,000 bootstrap.
Statistical Analysis
Calculations for descriptive statistics were carried out using the WINPEPI (PEPI-for-Windows) computer programs for epidemiologist V.11.30 (Abramson, 2004). All data were compared using the χ2 test with 95% CIs, and a p-value < 0.05 was considered statistically significant. In addition, a generalized linear regression model with mixed effects (GLM) was performed considering the type of sample collected (loin head, heart, ribs, diaphragm, liver, liver exudate, lean ham, bacon, kidney, feces, and blood) and slaughterhouse (nine establishments from A to I) as explanatory variables against the binary response variable (i.e., detection of HEV by using RT-qPCR method). A backward selection method was chosen, and mean estimated parameters together with goodness-of-fit indices were obtained. The latter corresponded to the log likelihood (logL), Akaike Information Criterion (AIC), and Bayesian Information Criterion (BIC). The model structure was defined as:
which models the log odds of probability of the presence of HEV RNA by RT-qPCR (yi) as a function of a set of explanatory variables (xi, 1 … xi, p–1). β0, β1, …βp–1 are the unknown regression parameters, and σ2 the unknown (constant) error variance. The logit link function (logit π) models the log odds of the mean (π), assuming a binomial distribution of yi. The software R v.3.5.11 was used, taking as a level of significance a p-value < 0.05.
Results
Detection of Anti-HEV Antibodies in Pig Sera
Thirty-three of the 45 pigs of the study showed IgG antibodies against HEV, which represents an overall seroprevalence of 73.3% (95% CI: 58.9–84.0). Seropositive animals were found in all of the nine slaughterhouses evaluated in this study.
Efficiencies of HEV Nucleic Acid Extraction
The mean virus extraction efficiency of the process was 50.3% with a standard error of 1.83%. Values ranged from 2.26 to 98.4%. Overall, 14.6% of the samples showed acceptable extraction efficiency (5–25%), and 42.2 and 42.6% showed good (25–50%) and very good (>50%) extraction efficiencies, respectively.
Detection of HEV RNA and Distribution Between the Type of Samples and Slaughterhouses
Table 1 summarizes the results obtained for the presence of HEV RNA in the 10 different types of samples tested in this study. Significant differences in the presence of HEV RNA were obtained (p < 0.001). Although only 26 out of 450 samples (5.78%; 95% CI: 3.97–8.33%) were positive by the HEV-specific RT-qPCR, those samples came from 20 pigs; that is, 20 pigs were positive for at least one of the 10 types of analyzed samples, which represents a 44.4% (95% CI: 30.9–58.8%) of the total of pigs tested (Table 2). However, only four pigs were HEV RNA-positive in two or more samples tested (one and three pigs with four and two HEV RNA-positive samples, respectively) (Table 2). Consequently, there were 16 animals (35.6%; 95% CI: 23.3–50.2%) with only one positive sample, three (6.67%; 95% CI: 2.29–17.9%) with two positive samples and one (2.22%; 95% CI: 0.39–11.6%) with four positive samples. In three of those four animals, one of the positive samples was feces (Table 2).
Table 2. Distribution of positive samples for hepatitis E virus (HEV)-specific ELISA and RT-qPCR tests according to the type of samples.
The mean Cq values were very low (37.3 ± 0.6 SE) regardless the type of sample analyzed, ranging from 27.8 to 39.9. Interestingly, the distribution of positive samples varied according to the type of samples analyzed; whereas the muscle type samples were all negative (except for a single positive sample in the case of the diaphragm), the number of positive samples was significantly higher in the case of samples from organs (liver: n = 7, 15.5%; kidney: n = 5, 11.1%; and heart: n = 4, 8.89%) or from stool samples (n = 6, 13.3%) and serum (n = 3, 5.57%) (Table 1 and Supplementary Table 1). Regarding slaughterhouses (from A to I), all of them apart from E had HEV-positive samples. The highest prevalence corresponded to F (15.2%; 95% CI: 4.8–25.6%) with 7 out of 46 positive samples whereas for A, B, C, D, G, H, and I, average prevalence of HEV ranged from 3.9 to 7.8%. Although some variability was observed between slaughterhouses, differences were not statistically significant (p = 0.186).
HEV Genotyping
According to the mean Cq values, sequencing yield of the samples for HEV genotyping was low; however, two sequences were obtained (from liver samples of pigs 1 and 28). After phylogenetic analysis (Figure 1), Hi1 sequence was identified as genotype 3f, and Hi28 sequence showed a high identity with the KU513561 Spanish sequence, previously described in humans, which is pending of subtype assignment (Munoz-Chimeno et al., 2016).
Figure 1. Phylogenetic tree of two sequences (Hi1 and Hi28) obtained from liver samples. Built with MEGA 7.0.
Correlation Between ELISA and RT-qPCR Results
Serology by ELISA and RT-qPCR for detection of HEV RNA in blood indicated that positive results obtained by both methods were not significantly correlated (Pearson χ2 = 0.073; p = 0.793). Similarly, seven animals tested negative for HEV detection both in ELISA and in the RT-qPCR, which indicates that only a 15.5% (95% CI: 7.75–28.7%) of the pigs had not been in contact with the HEV. Similarly, 18 pigs (40.0%; 95% CI: 27.0–54.5%) tested seropositive but were negative by the HEV-specific RT-qPCR (Table 2), which highlights that although the animals were not infectious at the moment of slaughter, they had been in contact with the virus previously. Fifteen pigs were seropositive and also tested positive by RT-qPCR in at least one of the samples analyzed (33.3%; 95% CI 21.3–47.9%); HEV RNA was detected in four liver, four kidney, and four heart samples and in three feces (Table 2); 12 animals tested RNA-positive in one of the samples, and the three remaining pigs were positive in two samples (feces and liver; feces and kidney; serum and heart). Finally, five animals tested seronegative but were positive by RT-qPCR in at least one of the samples analyzed for each animal (11.1%; 95% CI: 4.84–23.5%). From those animals, two excreted HEV in feces (RNA-positive), although they tested negative by both RT-qPCR and ELISA in the rest of the samples, two tested RNA-positive only in liver or kidney, respectively, and finally one seronegative animal was RNA-positive in diaphragm, liver, fecal, and serum samples.
GLM Model for Detection of HEV RNA
The estimations obtained by the GLM model are represented in Table 3. Results from loin head samples together with those from slaughterhouse E were considered controls because all samples were negative. Considering the control group, the odds of having one HEV RNA-positive sample were calculated for the interactions between type of sample and slaughterhouse. Overall, all single factors were considered non-significant terms in the GLM model (p > 0.05). Furthermore, the odds ratio (OR) was estimated, quantifying the strength of the association between individual explanatory variables (type of sample and slaughterhouse) and the response variable (presence of HEV RNA). The strongest association was found for the kidney samples from slaughterhouse B (OR = 1.64, 95% CI: 1.11–2.44), meaning that the odds of finding one HEV RNA-positive sample was 64% higher than the control group for this combination. Other significant interactions were obtained for liver and kidney samples from slaughterhouses F, C, and G and for feces and blood samples from slaughterhouses A and F (Table 3). The odds of finding one HEV RNA-positive sample were 49% higher than the control group for these combinations.
Table 3. Estimations of the generalized linear regression model with fixed effects for the presence of hepatitis E virus (HEV) RNA by RT-qPCR as a function of the type of sample and slaughterhouse.
The goodness of fit indices AIC, logL, and BIC were estimated as −14.31, 131.32, and 403.90, respectively. Regarding model predictions, considering a cutoff value of 0.10 for the probability of having one HEV RNA-positive sample in the collected samples, 84.06% of the cases observed as negative were correctly predicted by the model, whereas 15.94% of observed negative cases were misclassified as positive (fail-safe). However, all positive cases were correctly classified by the model as such, thus indicating the high discriminatory power of the GLM model according to the studied factors.
Discussion
In this study, anti-HEV antibodies detection in serum and molecular analyses of 10 different samples (serum, feces, liver, diaphragm, kidney, heart, bacon, head of loin, rib, and lean) was performed in 45 apparently healthy pigs from different farms and collected from nine slaughterhouses geographically widespread in Spain. The overall seroprevalence obtained by an ELISA test, 73.3% (33/45), was not unexpected because a high anti-HEV antibody prevalence has been observed in an apparently healthy swine population since the 1980s in Spain (Casas et al., 2009b, 2011; Jimenez de Oya et al., 2011). This result is also in agreement with reports from a variety of European countries, as the seroprevalence described in farmed pigs ranged from 30 to 98% (Salines et al., 2017). Besides, in swine abattoirs, results are similar as for example a 59% (55.5–61.4%) of HEV seroprevalence was recently found in France (Feurer et al., 2018). Interestingly, in our study, a total of seven animals resulted negative for HEV detection by both ELISA and RT-qPCR tests, which indicates that only a maximum of 15.56% (95% CI: 7.75–28.7) was not previously in contact or infected by HEV. Not surprisingly, the presence of anti-HEV IgG in serum was higher (73.3%) than the presence of the virus detected by RT-qPCR (44.4%) in concordance with previous studies (Di Bartolo et al., 2011), and more than a half of the seropositive animals (18/33) were negative to the RT-qPCR (54.5%, 95% CI: 37.9–70.1), indicating a previous contact with HEV, but not current infection. This fact can be explained by the decrease in the prevalence of HEV RNA detection between 3 and 6 months described by previous studies (Casas et al., 2011). Five seronegative animals were positive for the RT-qPCR in different samples analyzed (41.7%, 95% CI: 19.3–68.1). Similar results were observed in the study of Di Bartolo et al. (2011) in which 4/6 seronegative animals were positive for the presence of viral RNA in bile, feces, and/or liver samples (Di Bartolo et al., 2011). Among them, two animals excreted the HEV by feces, but they tested negative in other samples, which could indicate the pass of the virus through the intestinal system after oral ingestion without any replication of the virus, as other authors have suggested (Di Bartolo et al., 2011). Besides, one animal was positive only in liver and another one only positive in the kidney. Finally, one seronegative animal was HEV RNA-positive in diaphragm, feces, liver, and serum samples, although only IgG antibodies have been detected. Absence of anti-HEV antibodies (IgG, IgM, and IgA) in pigs with HEV RNA in muscle has been described before in an experimental study of coinfection with porcine reproductive and respiratory syndrome virus (PRRSV), hypothesizing that the cause could be HEV replication in muscle cells favored by PRRSV or an interaction between heparin sulfate expressed at the surface of muscle cells with HEV particles during a long-term viremia (Salines et al., 2019a). This also could be explained by recent viral infection in which no immune response is detectable or due to a chronic infection in which antibodies may disappear because they do not persist for a long time, as other studies had demonstrated (Kanai et al., 2010).
Twenty animals tested HEV RNA-positive for at least one of the 10 samples analyzed, which indicated that in a high percentage of the pigs tested (44.4%, 95% CI: 30.9–58.8), HEV had disseminated through the organism, similar to the results previously obtained by Di Bartolo et al. (2012), with 38% (15/39) positive samples in a Spanish slaughterhouse. This is also in accordance with other authors who also reported the presence of HEV RNA in many different samples, such as lymph nodes, bladder, liver, bile, or tonsils collected from pigs in abattoirs (Leblanc et al., 2010; Raspor Lainscek et al., 2017; Feurer et al., 2018). Similarly, a total of 26 samples (5.78%) tested RT-qPCR-positive for HEV in our study: liver (n = 7; 15.5%), feces (n = 6; 13.3%), kidney (n = 5; 11.1%), heart (n = 4; 8.89%), serum (n = 3; 6.67), and diaphragm (n = 1; 2.22%) (Table 1 and Supplementary Table 1). Cq values ranged from 27.8 to 39.5, indicating a different viral load, although no association was observed between the Cq value and the type of sample (Table 2). As expected, liver was the most frequently positive sample identified with 15.6% (95% CI: 7.75–28.8), as it is the target organ for HEV replication (Lee et al., 2009). In previous studies, liver, also with the bile, had been described as the highest infected tissue (de Deus et al., 2007; Leblanc et al., 2010). In contrast, other studies describe feces as the sample with the highest prevalence for HEV presence (Di Bartolo et al., 2012; Raspor Lainscek et al., 2017). The liver positivity rate evidenced in our study (15.6%) is similar to the one obtained in Italy (20.8%) (Di Bartolo et al., 2011). However, it is higher than results obtained by other researchers in different regions, including Europe, Africa, and South America, and summarized in 2017 by Salines et al. (2017) with a mean of 5.3% (from 0.8% in Cameroon to 10% in Canada). This must be explained by the fact that some risk factors have been associated with the presence of HEV RNA in pig liver such as coinfection with PRRSV (Salines et al., 2019b), age, genetic background, or lack of hygienic measures (Walachowski et al., 2014). Nevertheless, it must be taken into account that in the present study, as in most of the published works, the presence of the HEV genetic material (HEV RNA) had been demonstrated, but the potential infectivity was not evaluated due to the difficulty of systems, which determine the viability of the virus (reliable cell culture systems or animal experimental models). However, HEV-contaminated pig liver can enter the food chain with the consequent risks demonstrated by numerous studies. In the review of Salines et al. (2017), nine different studies conducted on market pork products (raw livers, sausages, paté, etc.) were analyzed, and contamination with HEV ranged from 1 to 50% depending on the product analyzed and on the country where the survey was conducted. The highest prevalence was in products made with raw liver such as figatelli from France (Pavio et al., 2014), confirming previous studies that indicated that liver is a risk product for HEV infection, especially if it is consumed raw or undercooked, not only from pigs but also from wild boars and deer (Yazaki et al., 2003; Mizuo et al., 2005; Feagins et al., 2007). Our findings highlight that liver could be contaminated with HEV and could represent a risk for the consumer if is not well-cooked and confirms that pig liver and liver-made products must be controlled.
Similarly, six of 45 animals were shedding HEV in feces at the time of slaughterhouse, reaching a positivity of 13.3% (95% CI: 6.26–26.2), similar (6/43) to that observed in France (Leblanc et al., 2010) but lower than those observed in other countries 20–32% in Portugal (Berto et al., 2012), 21.5% in the United Kingdom (McCreary et al., 2008), 33.3% in Italy (Di Bartolo et al., 2011), or 55% in Denmark (Breum et al., 2010). It is important to highlight that as pigs are excreting the HEV in the feces, it is indispensable to optimize hygienic measures in abattoirs to avoid cross contamination with materials that could enter the food chain.
A very interesting finding in our study was that only one muscle sample (1/225; 0.44%) was HEV RNA-positive. This animal was also RT-qPCR-positive in feces, liver, and serum but was negative to the ELISA test. HEV presence in pig muscle has been demonstrated in experimental studies (Williams et al., 2001; Bouwknegt et al., 2009; Salines et al., 2019a), but few studies have been performed to establish the presence of HEV in pig muscles in naturally infected animals. Our findings are in agreement with the studies conducted in French and Canadian slaughterhouses in which the HEV seroprevalence was high, but the virus was not found in muscle samples (Leblanc et al., 2010; Feurer et al., 2018) and with a longitudinal study performed in Denmark (Krog et al., 2019). Also, one study analyzing pork products at the market failed in the detection of HEV RNA in pork chops and fresh sausages (Boxman et al., 2019). Some studies have detected HEV RNA in meat samples, but the prevalence was very low; one study performed on lingual muscle revealed a prevalence of 3 and 6% in Czechia and in Italy, respectively, whereas no positive samples were found in Spain (Di Bartolo et al., 2012). Cross contamination with feces, bile, or utensils cannot be ruled out (Di Bartolo et al., 2012).
Studies linking boar meat (Matsuda et al., 2003; Tamada et al., 2004; Li et al., 2005; Masuda et al., 2005; Wichmann et al., 2008; Rivero-Juarez et al., 2017) and also deer meat (Tei et al., 2003, 2004) with hepatitis E cases have been described since many years ago. Reports were mainly from Asia probably due to consumption habits (eating raw or undercook food). However, cases directly associated with pork meat are limited to the study of Deest et al. (2007), who reported the case of two patients who had eaten undercooked pig meat 4 weeks before suffering hepatitis E (Deest et al., 2007). Besides, in a case control study in Germany, raw pig meat and sausage were not associated with HEV cases (Wichmann et al., 2008). Although pork meat is not usually consumed raw, it is more probably that pig meat is not frequently contaminated with HEV or the viral load is too low to cause infection, confirmed by previous studies (Leblanc et al., 2010; Feurer et al., 2018) and the present study. The few studies previously performed on pig muscle have failed in the detection of HEV (Leblanc et al., 2010; Crossan et al., 2015; Krog et al., 2019).
This study is the first description of HEV RNA in kidney or heart in naturally infected pigs with five and four positive for heart and four kidney samples, respectively. These organs are not usually evaluated for the presence of HEV in pigs, and only two experimental studies have tested these organs (Williams et al., 2001; Bouwknegt et al., 2009). Extrahepatic dissemination in pigs was confirmed (Williams et al., 2001; Bouwknegt et al., 2009; Thiry et al., 2016), although the viral load was lower in those localizations (Krog et al., 2019). This fact is corroborated in the present study as the mean Cq values were higher in heart and kidney (Table 1 and Supplementary Table 1). To explain the HEV distribution to other organs or tissues other than the liver, different causes could be invoked, such as other concomitant diseases (e.g., PRRSV), which could influence the infection dynamics (Salines et al., 2019b).
Conclusion
In conclusion, our findings highlight that the Spanish pigs were frequently in contact with HEV previously (high seropositivity rate). In addition, and at the moment of slaughter, HEV could be present in pig liver, the virus could be being actively excreted (HEV RNA found in feces), and even in some cases, pigs could display viremia (HEV found in serum). Unfortunately, we only obtained two HEV sequences that limit gaining a better understanding of HEV transmission between pigs and humans. However, our results demonstrate that HEV appears to be almost always absent in different pork muscles, so pork meat could be considered as a low risk material for HEV infections via foodborne route. To confirm this hypothesis, future studies that include a larger number of animals are needed.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding authors.
Author Contributions
NG, JG, DR-L, LD, and MH designed the study. CB collected the samples at slaughterhouses. MG-B, FE, IM, AF-M, and AN performed the laboratory analysis. SG and AA designed the tables and figures. AA and MM-C performed the genotype analysis, including ORF amplification, sequencing, and phylogenetic analysis. AV, GF, and JE performed statistical analysis of the data obtained. NG and DR-L analyzed the data. NG drafted the first version of the manuscript. DR-L and SG modified and adapted the draft version, which was subsequently revised by MH and JG.
Funding
This study was partially supported by the RTA2014-00024-C04 from the Spanish Ministry of Economy and Innovation. NG and DR-L received a research grant by INTERPORC.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank all those who participated in study development, sampling, and laboratory analyses, especially the VISAVET and ITACyl teams.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02990/full#supplementary-material
Footnotes
References
Abramson, J. H. (2004). WINPEPI (PEPI-for-Windows): computer programs for epidemiologists. Epidemiol. Perspect. Innov. 1:6. doi: 10.1186/1742-5573-1-6
Berto, A., Mesquita, J. R., Hakze-van der Honing, R., Nascimento, M. S., and van der Poel, W. H. (2012). Detection and characterization of hepatitis E virus in domestic pigs of different ages in Portugal. Zoonoses Public Health 59, 477–481. doi: 10.1111/j.1863-2378.2012.01488.x
Bouwknegt, M., Rutjes, S. A., Reusken, C. B., Stockhofe-Zurwieden, N., Frankena, K., de Jong, M. C., et al. (2009). The course of hepatitis E virus infection in pigs after contact-infection and intravenous inoculation. BMC Vet. Res. 5:7. doi: 10.1186/1746-6148-5-7
Boxman, I. L. A., Jansen, C. C. C., Hagele, G., Zwartkruis-Nahuis, A., Tijsma, A. S. L., and Vennema, H. (2019). Monitoring of pork liver and meat products on the Dutch market for the presence of HEV RNA. Int. J. Food Microbiol. 296, 58–64. doi: 10.1016/j.ijfoodmicro.2019.02.018
Breum, S. O., Hjulsager, C. K., de Deus, N., Segales, J., and Larsen, L. E. (2010). Hepatitis E virus is highly prevalent in the Danish pig population. Vet. Microbiol. 146, 144–149. doi: 10.1016/j.vetmic.2010.05.002
Casas, M., Cortes, R., Pina, S., Peralta, B., Allepuz, A., Cortey, M., et al. (2011). Longitudinal study of hepatitis E virus infection in Spanish farrow-to-finish swine herds. Vet. Microbiol. 148, 27–34. doi: 10.1016/j.vetmic.2010.08.010
Casas, M., Pina, S., de Deus, N., Peralta, B., Martin, M., and Segales, J. (2009a). Pigs orally inoculated with swine hepatitis E virus are able to infect contact sentinels. Vet. Microbiol. 138, 78–84. doi: 10.1016/j.vetmic.2009.03.008
Casas, M., Pujols, J., Rosell, R., de Deus, N., Peralta, B., Pina, S., et al. (2009b). Retrospective serological study on hepatitis E infection in pigs from 1985 to 1997 in Spain. Vet. Microbiol. 135, 248–252. doi: 10.1016/j.vetmic.2008.09.075
Crossan, C., Grierson, S., Thomson, J., Ward, A., Nunez-Garcia, J., Banks, M., et al. (2015). Prevalence of hepatitis E virus in slaughter-age pigs in Scotland. Epidemiol. Infect. 143, 2237–2240. doi: 10.1017/S0950268814003100
D’Agostino, M., Cook, N., Rodriguez-Lazaro, D., and Rutjes, S. (2011). Nucleic acid amplification-based methods for detection of enteric viruses: definition of controls and interpretation of results. Food Environ. Virol. 3, 55–60. doi: 10.1007/s12560-011-9063-9068
de Deus, N., Seminati, C., Pina, S., Mateu, E., MartIn, M., and Segales, J. (2007). Detection of hepatitis E virus in liver, mesenteric lymph node, serum, bile and faeces of naturally infected pigs affected by different pathological conditions. Vet. Microbiol. 119, 105–114. doi: 10.1016/j.vetmic.2006.08.027
de Deus, N., Casas, M., Peralta, B., Nofrarias, M., Pina, S., Martin, M., et al. (2008). Hepatitis E virus infection dynamics and organic distribution in naturally infected pigs in a farrow-to-finish farm. Vet. Microbiol. 132, 19–28. doi: 10.1016/j.vetmic.2008.04.036
Deest, G., Zehner, L., Nicand, E., Gaudy-Graffin, C., Goudeau, A., and Bacq, Y. (2007). [Autochthonous hepatitis E in France and consumption of raw pig meat]. Gastroenterol. Clin. Biol. 31, 1095–1097.
Di Bartolo, I, Diez-Valcarce, M., Vasickova, P., Kralik, P., Hernandez, M., Angeloni, G., et al. (2012). Hepatitis E virus in pork production chain in Czech republic, Italy, and Spain, 2010. Emerg. Infect. Dis. 18, 1282–1289.
Di Bartolo, I., Ponterio, E., Castellini, L., Ostanello, F., and Ruggeri, F. M. (2011). Viral and antibody HEV prevalence in swine at slaughterhouse in Italy. Vet. Microbiol. 149, 330–338. doi: 10.1016/j.vetmic.2010.12.007
Diez-Valcarce, M., Cook, N., Hernandez, M., and Rodriguez-Lazaro, D. (2011a). Analytical application of a sample process control in detection of foodborne viruses. Food Anal. Methods 4, 614–618. doi: 10.1007/s12161-011-9262-9269
Diez-Valcarce, M., Kokkinos, P., Soderberg, K., Bouwknegt, M., Willems, K., de Roda-Husman, A. M., et al. (2012). Occurrence of human enteric viruses in commercial mussels at retail level in three European countries. Food Environ. Virol. 4, 73–80. doi: 10.1007/s12560-012-9078-9079
Diez-Valcarce, M., Kovac, K., Cook, N., Rodriguez-Lazaro, D., and Hernandez, M. (2011b). Construction and analytical application of Internal Amplification Controls (IAC) for detection of food supply chain-relevant viruses by real-time PCR-based assays. Food Anal. Methods 4, 437–445. doi: 10.1007/s12161-011-9224-9222
Echevarria, J. M., Fogeda, M., and Avellon, A. (2015). [Epidemiology of hepatitis E virus infection in Spain]. Enferm. Infecc. Microbiol. Clin. 33, 281–286. doi: 10.1016/j.eimc.2013.11.009
EFSA, E. B. (2017). Scientific opinion on the public health risks associated with hepatitis E virus (HEV) as a food-borne pathogen. EFSA J. 15, 4886–4889. doi: 10.2903/j.efsa.2017.4886
Emerson, S. U., and Purcell, R. H. (2003). Hepatitis E virus. Rev. Med. Virol. 13, 145–154. doi: 10.1002/rmv.384
Feagins, A. R., Opriessnig, T., Guenette, D. K., Halbur, P. G., and Meng, X. J. (2007). Detection and characterization of infectious hepatitis E virus from commercial pig livers sold in local grocery stores in the USA. J. Gen. Virol. 88(Pt 3), 912–917. doi: 10.1099/vir.0.82613-82610
Feurer, C., Le Roux, A., Rossel, R., Barnaud, E., Dumarest, M., Garry, P., et al. (2018). High load of hepatitis E viral RNA in pork livers but absence in pork muscle at French slaughterhouses. Int. J. Food Microbiol. 264, 25–30. doi: 10.1016/j.ijfoodmicro.2017.10.013
Galiana, C., Fernandez-Barredo, S., Garcia, A., Gomez, M. T., and Perez-Gracia, M. T. (2008). Occupational exposure to hepatitis E virus (HEV) in swine workers. Am. J. Trop. Med. Hyg. 78, 1012–1015. doi: 10.4269/ajtmh.2008.78.1012
Jimenez de Oya, O. N., Escribano-Romero, E., Blazquez, A. B., and Saiz, J. C. (2007). [Hepatitis E virus: zoonotic implications]. Gastroenterol. Hepatol. 30, 408–418.
Jimenez de Oya, N., de Blas, I., Blazquez, A. B., Martin-Acebes, M. A., Halaihel, N., Girones, O., et al. (2011). Widespread distribution of hepatitis E virus in Spanish pig herds. BMC Res. Notes 4:412. doi: 10.1186/1756-0500-4-412
Kanai, Y., Tsujikawa, M., Yunoki, M., Nishiyama, S., Ikuta, K., and Hagiwara, K. (2010). Long-term shedding of hepatitis E virus in the feces of pigs infected naturally, born to sows with and without maternal antibodies. J. Med. Virol. 82, 69–76. doi: 10.1002/jmv.21647
Krog, J. S., Larsen, L. E., and Breum, S. O. (2019). Tracing hepatitis E virus in pigs from birth to slaughter. Front. Vet. Sci. 6:50. doi: 10.3389/fvets.2019.00050
Lapa, D., Capobianchi, M. R., and Garbuglia, A. R. (2015). Epidemiology of hepatitis E virus in European countries. Int. J. Mol. Sci. 16, 25711–25743. doi: 10.3390/ijms161025711
Leblanc, D., Poitras, E., Gagne, M. J., Ward, P., and Houde, A. (2010). Hepatitis E virus load in swine organs and tissues at slaughterhouse determined by real-time RT-PCR. Int. J. Food Microbiol. 139, 206–209. doi: 10.1016/j.ijfoodmicro.2010.02.016
Lee, Y. H., Ha, Y., Ahn, K. K., and Chae, C. (2009). Localisation of swine hepatitis E virus in experimentally infected pigs. Vet. J. 179, 417–421. doi: 10.1016/j.tvjl.2007.10.028
Li, T. C., Chijiwa, K., Sera, N., Ishibashi, T., Etoh, Y., Shinohara, Y., et al. (2005). Hepatitis E virus transmission from wild boar meat. Emerg. Infect. Dis. 11, 1958–1960. doi: 10.3201/eid1112.051041
MAPAMA (2018). Subdirección General de Productos Ganderos, Dirección General de Producciones y Mercados Agrarios. El Sector de la Carne De Cerdo en Cifras: Principales Indicadores Económicos. Madrid: CPAGE.
Martinez-Martinez, M., Diez-Valcarce, M., Hernandez, M., and Rodriguez-Lazaro, D. (2011). Design and application of nucleic acid standards for quantitative detection of enteric viruses by real-time PCR. Food Environ. Virol. 3, 92–98. doi: 10.1007/s12560-011-9062-9069
Masuda, J., Yano, K., Tamada, Y., Takii, Y., Ito, M., Omagari, K., et al. (2005). Acute hepatitis E of a man who consumed wild boar meat prior to the onset of illness in Nagasaki. Japan. Hepatol. Res. 31, 178–183. doi: 10.1016/j.hepres.2005.01.008
Matsuda, H., Okada, K., Takahashi, K., and Mishiro, S. (2003). Severe hepatitis E virus infection after ingestion of uncooked liver from a wild boar. J. Infect. Dis. 188:944. doi: 10.1086/378074
McCreary, C., Martelli, F., Grierson, S., Ostanello, F., Nevel, A., and Banks, M. (2008). Excretion of hepatitis E virus by pigs of different ages and its presence in slurry stores in the United Kingdom. Vet. Rec. 163, 261–265. doi: 10.1136/vr.163.9.261
Meng, X. J., Purcell, R. H., Halbur, P. G., Lehman, J. R., Webb, D. M., Tsareva, T. S., et al. (1997). A novel virus in swine is closely related to the human hepatitis E virus. Proc. Natl. Acad. Sci. U.S.A. 94, 9860–9865. doi: 10.1073/pnas.94.18.9860
Mizuo, H., Yazaki, Y., Sugawara, K., Tsuda, F., Takahashi, M., Nishizawa, T., et al. (2005). Possible risk factors for the transmission of hepatitis E virus and for the severe form of hepatitis E acquired locally in Hokkaido. Japan. J. Med. Virol. 76, 341–349. doi: 10.1002/jmv.20364
Munoz-Chimeno, M., Forero, J. E., Echevarria, J. M., Munoz-Bellido, J. L., Vazquez-Lopez, L., Morago, L., et al. (2016). Full coding hepatitis E virus genotype 3 genome amplification method. J. Virol. Methods 230, 18–23. doi: 10.1016/j.jviromet.2016.01.004
Pavio, N., Meng, X. J., and Renou, C. (2010). Zoonotic hepatitis E: animal reservoirs and emerging risks. Vet. Res. 41:46. doi: 10.1051/vetres/2010018
Pavio, N., Merbah, T., and Thebault, A. (2014). Frequent hepatitis E virus contamination in food containing raw pork liver. France. Emerg. Infect. Dis. 20, 1925–1927. doi: 10.3201/eid2011.140891
Raspor Lainscek, P., Toplak, I., and Kirbis, A. (2017). A comprehensive study of hepatitis E virus infection in pigs entering a slaughterhouse in Slovenia. Vet. Microbiol. 212, 52–58. doi: 10.1016/j.vetmic.2017.11.002
Rivero-Juarez, A., Frias, M., Martinez-Peinado, A., Risalde, M. A., Rodriguez-Cano, D., Camacho, A., et al. (2017). Familial hepatitis E outbreak linked to wild boar meat consumption. Zoonoses Public Health 64, 561–565. doi: 10.1111/zph.12343
Rodriguez-Lazaro, D., Cook, N., Ruggeri, F. M., Sellwood, J., Nasser, A., Nascimento, M. S., et al. (2012). Virus hazards from food, water and other contaminated environments. FEMS Microbiol. Rev. 36, 786–814. doi: 10.1111/j.1574-6976.2011.00306.x
Rodriguez-Lazaro, D., Diez-Valcarce, M., Montes-Briones, R., Gallego, D., Hernandez, M., and Rovira, J. (2015). Presence of pathogenic enteric viruses in illegally imported meat and meat products to EU by international air travelers. Int. J. Food Microbiol. 209, 39–43. doi: 10.1016/j.ijfoodmicro.2015.04.036
Salines, M., Andraud, M., and Rose, N. (2017). From the epidemiology of hepatitis E virus (HEV) within the swine reservoir to public health risk mitigation strategies: a comprehensive review. Vet. Res. 48:31. doi: 10.1186/s13567-017-0436-433
Salines, M., Demange, A., Stephant, G., Renson, P., Bourry, O., Andraud, M., et al. (2019a). Persistent viremia and presence of hepatitis E virus RNA in pig muscle meat after experimental co-infection with porcine reproductive and respiratory syndrome virus. Int. J. Food Microbiol. 292, 144–149. doi: 10.1016/j.ijfoodmicro.2018.12.023
Salines, M., Dumarest, M., Andraud, M., Mahe, S., Barnaud, E., Cineux, M., et al. (2019b). Natural viral co-infections in pig herds affect hepatitis E virus (HEV) infection dynamics and increase the risk of contaminated livers at slaughter. Transbound. Emerg. Dis. 66, 1930–1945. doi: 10.1111/tbed.13224
Seminati, C., Mateu, E., Peralta, B., de Deus, N., and Martin, M. (2008). Distribution of hepatitis E virus infection and its prevalence in pigs on commercial farms in Spain. Vet. J. 175, 130–132. doi: 10.1016/j.tvjl.2006.11.018
Tamada, Y., Yano, K., Yatsuhashi, H., Inoue, O., Mawatari, F., and Ishibashi, H. (2004). Consumption of wild boar linked to cases of hepatitis E. J. Hepatol. 40, 869–870. doi: 10.1016/j.jhep.2003.12.026
Tei, S., Kitajima, N., Ohara, S., Inoue, Y., Miki, M., Yamatani, T., et al. (2004). Consumption of uncooked deer meat as a risk factor for hepatitis E virus infection: an age- and sex-matched case-control study. J. Med. Virol. 74, 67–70. doi: 10.1002/jmv.20147
Tei, S., Kitajima, N., Takahashi, K., and Mishiro, S. (2003). Zoonotic transmission of hepatitis E virus from deer to human beings. Lancet 362, 371–373. doi: 10.1016/S0140-6736(03)14025-14021
Thiry, D., Rose, N., Mauroy, A., Paboeuf, F., Dams, L., Roels, S., et al. (2016). Susceptibility of pigs to zoonotic hepatitis E virus genotype 3 isolated from a wild boar. Transbound. Emerg. Dis. 64, 1589–1597. doi: 10.1111/tbed.12550
United States Department of Agriculture Foreign Agricultural Service [USDA]. (2019). Livestock and Poultry: World Markets and Trade. Washington, D.C: USDA, 1–21.
van der Poel, W. H., Verschoor, F., van der Heide, R., Herrera, M. I., Vivo, A., Kooreman, M., et al. (2001). Hepatitis E virus sequences in swine related to sequences in humans, the Netherlands. Emerg. Infect. Dis. 7, 970–976. doi: 10.3201/eid0706.010608
Walachowski, S., Dorenlor, V., Lefevre, J., Lunazzi, A., Eono, F., Merbah, T., et al. (2014). Risk factors associated with the presence of hepatitis E virus in livers and seroprevalence in slaughter-age pigs: a retrospective study of 90 swine farms in France. Epidemiol. Infect. 142, 1934–1944. doi: 10.1017/S0950268813003063
Wichmann, O., Schimanski, S., Koch, J., Kohler, M., Rothe, C., Plentz, A., et al. (2008). Phylogenetic and case-control study on hepatitis E virus infection in Germany. J. Infect. Dis. 198, 1732–1741. doi: 10.1086/593211
Williams, T. P., Kasorndorkbua, C., Halbur, P. G., Haqshenas, G., Guenette, D. K., Toth, T. E., et al. (2001). Evidence of extrahepatic sites of replication of the hepatitis E virus in a swine model. J. Clin. Microbiol. 39, 3040–3046. doi: 10.1128/jcm.39.9.3040-3046.2001
Yazaki, Y., Mizuo, H., Takahashi, M., Nishizawa, T., Sasaki, N., Gotanda, Y., et al. (2003). Sporadic acute or fulminant hepatitis E in Hokkaido, Japan, may be food-borne, as suggested by the presence of hepatitis E virus in pig liver as food. J. Gen. Virol. 84(Pt 9), 2351–2357. doi: 10.1099/vir.0.19242-19240
Keywords: hepatitis E virus, prevalence, seroprevalence, pigs, pork cuttings, pig organs, food safety, slaughterhouse
Citation: García N, Hernández M, Gutierrez-Boada M, Valero A, Navarro A, Muñoz-Chimeno M, Fernández-Manzano A, Escobar FM, Martínez I, Bárcena C, González S, Avellón A, Eiros JM, Fongaro G, Domínguez L, Goyache J and Rodríguez-Lázaro D (2020) Occurrence of Hepatitis E Virus in Pigs and Pork Cuts and Organs at the Time of Slaughter, Spain, 2017. Front. Microbiol. 10:2990. doi: 10.3389/fmicb.2019.02990
Received: 16 September 2019; Accepted: 10 December 2019;
Published: 28 January 2020.
Edited by:
Maria Teresa Pérez-Gracia, Universidad CEU Cardenal Herrera, SpainReviewed by:
Petra Vasickova, Veterinary Research Institute (VRI), CzechiaPhilippe Colson, IHU Mediterranee Infection, France
Ilaria Di Bartolo, Istituto Superiore di Sanità, Italy
Copyright © 2020 García, Hernández, Gutierrez-Boada, Valero, Navarro, Muñoz-Chimeno, Fernández-Manzano, Escobar, Martínez, Bárcena, González, Avellón, Eiros, Fongaro, Domínguez, Goyache and Rodríguez-Lázaro. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nerea García, bmdhcmNpYWJAdmlzYXZldC51Y20uZXM=; David Rodríguez-Lázaro, ZHJsYXphcm9AdWJ1LmVz; cm9kbGF6ZGFAZ21haWwuY29t
†These authors have contributed equally to this work