Innate immune response to a H3N2 subtype swine influenza virus in newborn porcine trachea cells, alveolar macrophages, and precision-cut lung slices
© Delgado-Ortega et al.; licensee BioMed Central Ltd. 2014
Received: 2 December 2013
Accepted: 12 March 2014
Published: 9 April 2014
Viral respiratory diseases remain of major importance in swine breeding units. Swine influenza virus (SIV) is one of the main known contributors to infectious respiratory diseases. The innate immune response to swine influenza viruses has been assessed in many previous studies. However most of these studies were carried out in a single-cell population or directly in the live animal, in all its complexity. In the current study we report the use of a trachea epithelial cell line (newborn pig trachea cells – NPTr) in comparison with alveolar macrophages and lung slices for the characterization of innate immune response to an infection by a European SIV of the H3N2 subtype. The expression pattern of transcripts involved in the recognition of the virus, interferon type I and III responses, and the host-response regulation were assessed by quantitative PCR in response to infection. Some significant differences were observed between the three systems, notably in the expression of type III interferon mRNA. Then, results show a clear induction of JAK/STAT and MAPK signaling pathways in infected NPTr cells. Conversely, PI3K/Akt signaling pathways was not activated. The inhibition of the JAK/STAT pathway clearly reduced interferon type I and III responses and the induction of SOCS1 at the transcript level in infected NPTr cells. Similarly, the inhibition of MAPK pathway reduced viral replication and interferon response. All together, these results contribute to an increased understanding of the innate immune response to H3N2 SIV and may help identify strategies to effectively control SIV infection.
Viral respiratory diseases are still a major health issue in pigs reared under confined conditions on intensive breeding farms worldwide. Currently the most common viral pathogens are porcine reproductive and respiratory syndrome virus (PRRSV), swine influenza virus (SIV), pseudorabies virus, and porcine circovirus type 2 [1–3]. In the field, these viruses are usually found in association with each other or with bacteria such as Actinobacillus spp., Bordetella bronchiseptica, Haemophilus parasuis, Mycoplasma hyopneumoniae, Pasteurella multocida, and Streptococcus suis[1, 3–5]. In particular, influenza A viruses are a major cause of acute respiratory disease on pig breeding farms . Usually, the disease is characterized by depression, loss of appetite, abdominal breathing, tachypnea, fever, and less often, coughing. Morbidity associated with the virus is high and mortality is often very low . Influenza A viruses are classified into subtypes based on the antigenicity of their hemagglutinin (HA) and their neuraminidase (NA) surface proteins . In pigs, three subtypes are described worldwide: H1N1, H1N2, and H3N2 [7, 8], however genetic lineage may vary within each subtype depending on the geographical location (North America, Europe, and Asia). Viruses of the three subtypes have been reported frequently in European pigs, often associated with clinical disease [9, 10].
The innate immune response is crucial in the fight against respiratory viruses. It controls viral invasion and replication while the adaptive response is effective for viral clearance through the activity of adaptive immune cells such as lymphocytes [11, 12]. Among innate immune cells, epithelial cells, alveolar macrophages, and dendritic cells provide the first line of defense [11–15]. They rapidly recognize pathogens through pattern recognition receptors (PRRs) such as Toll-like receptors (TLRs), intracellular viral sensors such as retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5), and the nucleotide-binding domain, leucine-rich repeat-containing proteins (NLRs) [16, 17]. The innate immune response to SIV has been assessed in many previous studies (for review see [6, 16]). However, most of these studies were carried out in a single-cell population, sometimes isolated from other species, or directly in the pig in all its complexity making the establishment of clear conclusions difficult.
Newborn porcine trachea (NPTr) cells , porcine alveolar macrophages (PAMs), and precision-cut lung slices (PCLS) may constitute informative and complementary models to study the innate immune response to SIV. Respiratory epithelial cells and PAMs are the main target of the virus in vivo and have been successfully used for in vitro infection studies [6, 15, 16, 19]. PCLS have previously been used for infection studies in birds , cattle , and pigs [22, 23]. The PCLS culture system has several advantages over other systems: the slices can be obtained in large numbers and the general architecture of the tissue is preserved. Thus, differentiated epithelial cells, which are the main target cells of SIV, are maintained in situ and the slices remain viable for more than 7 days . Thus, NPTr cells, PAMs and PCLS together enable the study of host/pathogen interactions in single-cell type populations as well as a multi-cellular tissue, granting more accurate analysis of the contribution of epithelial cells and macrophages to the global disease response. Use of lung explants (ie PCLS) to study SIV infection has only been reported in a few cases [23–27]. However, the focus of these studies was not the characterization of the innate immune response to the virus. Here, we report the use of NPTr cells in comparison to PAMs and PCLS for the study of the innate immune response to a European H3N2 SIV. Transcripts involved in the recognition of viral patterns, interferon (IFN) type I and III responses, and regulation of the host response (suppressor of cytokine signaling, SOCS) were assessed by qPCR for their expression in response to the infection at different time points. Some significant differences were observed between the three systems, notably in the expression of type III IFN mRNA. For instance, in NPTr cells and PCLS we observed mostly IFNβ and IFNλ1 transcript expression in response to the virus while in PAMs, IFN type III was not significantly induced. The involvement of different signaling pathways such as janus kinase (JAK)/signal transducer and activator of transcription (STAT), mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK1/2), and phosphatidylinositide 3-kinase (PI3K)/protein kinase B (Akt) was also evaluated in NPTr cells. Our results show a clear induction of JAK/STAT and MAPK (ERK1/2, p38) signaling pathways while the H3N2 SIV did not induce PI3K/Akt activation in infected NPTr cells. The inhibition of the JAK/STAT pathway using JAK Inhibitor I (420099) in infected NPTR cells had a clear impact on the response of IFN types I and III and the induction of SOCS1 at the transcript level. All together, our data contribute to an increased understanding of the innate immune response of pigs to SIV.
Materials and methods
Pigs used for this research were kept in the clinic for swine and small ruminants for demonstration and veterinary student training (approval number 33.9-42502-05-09A627) or were obtained from the INRA experimental unit (Nouzilly, France). A total of 12 eight-week-old pigs (German Landrace and Large White) were used. Pigs were healthy and showed no clinical symptoms or serological evidence of influenza and other respiratory or systemic diseases. All studies were carried out in accordance with the recommendations of the European Convention for the Protection of Vertebrate Animals used for Experimental and Other Scientific Purposes (European Treaty Series, nos. 123  and 170  and in accordance with the guidelines of the Institutional Animal Care and Use committee at INRA (France). The protocol was approved by the national and local permitting authorities (animal welfare officer of the University of Veterinary Medicine, Lower Saxony State Office for Consumer Protection and Food Safety). All experimental measurements were in accordance with the requirements of the national animal welfare law. Euthanasia and tissue sampling were performed under sodium pentobarbital anesthesia, and all efforts were made to minimize suffering of the animals.
Epithelial cell line culture
The newborn pig trachea (NPTr, purchased from Istituto Zooprofilattico Sperimentale, della Lombardia e dell’Emilia Romagna, Brescia, Italy) cells  (between 30 and 50 passages) were cultured in Dulbecco’s modified Eagle medium (DMEM) (Invitrogen, Cergy Pontoise, France) supplemented with 10% fetal calf serum (FCS) (Sigma-Aldrich, Saint-Quentin, France), 20 IU/mL of penicillin, and 20 mg/mL of streptomycin (Invitrogen). Cells were plated onto sterile 24-well plastic plates (Greiner bio-one, Courtabœuf, France) and incubated at 37 °C in a 5% CO2 humidified atmosphere. Sub-passages were made when cells reached 100% confluence.
Porcine alveolar macrophages
Porcine alveolar macrophages (PAMs) were obtained from broncho-alveolar lavage (BAL) fluid from a total of 7 eight-week-old pigs and maintained in DMEM supplemented with 2% FCS, 14 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) (Gibco-Life Technologies SAS, Saint-Aubin, France), 0.1 mM non-essential amino acids (Gibco), 100 μg/mL gentamycin (Gibco), 20 IU/mL of penicillin, and 20 mg/mL of streptomycin (Invitrogen). Red blood cells were eliminated and PAMs were isolated by multiple washes and low-speed centrifugations. PAMs were plated onto sterile 24-well plastic plates (Greiner bio-one) and incubated at 37 °C and 5% CO2 for 3 h to allow macrophages to adhere. Non-adherent cells were eliminated. PAMs represented > 90% of cells in broncho-alveolar lavage fluid as previously reported .
Precision-cut lung slices
Precision-cut lung slices (PCLS) were prepared from the lungs of 5 eight-week-old pigs (minimum one slice/pig for each time point). Immediately after euthanasia, lungs were carefully removed and the left cranial, middle, and caudal lobes were filled with 37 °C low-gelling temperature agarose (GERBU Biotechnik GmbH, Gaiberg, Germany) followed by polymerization on ice. Tissue was excised in cylindrical portions (8-mm tissue-coring tool) and roughly 200 slices/pig approximately 250 μm thick were prepared by using a Krumdieck tissue slicer (model MD4000-01, TSE systems, Chesterfield, MO, USA) with a cycle speed of 60 slices/min. PCLS were incubated in 1 mL of Roswell Park Memorial Institute medium (RPMI) 1640 medium (Gibco), supplemented with 1% antibiotic/antimycotic liquid (100X Antibiotic-Antimycotic, Gibco), 1 μg/mL clotrimazole (Sigma–Aldrich), 10 μg/mL enrofloxacin (Bayer Animal Health, Germany), and 80 μg/mL kanamycin (Gibco) in a sterile 24-well plate at 37 °C and 5% CO2. The medium was changed every hour during the first 4 h and once after 24 h, prior to infection. Viability was analyzed by observing ciliary activity under a light microscope (Olympus CKX31, Tokyo, Japan). In selected samples, slices were analyzed for bronchoconstriction by addition of 10-4 M methacholine (acetyl-ß-methylcholine chloride, Sigma-Aldrich), as previously described .
Virus strain and propagation
The swine influenza virus strain A/Swine/Bissendorf/IDT1864/2003 of the H3N2 subtype was isolated from a pig with acute respiratory syndrome in a German herd and kindly provided by Ralf Dürrwald, IDT Biologika GmbH (Dessau-Rosslau, Germany). It was further propagated in the porcine epithelial NPTr cell line by infection at a low multiplicity of infection (MOI) (0.001 plaque forming unit/cell). NPTr cell line has been demonstrated suitable for the culture of SIV . Forty-eight hours post-infection, supernatants were clarified by low-speed centrifugation (2000 × g, 5 min), then stored at -80 °C until use. The viral titer of this stock reached 3.3 × 108 plaque forming units (pfu)/mL as determined by a plaque assay on NPTr cells, as described previously .
NPTr cells and PAMs were seeded onto sterile 24-well plates at 2-4 × 105 cells per well and infected with H3N2 at an MOI of 1 (enabling, according to the Poisson distribution, the infection of 63.2% of the cells with at least one single particle). After 1 h of incubation at 37 °C and 5% CO2 to allow virus adsorption, cells were washed once with phosphate buffered saline (PBS) and further maintained at 37 °C and 5% CO2 in 1 mL of medium. For PCLS infection, the procedure was nearly identical except that 106 pfu of virus/slice were used since it was not possible to determine the number of target cells in a single slice. The cells and lung slices were incubated at 37 °C and 5% CO2 for the different time points before collection for RNA extraction or staining.
Real-time PCR assays and validation of reference genes
Primer sequences, annealing temperatures of primer sets (°C), expected PCR fragment sizes (bp) and accession numbers or references
Annealing temperature (s) (°C)
PCR product (bp)
Accession number or reference
Cytokine-inducible SH2-containing protein
Hydroxymethylbilane synthase 2
Hypoxanthine phosphoribosyltransferase 1
Interleukine 1 beta
Interferon alpha (Type I)
Interferon beta (Type I)
Interferon gamma (Type II)
Interferon lambda1 (Type III)
Interferon lambda3 (Type III)
Inducible nitric oxide synthase
Myxovirus resistance 1
Myxovirus resistance 2
2′-5′-Oligoadenylate synthetase 1
SIV M protein
Protein Kinase RNA-dependent
Ribosomal protein L19
Retinoic acid-inducible gene I
Succinate dehydrogenase complex subunit A
Suppressor of cytokine signaling 1
Suppressor of cytokine signaling 3
TATA box binding protein 1
Toll Like Receptor 3
Toll Like Receptor 7
Toll Like Receptor 8
Tumor Necrosis Factor alpha
Cryosections and immunofluorescence analysis
Infected and non-infected PCLS were mounted on small pieces of filter paper with tissue-freezing medium (Jung, Heidelberg, Germany), then frozen in liquid nitrogen and kept at -80 °C prior to cutting. Ten μm-thick slices were cut by a cryotome (Reichert-Jung, Nußloch, Germany). The sections were dried overnight at room temperature and then kept frozen at -20 °C until staining.
The sections were fixed with 3% paraformaldehyde for 20 min and permeabilized with 0.2% Triton X-100 for 5 min followed by three washing steps with PBS. All antibodies were diluted in 1% bovine serum albumin (Sigma-Aldrich) and incubated with the sections for 1 h at room temperature (RT) in a humid incubation chamber. After the final incubation step, the sections were washed three times with PBS and once with distilled water. The slices were embedded in Mowiol 4-88 resin (Sigma-Aldrich), covered by no. 1½ circular micro-cover glass (12 mm) (Electron Microscopy Sciences, Hatfield, PA, USA), and stored at 4 °C until examination under the confocal microscope. For detection of infected cells, a monoclonal antibody (IgG2a) against the influenza A virus nucleoprotein (NP) (Clone AA5H, AbDSeroTec MCA400, Düsseldorf, Germany) was used at a 1:750 dilution followed by incubation with an anti-mouse IgG (Sigma-Aldrich) secondary antibody. To visualize cilia, samples were treated with a Cy3-labeled monoclonal antibody recognizing beta-tubulin (dilution 1/600) (Sigma-Aldrich). Nuclei were stained by incubating sections for 15 min at 37 °C in 4′,6′-diamidino-2-phenylindole (DAPI) (Life Technologies Inc., Darmstadt, Germany).
Antibodies used for Western blotting
Rabbit polyclonal anti-phospho-Akt (Ser473)
Rabbit monoclonal anti-phospho-p44/42 MAPK (Erk1/2)
(Thr202/Tyr204) (D13.14.4E) #4370 (Ozyme)
Rabbit polyclonal anti-phospho-JAK2 (Tyr1007/1008)
Rabbit monoclonal anti-phospho-p38 MAPK (Thr180/Tyr182)
(12 F8) #4631 (Ozyme)
Rabbit monoclonal anti-Akt (11E7) #4685 (Ozyme)
Rabbit polyclonal anti-ERK2 (GTX27948) (tebu-bio)
Rabbit monoclonal anti-JAK2 (D2E12) #3230 (Ozyme)
Rabbit polyclonal anti-p38 MAPK (GTX50566) (tebu-bio)
Antibodies and chemical inhibitors
Antibodies to phospho-Akt (Ser 473), phospho-ERK1/2 (Thr 202/Tyr 204), phospho-JAK2 (Tyr 1007/1008), and phospho-p38 (Thr180/Tyr182) were purchased from Ozyme. Monoclonal and polyclonal antibodies to Akt, ERK2, p38, and JAK2 were obtained from Tebu Bio (Le Perray-en-Yvelines, France) and Ozyme. All antibodies were used at a 1:1000 dilution in assays. Stock solutions of pharmacological inhibitors such as inhibitor U0126 (inhibiting MEK1/2 and ERK1/2), p38/SAPK2-specific inhibitor SB 202190, and JAK Inhibitor I (420099) (Millipore, Molsheim, France) were all prepared as 1000-× concentrated stocks in dimethyl sulfoxide (DMSO), in order to ensure that the final concentration of DMSO in the culture medium did not exceed 0.1%. Starting one hour before infection, NPTr cells were treated for the whole procedure with each inhibitor at a concentration of 10 μM to block the signaling. The blockage of the cascades was verified at 30 min post-infection.
Expression of mRNA in cells and tissues was expressed as relative values. All statistical analyses were done using Prism 5 computer software (Prism 5 for Windows; GraphPad Software, San Diego, CA, USA). One-Way ANOVA was used to detect differences between groups. To account for the non-normal distribution, all data were sorted by rank prior to performing the ANOVA. Tukey’s test was used to compare the means of the ranks among the groups. P values less than 0.05 were considered significant.
Reference gene selection
Viral replication in newborn pig trachea cells and alveolar macrophages
Innate immune response evaluation in epithelial cells and alveolar macrophages
Respiratory epithelial cells and PAMs are the first cells encountered by the influenza A virus during an infection in the pig. To explore the innate cellular response, the expression of various transcripts was assessed (Table 1 and Figures 2 and 3). The virus induced the expression of RIG-I mRNA as early as 3 hpi reaching the highest levels of expression at 24 h in NPTr cells (Figure 2). In PAMs, a significant increase (P < 0.001) of RIG-I mRNA was observed only after 24 h of infection (Figure 3). Regarding TLR3, TLR7, and TLR8 no statistically significant differences were observed between control and infected NPTr cells (data not shown). However in PAMs, TLR3, TLR7, and TLR8 mRNA expressions were significantly up-regulated by the viral infection (Figure 3). Expression of IFNβ mRNA was clearly increased from 3 hpi in NPTr cells (P < 0.05), and from 8 hpi in PAMs (P < 0.001) (Figures 2 and 3), and reached its highest level at 24 hpi in both cell types (P < 0.001). No significant differences in IFNα mRNA were observed between controls and infected cells, whatever the cell type (Additional file 1). For IFN type III such as IFNλ1, statistically significant increases were observed after 8 and 24 h in NPTr cells but not in PAMs (Figures 2 and 3) while no significant differences were observed for IFNλ3 (data not shown). Regarding the interferon-stimulated genes (ISGs) (Figures 2 and 3), virus-induced mRNA over-expressions (P < 0.05 and 0.001) were observed for myxovirus resistance 1 (Mx1), and Mx2 (data not shown), 2′-5′ oligoadenylate synthetase 1 (OAS1), and protein kinase R (PKR) from 3 hpi in NPTr cells and from 8 (Mx1) or 24 hpi (Mx2 – data not shown) in PAMs (P < 0.001). The pattern of expression for these genes was similar to that observed for IFNβ transcripts. Among all analyzed SOCS transcripts, only SOCS1 mRNA showed a statistically increased expression after 24 h in both NPTr cells and PAMs (P < 0.001). IL-6 (Figure 2) and IL-8 (data not shown) mRNAs were also up-regulated in response to the infection in NPTr cells at 8 and 24 hpi but not in PAMs. Regarding IL-1β and TNFα transcripts, we did not detect any significant differences between conditions in NPTr cells even though some trends of induction by the virus were observed (data not shown), while TNFα mRNA expression was significantly increased (P < 0.001) in response to the virus in PAMs after only 3 h of infection (Figure 3).
Viral infection of precision-cut lung slices
Host transcript expression in precision-cut lung slices
Lung slices offer a more realistic system than cell lines to study the interaction between host and pathogen. To better characterize the innate immune response toward the SIV strain, PCLS were infected and the expression of several host genes involved in the antiviral response was assessed at different time points using qPCR assays (Table 1 and Figure 5). An increase in RIG-I transcripts was observed at 8 h, but it did not become significant until 24 hpi (P < 0.001). We did not detect any statistically significant increase in expression of TLR3, TLR7, and TLR8 transcripts post-infection, however some trends potentially indicating a viral induction were observed (data not shown). In response to the virus, IFNβ and IFNλ1 mRNA expression were significantly increased (P < 0.001) at 24 hpi (Figure 5) while no significant increase was observed for IFNα and IFNλ3 (Additional file 1 and data not shown). Along with the increased expression of IFN types I and III transcripts, mRNA expression of ISGs Mx1, Mx2, and PKR was also significantly elevated relative to controls after 24 h of infection (P < 0.001) (Figure 5 and data not shown). The mRNA expression of OAS1 was significantly increased (P < 0.05) after only 8 h of infection (Figure 5). Unlike IL-1β, IL-8, and TNFα, IL-6 transcription was significantly up-regulated 8 h after infection (Figure 5). Overall, except for IFNλ1 and IL-6, similar patterns of transcript expression were observed in vitro in the epithelial cells and alveolar macrophages and ex vivo in the PCLS. Regarding the mRNA expression of SOCS (CISH, SOCS1, and SOCS3), only SOCS1 transcripts were significantly increased at 24 hpi in PCLS (P < 0.01) (Figure 5). These findings are similar to those observed with NPTr cells and PAMs and suggest an involvement of SOCS1 in the immune response against SIV.
Activation of MAP kinases and JAK/STAT pathways in virus-infected epithelial cells
Specific inhibition of pathways activated in epithelial cells
The innate immune response to swine influenza viruses has been assessed in many previous studies [6, 16]. However, to our knowledge most of these studies were carried out in single-cell populations or directly in an animal host. In this study we compared an epithelial cell line (NPTr), alveolar macrophages (PAMs), and precision cut lung slices (PCLS) for their innate immune response to strain A/Swine/Bissendorf/IDT1864/2003 (H3N2). Then, we assessed several signaling pathways and their relation to the innate immune response observed in the NPTr cells.
Viral replication was higher in NPTr cells than in PCLS and PAMs. This is in agreement with the current knowledge of influenza pathogenesis. Indeed, it is well known that the main target of the virus is epithelial cells [12, 15] and NPTr is a cell line consisting of tracheal epithelial cells . This was also confirmed by the immune-staining of the PCLS, where many ciliated epithelial cells were observed in the associated bronchi. In the three systems used in this study, a strong and progressive innate immune response was triggered by the virus. The response in NPTr cells and PCLS was very similar except for a trend toward TLR induction in the latter. This could be explained by the presence of various cell types such as macrophages and/or dendritic cells in this tissue [6, 14]. Moreover, observed discrepancies between NPTr cells and PCLS may also be related to the less controlled percentage of cells infected in the PCLS (MOI cannot be determined) and to the potential resistance of a significant proportion of the various cell types to SIV infection. Our study further demonstrated the value of using PCLS in the study of pig/SIV interactions. RIG-I transcription was induced by the virus in NPTr cells, PAMs, and PCLS while TLR3, TLR7 as well as TLR8 transcription were significantly induced only in PAMs. This result was a little bit surprising since it has been observed recently that airway epithelial cells mostly use surface and endosomal TLR3 to detect influenza virus and to initiate IFN production . Following recognition of the virus, interferon responses were triggered. In NPTr cells and PCLS we observed mostly IFNβ and IFNλ1 transcript expression in response to the virus while in PAMs, IFN type III was not significantly induced. Regarding the expression of IFNα transcripts, no significant differences between infected and non-infected cells/tissues were observed. This result is surprising as IFNα transcripts were previously detected in PAMs infected by another strain of SIV . In that study IFNα mRNA was detected earlier than IFNβ mRNA. The use of qPCR primers targeting IFNα 1-17 may have masked differences in the expression of transcripts associated to a specific IFNα gene. It is known that significant variations in the expression of the various IFNα transcripts upon infection can be observed . Recently, it has been demonstrated that type III IFN is the predominant IFN produced by the airway epithelium . Our data showed the increased production of mRNAs encoding both type I and type III IFNs in epithelial cells. The absence of type III IFN induction in PAMs supports the idea of a predominant role of type III IFN in a specific cell type such as epithelial cells as suggested by Ioannidis et al. . In response to IFN, ISG (Mx1, Mx2, OAS1, PKR) transcripts were identified in the three systems. This observation was particularly clear with the NPTr cells probably because of their nature (epithelial versus alveolar macrophages) and their high purity in comparison to primary epithelial cells located in the slices. These two features also probably account for the quicker and higher response observed in NPTr cells when compared to PAMs and PCLs. Because of their role in the regulation of immune response [47, 48] and in influenza pathogenesis , we looked at the expression of SOCS transcripts in response to SIV in the three systems. Twenty-four hpi SOCS1 transcripts were up-regulated in the three systems. This observation confirmed our previous study postulating the inducible nature of SOCS1 in porcine lungs . Instead of SOCS1, a role for SOCS3 has been recently identified in pigs . The authors showed that the lower susceptibility of pigs to contemporary Eurasian highly pathogenic avian influenza (HPAI) H5N1 virus infection is linked to SOCS3 induction. Under their conditions, swine epithelial cells were expressing SOCS3 proteins in response to both human H1N1 virus and HPAI H5N1 virus . In another study, it was demonstrated that influenza A virus induced NF-κB dependent SOCS3 gene expression in human epithelial cells, suppressing type I IFN response through the JAK/STAT pathway . In our conditions, SOCS3 was never significantly induced in response to H3N2 SIV. It remains to be shown whether these differences are due to different experimental parameters (e.g. virus strain and timing of measurements) or potential discrepancies between RNA levels and protein expression.
In NPTr cells, the Akt signaling pathway was not significantly induced, while we did observe the triggering of MAPK (ERK1/2 and p38) and JAK/STAT pathways. The absence of Akt activation was unexpected, since this signaling pathway was shown to be activated by influenza A virus, probably favoring viral replication [32, 52]. However, we assessed different pathway activations early in the infection process, and the protocol used in our study did not allow for the observation of Akt activation. The activation of JAK/STAT pathway has been reported in other species . This pathway controls type I IFN response. JAK/STAT signaling induced formation of a trimeric transcription factor ISGF3, consisting of STAT1, STAT2, and interferon regulatory factor 9 (IRF-9), which regulates the expression of several ISGs . Four MAPK cascades have been described including two isoforms of the ERK1/2, the Jun-terminal kinase (JNK), p38, and ERK5 [52, 53]. The signaling pathways convert extra- and intracellular signals into cellular responses, regulating cell activation, differentiation, immune responses, and proliferation [52, 53]. In our study we assessed the activation of ERK1/2 and p38, and activation of both cascades was observed. p38 predominantly responds to stress conditions and pro-inflammatory cytokines, whereas ERK1/2 senses mitogenic stimuli [52, 53]. MAPK signaling pathways and their role in the regulation of innate response and viral replication have not been well studied in porcine cells infected by influenza virus. However, a recent study showed a robust activation of ERK1/2 in influenza virus-infected swine macrophages  and demonstrated that the induction of RIG-I and MDA-5 depends on the activation of ERK1/2 and JNK in these cells. Interestingly and similar to our findings, these authors did not observe any IL-1β induction in response to the tested influenza virus . Unlike IL-1β expression, TNFα was clearly induced early at 3 hpi.
We also investigated to what extent inhibitors of the different pathways could impact SOCS1 transcription. Only JAK/STAT pathway inhibition was accompanied by a significant decrease in the expression of SOCS1 transcripts. This observation confirms the strong association of SOCS (particularly CISH, SOCS1, and SOCS3) with the JAK/STAT pathway. Nevertheless, in our conditions the association was only clear for SOCS1 and not for the two other SOCS family members. To further assess the impact of JAK/STAT pathway inhibition, we then looked at the expression of transcripts accounting for viral replication, viral recognition and interferon response. The treatment with JAK Inhibitor I was accompanied by a significant decrease in viral replication (as assessed by quantification of M-vRNA) and by a diminished expression of RIG-I, IFNβ, IFNλ1, PKR, and Mx1 transcripts 24 h post-infection. The decrease in viral replication remains unexplained, although it could be linked to the presence of more IFN in the culture supernatant consecutive to an inhibition of SOCS1. However, the lower level of expression of IFN transcripts argues against this hypothesis. In another set of experiments using another SIV strain (unpublished results) we used pJAK2(1001-1013) inhibitor [55, 56] and obtained similar results in both epithelial cells and macrophages. The viral replication was significantly reduced, as were the transcriptions of IFN and ISGs. Regarding the inhibition of MAPK pathways, we observed a decrease in viral replication as well as a decrease in the interferon response. Similarly, a study by Pinto et al.  demonstrated in human epithelial cells that the inhibition of MEK in the Raf/MEK/ERK pathway resulted in reduced virus titers and lower expression of pro-inflammatory cytokines [57, 58]. In another study on HPAI H5N1 virus in mice  it was also shown that the inhibition of p38 greatly reduced virus-induced cytokine expression concomitant with reduced viral titers.
In conclusion, our study showed the advantage of using various cellular and tissue systems in the study of the innate immune response to influenza virus. The kinetic approach enabled a comparative assessment of the innate immune response in the main targets of the virus and brought insights to viral triggering of major signaling pathways in NPTR cells. Our data suggests a link between SOCS1 and the JAK/STAT pathway in porcine respiratory epithelial cells. The inhibition of the JAK/STAT pathway using JAK Inhibitor I clearly reduced the response of IFN types I and III and the induction of SOCS1 at the transcript level in NPTr cells. Similarly the inhibition of MAPK pathway reduced viral replication and interferon response. Further investigations are required to understand precisely the role of SOCS in pigs. A better understanding of the host innate immune response and its underlying regulation mechanisms will help identify strategies for effective control of SIV.
This work was supported by grants from the Institut National de la Recherche Agronomique (INRA), Conseil Régional du Centre (France) (Mario Delgado-Ortega), and the Natural Science and Engineering Research Council of Canada (NSERC, grant 435887-2013). Dr Georg Herrler was supported by a grant from the German FluResearchNet, a nationwide research network on zoonotic influenza sponsored by the Ministry of Education and Research. We are thankful to UEPAO and UE PRC staff (INRA) for their invaluable help with broncho-alveolar lavages. We are also very grateful to Dr Hugh Townsend for his valuable advices to carry out the statistical analysis and to Dr Colette Wheler for the careful revision of the manuscript. We thank Ralf Dürrwald (IDT Biologika GmbH, Dessau-Rosslau, Germany) for the gift of influenza virus strain A/Swine/ Bissendorf/IDT1864/2003 and Sabine Uhlenbruck for the training in the preparation of the lung slices. The manuscript was published with permission of the Director of VIDO as manuscript # 689.
- Fablet C, Marois C, Kuntz-Simon G, Rose N, Dorenlor V, Eono F, Eveno E, Jolly JP, Le Devendec L, Tocqueville V, Quéguiner S, Gorin S, Kobisch M, Madec F: Longitudinal study of respiratory infection patterns of breeding sows in five farrow-to-finish herds. Vet Microbiol. 2011, 147: 329-339. 10.1016/j.vetmic.2010.07.005.View ArticlePubMedGoogle Scholar
- Fablet C, Marois-Crehan C, Simon G, Grasland B, Jestin A, Kobisch M, Madec F, Rose N: Infectious agents associated with respiratory diseases in 125 farrow-to-finish pig herds: a cross-sectional study. Vet Microbiol. 2012, 157: 152-163. 10.1016/j.vetmic.2011.12.015.View ArticlePubMedGoogle Scholar
- Opriessnig T, Gimenez-Lirola LG, Halbur PG: Polymicrobial respiratory disease in pigs. Anim Health Res Rev. 2011, 12: 133-148. 10.1017/S1466252311000120.View ArticlePubMedGoogle Scholar
- Choi YK, Goyal SM, Joo HS: Retrospective analysis of etiologic agents associated with respiratory diseases in pigs. Can Vet J. 2003, 44: 735-737.PubMed CentralPubMedGoogle Scholar
- Fablet C, Marois C, Dorenlor V, Eono F, Eveno E, Jolly JP, Le Devendec L, Kobisch M, Madec F, Rose N: Bacterial pathogens associated with lung lesions in slaughter pigs from 125 herds. Res Vet Sci. 2012, 93: 627-630. 10.1016/j.rvsc.2011.11.002.View ArticlePubMedGoogle Scholar
- Crisci E, Mussa T, Fraile L, Montoya M: Review: influenza virus in pigs. Mol Immunol. 2013, 55: 200-211. 10.1016/j.molimm.2013.02.008.View ArticlePubMedGoogle Scholar
- Brown IH: History and epidemiology of Swine influenza in Europe. Curr Top Microbiol Immunol. 2013, 370: 133-146.PubMedGoogle Scholar
- Kuntz-Simon G, Madec F: Genetic and antigenic evolution of swine influenza viruses in Europe and evaluation of their zoonotic potential. Zoonoses Public Health. 2009, 56: 310-325. 10.1111/j.1863-2378.2009.01236.x.View ArticlePubMedGoogle Scholar
- Kyriakis CS, Brown IH, Foni E, Kuntz-Simon G, Maldonado J, Madec F, Essen SC, Chiapponi C, Van Reeth K: Virological surveillance and preliminary antigenic characterization of influenza viruses in pigs in five European countries from 2006 to 2008. Zoonoses Public Health. 2011, 58: 93-101. 10.1111/j.1863-2378.2009.01301.x.View ArticlePubMedGoogle Scholar
- Kyriakis CS, Rose N, Foni E, Maldonado J, Loeffen WL, Madec F, Simon G, Van Reeth K: Influenza A virus infection dynamics in swine farms in Belgium, France, Italy and Spain, 2006-2008. Vet Microbiol. 2013, 162: 543-550. 10.1016/j.vetmic.2012.11.014.View ArticlePubMedGoogle Scholar
- Evans SE, Xu Y, Tuvim MJ, Dickey BF: Inducible innate resistance of lung epithelium to infection. Annu Rev Physiol. 2010, 72: 413-435. 10.1146/annurev-physiol-021909-135909.PubMed CentralView ArticlePubMedGoogle Scholar
- Vareille M, Kieninger E, Edwards MR, Regamey N: The airway epithelium: soldier in the fight against respiratory viruses. Clin Microbiol Rev. 2011, 24: 210-229. 10.1128/CMR.00014-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Ehrhardt C, Seyer R, Hrincius ER, Eierhoff T, Wolff T, Ludwig S: Interplay between influenza A virus and the innate immune signaling. Microbes Infect. 2010, 12: 81-87. 10.1016/j.micinf.2009.09.007.View ArticlePubMedGoogle Scholar
- Mussa T, Ballester M, Silva-Campa E, Baratelli M, Busquets N, Lecours MP, Dominguez J, Amadori M, Fraile L, Hernandez J, Montoya M: Swine, human or avian influenza viruses differentially activates porcine dendritic cells cytokine profile. Vet Immunol Immunopathol. 2013, 154: 25-35. 10.1016/j.vetimm.2013.04.004.View ArticlePubMedGoogle Scholar
- Taubenberger JK, Morens DM: The pathology of influenza virus infections. Annu Rev Pathol. 2008, 3: 499-522. 10.1146/annurev.pathmechdis.3.121806.154316.PubMed CentralView ArticlePubMedGoogle Scholar
- Kreijtz JH, Fouchier RA, Rimmelzwaan GF: Immune responses to influenza virus infection. Virus Res. 2011, 162: 19-30. 10.1016/j.virusres.2011.09.022.View ArticlePubMedGoogle Scholar
- Pang IK, Iwasaki A: Inflammasomes as mediators of immunity against influenza virus. Trends Immunol. 2011, 32: 34-41. 10.1016/j.it.2010.11.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Ferrari M, Scalvini A, Losio MN, Corradi A, Soncini M, Bignotti E, Milanesi E, Ajmone-Marsan P, Barlati S, Bellotti D, Tonelli M: Establishment and characterization of two new pig cell lines for use in virological diagnostic laboratories. J Virol Methods. 2003, 107: 205-212. 10.1016/S0166-0934(02)00236-7.View ArticlePubMedGoogle Scholar
- Massin P, Kuntz-Simon G, Barbezange C, Deblanc C, Oger A, Marquet-Blouin E, Bougeard S, van der Werf S, Jestin V: Temperature sensitivity on growth and/or replication of H1N1, H1N2 and H3N2 influenza A viruses isolated from pigs and birds in mammalian cells. Vet Microbiol. 2010, 142: 232-241. 10.1016/j.vetmic.2009.10.012.View ArticlePubMedGoogle Scholar
- Abd El Rahman S, Winter C, El-Kenawy A, Neumann U, Herrler G: Differential sensitivity of well-differentiated avian respiratory epithelial cells to infection by different strains of infectious bronchitis virus. J Virol. 2010, 84: 8949-8952. 10.1128/JVI.00463-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Goris K, Uhlenbruck S, Schwegmann-Wessels C, Kohl W, Niedorf F, Stern M, Hewicker-Trautwein M, Bals R, Taylor G, Braun A, Bicker G, Kietzmann M, Herrler G: Differential sensitivity of differentiated epithelial cells to respiratory viruses reveals different viral strategies of host infection. J Virol. 2009, 83: 1962-1968. 10.1128/JVI.01271-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Meng F, Punyadarsaniya D, Uhlenbruck S, Hennig-Pauka I, Schwegmann-Wessels C, Ren X, Durrwald R, Herrler G: Replication characteristics of swine influenza viruses in precision-cut lung slices reflect the virulence properties of the viruses. Vet Res. 2013, 44: 110-10.1186/1297-9716-44-110.PubMed CentralView ArticlePubMedGoogle Scholar
- Punyadarsaniya D, Liang CH, Winter C, Petersen H, Rautenschlein S, Hennig-Pauka I, Schwegmann-Wessels C, Wu CY, Wong CH, Herrler G: Infection of differentiated porcine airway epithelial cells by influenza virus: differential susceptibility to infection by porcine and avian viruses. PLoS One. 2011, 6: e28429-10.1371/journal.pone.0028429.PubMed CentralView ArticlePubMedGoogle Scholar
- Londt BZ, Brookes SM, Nash BJ, Nunez A, Kelly MD, Garcon F, Graham SP, Brown IH: Enhanced infectivity of H5N1 highly pathogenic avian influenza (HPAI) virus in pig ex vivo respiratory tract organ cultures following adaptation by in vitro passage. Virus Res. 2013, 178: 383-391. 10.1016/j.virusres.2013.09.015.View ArticlePubMedGoogle Scholar
- Londt BZ, Brookes SM, Nash BJ, Nunez A, Stagg DA, Brown IH: The infectivity of pandemic 2009 H1N1 and avian influenza viruses for pigs: an assessment by ex vivo respiratory tract organ culture. Influenza Other Respir Viruses. 2013, 7: 393-402. 10.1111/j.1750-2659.2012.00397.x.View ArticlePubMedGoogle Scholar
- Nunes SF, Murcia PR, Tiley LS, Brown IH, Tucker AW, Maskell DJ, Wood JL: An ex vivo swine tracheal organ culture for the study of influenza infection. Influenza Other Respir Viruses. 2010, 4: 7-15. 10.1111/j.1750-2659.2009.00119.x.View ArticlePubMedGoogle Scholar
- Van Poucke SG, Nicholls JM, Nauwynck HJ, Van Reeth K: Replication of avian, human and swine influenza viruses in porcine respiratory explants and association with sialic acid distribution. Virol J. 2010, 7: 38-10.1186/1743-422X-7-38.PubMed CentralView ArticlePubMedGoogle Scholar
- European Convention for the Protection of Vertebrate Animals used for Experimental and Other Scientific Purposes. [http://conventions.coe.int/Treaty/en/Treaties/Html/123.htm]
- Protocol of amendment to the European Convention for the Protection of Vertebrate Animals used for Experimental and other Scientific Purposes. [http://conventions.coe.int/Treaty/en/Treaties/Html/170.htm]
- White MR, Tecle T, Crouch EC, Hartshorn KL: Impact of neutrophils on antiviral activity of human bronchoalveolar lavage fluid. Am J Physiol Lung Cell Mol Physiol. 2007, 293: L1293-L1299. 10.1152/ajplung.00266.2007.View ArticlePubMedGoogle Scholar
- Vietmeier J, Niedorf F, Baumer W, Martin C, Deegen E, Ohnesorge B, Kietzmann M: Reactivity of equine airways–a study on precision-cut lung slices. Vet Res Commun. 2007, 31: 611-619. 10.1007/s11259-007-3501-y.View ArticlePubMedGoogle Scholar
- Shin YK, Liu Q, Tikoo SK, Babiuk LA, Zhou Y: Influenza A virus NS1 protein activates the phosphatidylinositol 3-kinase (PI3K)/Akt pathway by direct interaction with the p85 subunit of PI3K. J Gen Virol. 2007, 88: 13-18. 10.1099/vir.0.82419-0.View ArticlePubMedGoogle Scholar
- Delgado-Ortega M, Melo S, Meurens F: Expression of SOCS1-7 and CIS mRNA in porcine tissues. Vet Immunol Immunopathol. 2011, 144: 493-498. 10.1016/j.vetimm.2011.08.002.View ArticlePubMedGoogle Scholar
- Meurens F, Berri M, Siggers RH, Willing BP, Salmon H, Van Kessel AG, Gerdts V: Commensal bacteria and expression of two major intestinal chemokines, TECK/CCL25 and MEC/CCL28, and their receptors. PLoS One. 2007, 2: e677-10.1371/journal.pone.0000677.PubMed CentralView ArticlePubMedGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-PubMed CentralView ArticlePubMedGoogle Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8: R19-10.1186/gb-2007-8-2-r19.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao S, Fernald RD: Comprehensive algorithm for quantitative real-time polymerase chain reaction. J Comput Biol. 2005, 12: 1047-1064. 10.1089/cmb.2005.12.1047.PubMed CentralView ArticlePubMedGoogle Scholar
- Nygard AB, Jorgensen CB, Cirera S, Fredholm M: Selection of reference genes for gene expression studies in pig tissues using SYBR green qPCR. BMC Mol Biol. 2007, 8: 67-10.1186/1471-2199-8-67.PubMed CentralView ArticlePubMedGoogle Scholar
- Sang Y, Rowland RR, Hesse RA, Blecha F: Differential expression and activity of the porcine type I interferon family. Physiol Genomics. 2010, 42: 248-258. 10.1152/physiolgenomics.00198.2009.View ArticlePubMedGoogle Scholar
- Richt JA, Lager KM, Clouser DF, Spackman E, Suarez DL, Yoon KJ: Real-time reverse transcription-polymerase chain reaction assays for the detection and differentiation of North American swine influenza viruses. J Vet Diagn Invest. 2004, 16: 367-373. 10.1177/104063870401600501.View ArticlePubMedGoogle Scholar
- Meurens F, Berri M, Auray G, Melo S, Levast B, Virlogeux-Payant I, Chevaleyre C, Gerdts V, Salmon H: Early immune response following Salmonella enterica subspecies enterica serovar Typhimurium infection in porcine jejunal gut loops. Vet Res. 2009, 40: 5-10.1051/vetres:2008043.View ArticlePubMedGoogle Scholar
- Erkens T, Van Poucke M, Vandesompele J, Goossens K, Van Zeveren A, Peelman LJ: Development of a new set of reference genes for normalization of real-time RT-PCR data of porcine backfat and longissimus dorsi muscle, and evaluation with PPARGC1A. BMC Biotechnol. 2006, 6: 41-10.1186/1472-6750-6-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Knosp CA, Johnston JA: Regulation of CD4+ T-cell polarization by suppressor of cytokine signalling proteins. Immunology. 2012, 135: 101-111. 10.1111/j.1365-2567.2011.03520.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Ioannidis I, Ye F, McNally B, Willette M, Flano E: Toll-like receptor expression and induction of type I and type III interferons in primary airway epithelial cells. J Virol. 2013, 87: 3261-3270. 10.1128/JVI.01956-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Genini S, Delputte PL, Malinverni R, Cecere M, Stella A, Nauwynck HJ, Giuffra E: Genome-wide transcriptional response of primary alveolar macrophages following infection with porcine reproductive and respiratory syndrome virus. J Gen Virol. 2008, 89: 2550-2564. 10.1099/vir.0.2008/003244-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Razzuoli E, Villa R, Sossi E, Amadori M: Reverse transcription real-time PCR for detection of porcine interferon alpha and beta genes. Scand J Immunol. 2011, 74: 412-418. 10.1111/j.1365-3083.2011.02586.x.View ArticlePubMedGoogle Scholar
- Delgado-Ortega M, Marc D, Dupont J, Trapp S, Berri M, Meurens F: SOCS proteins in infectious diseases of mammals. Vet Immunol Immunopathol. 2013, 151: 1-19. 10.1016/j.vetimm.2012.11.008.View ArticlePubMedGoogle Scholar
- Yoshimura A, Naka T, Kubo M: SOCS proteins, cytokine signalling and immune regulation. Nat Rev Immunol. 2007, 7: 454-465. 10.1038/nri2093.View ArticlePubMedGoogle Scholar
- Pothlichet J, Chignard M, Si-Tahar M: Cutting edge: innate immune response triggered by influenza A virus is negatively regulated by SOCS1 and SOCS3 through a RIG-I/IFNAR1-dependent pathway. J Immunol. 2008, 180: 2034-2038.View ArticlePubMedGoogle Scholar
- Nelli RK, Dunham SP, Kuchipudi SV, White GA, Baquero-Perez B, Chang P, Ghaemmaghami A, Brookes SM, Brown IH, Chang KC: Mammalian innate resistance to highly pathogenic avian influenza H5N1 virus infection is mediated through reduced proinflammation and infectious virus release. J Virol. 2012, 86: 9201-9210. 10.1128/JVI.00244-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Pauli EK, Schmolke M, Wolff T, Viemann D, Roth J, Bode JG, Ludwig S: Influenza A virus inhibits type I IFN signaling via NF-kappaB-dependent induction of SOCS-3 expression. PLoS Pathog. 2008, 4: e1000196-10.1371/journal.ppat.1000196.PubMed CentralView ArticlePubMedGoogle Scholar
- Ehrhardt C, Marjuki H, Wolff T, Nurnberg B, Planz O, Pleschka S, Ludwig S: Bivalent role of the phosphatidylinositol-3-kinase (PI3K) during influenza virus infection and host cell defence. Cell Microbiol. 2006, 8: 1336-1348. 10.1111/j.1462-5822.2006.00713.x.View ArticlePubMedGoogle Scholar
- Ludwig S, Pleschka S, Planz O, Wolff T: Ringing the alarm bells: signalling and apoptosis in influenza virus infected cells. Cell Microbiol. 2006, 8: 375-386. 10.1111/j.1462-5822.2005.00678.x.View ArticlePubMedGoogle Scholar
- Gao W, Sun W, Qu B, Cardona CJ, Powell K, Wegner M, Shi Y, Xing Z: Distinct regulation of host responses by ERK and JNK MAP kinases in swine macrophages infected with pandemic (H1N1) 2009 influenza virus. PLoS One. 2012, 7: e30328-10.1371/journal.pone.0030328.PubMed CentralView ArticlePubMedGoogle Scholar
- Flowers LO, Johnson HM, Mujtaba MG, Ellis MR, Haider SM, Subramaniam PS: Characterization of a peptide inhibitor of Janus kinase 2 that mimics suppressor of cytokine signaling 1 function. J Immunol. 2004, 172: 7510-7518.View ArticlePubMedGoogle Scholar
- Waiboci LW, Ahmed CM, Mujtaba MG, Flowers LO, Martin JP, Haider MI, Johnson HM: Both the suppressor of cytokine signaling 1 (SOCS-1) kinase inhibitory region and SOCS-1 mimetic bind to JAK2 autophosphorylation site: implications for the development of a SOCS-1 antagonist. J Immunol. 2007, 178: 5058-5068.View ArticlePubMedGoogle Scholar
- Pinto R, Herold S, Cakarova L, Hoegner K, Lohmeyer J, Planz O, Pleschka S: Inhibition of influenza virus-induced NF-kappaB and Raf/MEK/ERK activation can reduce both virus titers and cytokine expression simultaneously in vitro and in vivo. Antiviral Res. 2011, 92: 45-56. 10.1016/j.antiviral.2011.05.009.View ArticlePubMedGoogle Scholar
- Planz O: Development of cellular signaling pathway inhibitors as new antivirals against influenza. Antiviral Res. 2013, 98: 457-468. 10.1016/j.antiviral.2013.04.008.View ArticlePubMedGoogle Scholar
- Borgeling Y, Schmolke M, Viemann D, Nordhoff C, Roth J, Ludwig S: Inhibition of p38 mitogen-activated protein kinase impairs influenza virus-induced primary and secondary host gene responses and protects mice from lethal H5N1 infection. J Biol Chem. 2013, 289: 13-27.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.