669
Views
2
CrossRef citations to date
0
Altmetric
Research Article

Serogrouping and seroepidemiology of North European hantaviruses using a novel broadly targeted synthetic nucleoprotein antigen array

ORCID Icon, , ORCID Icon, , ORCID Icon, ORCID Icon & ORCID Icon show all
Article: 1350086 | Received 04 Jan 2017, Accepted 24 Jun 2017, Published online: 26 Jul 2017

ABSTRACT

Introduction: Hantaviruses are globally distributed zoonotic pathogens. Great diversity and high antigenic cross-reactivity makes diagnosis by traditional methods cumbersome.

Materials and methods: ‘Megapeptides’, 119–120-mers from the amino terminus of the nucleoprotein of 16 hantaviruses, representing the four major branches of the hantavirus phylogenetic tree, were utilized in a novel IgG-based hantavirus suspension multiplex immunoassay (HSMIA) for detection of past hantavirus infections in 155 North European human samples. We compared HSMIA with established EIAs and focus reduction neutralization test (FRNT).

Results and discussion: The Puumala hantavirus (PUUV) component in the HSMIA gave concordant results with a PUUV IgG EIA in 142 sera from Northern Sweden (of which 31 were EIA positive, 7 borderline and 104 EIA negative, sensitivity 30/31 = 97%, specificity 104/ 104 = 100%, 134/135 = 99% concordance), with another immunoassay in 40 PUUV IgG positive sera from Finland (36/40 = 90% sensitivity), and was concordant in 8 of 11 cases with PUUV and DOBV neutralization titers, respectively. Two major IgG reactivity patterns were found: (i) a PUUV-specific pattern covering phylogroup IV and its serogroups B and C; and (ii) a Dobrava virus (DOBV)-specific pattern, covering the serogroup A portion of phylogroup III. In addition, we found several minor patterns with reactivity to only one or two megapeptides indicating additional hantaviruses infecting humans in the Swedish and Finnish populations.

Conclusion: The broadly reactive and rational HSMIA yielded results highly correlated with the established PUUV EIAs and the NT results. It is a sensitive and specific assay, which will be suited for efficient serosurveillance of hantaviruses in humans. Its use in animals should be further investigated.

Introduction

Hantaviruses are globally distributed and can cause serious human infections – hemorrhagic fever with renal syndrome (HFRS) or hantavirus cardiopulmonary syndrome (HCPS) [Citation1Citation5]. At present, 24 unique hantaviruses have been described and around 15 additional viruses await taxonomic acceptance [Citation1]. New hantaviruses are currently discovered at a high rate in different mammals. Only a minority of the hantaviruses have been isolated from humans [Citation2]. Hantaviral geographic distribution and host range are rather strict [Citation3Citation14]. Host mobility, including human-mediated transportation, may be the reason why Seoul virus (SEOV), originally found in South East Asia, now has been found in rats and humans in Belgium, UK, France, Sweden, and the Netherlands [Citation11]. These emerging zoonotic infections [Citation15Citation18] need to be monitored. Current diagnostic methods [Citation19] are mainly directed to a few viruses at a time, and thereby do not cover the great diversity of the hantaviruses.

We previously used synthetic, up to 30-mer, peptides in hantavirus serology [Citation20,Citation21] and multiplex serology with other microbes [Citation22,Citation23], later extended to 50- to 220-mers [Citation24]. Long synthetic peptides (‘megapeptides’) have occasional defects in their sequences due to incomplete amino acid coupling, but most molecules will have a conformation that is similar enough to that of the complete native molecule, thus mimicking conformational epitopes. Posttranslational modifications are not included, but simplicity of creation and absence of unwanted cellular antigens partially outweigh this disadvantage. We here used megapeptides longer than 100 amino acids for hantavirus serology, using our Suspension Multiplex Immuno Assay (SMIA) technique [Citation22Citation24]. It has earlier been shown that the first 107 amino acids of the amino terminus of hantavirus N proteins contain broadly cross-reactive immunodominant epitopes [Citation20,Citation25Citation27], dependent on two antiparallel beta helices [Citation25Citation27]. The SMIA array contained synthetic 119–120-amino acid long peptides from the amino terminus of the nucleoprotein from all known major branches of the hantavirus phylogenetic tree [Citation2], from phylogroup I/serogroup D; Thottapalayam virus (TPMV), from phylogroup II/serogroup D; Nova virus (NVAV), Laibin virus (LBV), from phylogroup III/Serogroups A and E; Dobrava virus (DOBV), Hantaan virus (HTNV), Andes virus (ANDV), Seoul virus (SEOV), Thailand virus (THAIV), Asama virus (ASAV), Oxbow virus (OXBV), Seewis virus (SWSV) and from phylogroup IV/Serogroups B and C; Puumala virus (PUUV), Tula virus (TULV), Sin Nombre virus (SNV), Laguna Negra virus (LANV), and Rockport virus (RKPV). The intention was to cover as many as possible of the known hantaviruses. We tested this array with sets of known hantavirus positive sera, exploring its capability to detect and type by using multivariate analysis. Using known PUUV and DOBV antibody positive sera we found that this indeed was possible. They displayed two major patterns of reactivity, conforming to either DOBV or PUUV specificity. We also got indications of infections with yet-to-be discovered hantaviruses in blood donors and other persons from the Nordic countries.

Materials and methods

Synthetic peptides

The peptide sequences are shown in . All had the spacer NH2-PEG6-His6-PEG6 at their amino terminus, where PEG6 is hexaethylene glycol and His6 is hexahistidine (enabling coupling control via anti-His antibodies). The peptides were synthesized under conditions optimized for maximum coupling yield of long peptides, see e.g. [Citation24]. Lyophilized peptides were dissolved in dimethyl sulfoxide (DMSO) (Sigma D2650).

Figure 1. Megapeptides from the N terminus of hantavirus N proteins used in the present work. Peptide sequences were aligned by Muscle. A majority consensus is shown above the alignment. Positions with identity to the consensus are shown as dots. Conservation is also shown with symbols below the alignment, i.e. * = Fully conserved, : = semiconserved and . = conservation with exchange of similar amino acids. GenBank Ids are shown for each sequence. The tree (phylogram) of was derived from this figure.

Figure 1. Megapeptides from the N terminus of hantavirus N proteins used in the present work. Peptide sequences were aligned by Muscle. A majority consensus is shown above the alignment. Positions with identity to the consensus are shown as dots. Conservation is also shown with symbols below the alignment, i.e. * = Fully conserved, : = semiconserved and . = conservation with exchange of similar amino acids. GenBank Ids are shown for each sequence. The tree (phylogram) of Figure 2 was derived from this figure.

Figure 2. Part of the total result table. Each row represents one serum. The first 11 rows are from the FRNT collection, the subsequent 28 are a portion from the North Swedish collection (totally 104). Megapeptides (heading) and their results (MFI) are arranged in phylogenetic order, according to a dendrogram based on megapeptide protein similarity obtained after Muscle alignment at the European Bioinformatics Centre (EBC) home page. The megapeptide components which contributed most to principal components 1 and 2, ‘Pc1’ and ‘Pc2’, respectively are also indicated. Results were automatically classified (shown in column ‘seroclass’) as described in the text. All results are averages of two determinations and subtracted with the MFI of the ‘naked’ bead. MFI values higher than 400 MFI and 1000 MFI have been highlighted, in a graded fashion. Nova peptide reactions were relatively strong in the FRNT sera. We assume this was due to heat inactivation.

Figure 2. Part of the total result table. Each row represents one serum. The first 11 rows are from the FRNT collection, the subsequent 28 are a portion from the North Swedish collection (totally 104). Megapeptides (heading) and their results (MFI) are arranged in phylogenetic order, according to a dendrogram based on megapeptide protein similarity obtained after Muscle alignment at the European Bioinformatics Centre (EBC) home page. The megapeptide components which contributed most to principal components 1 and 2, ‘Pc1’ and ‘Pc2’, respectively are also indicated. Results were automatically classified (shown in column ‘seroclass’) as described in the text. All results are averages of two determinations and subtracted with the MFI of the ‘naked’ bead. MFI values higher than 400 MFI and 1000 MFI have been highlighted, in a graded fashion. Nova peptide reactions were relatively strong in the FRNT sera. We assume this was due to heat inactivation.

The epitope prediction tools at the Immune Epitope Database [Citation28] and 3D structures found in the Protein structure DataBase (see e.g. PDB id 2K48) [Citation29] were utilized for optimization of megapeptide design.

Coupling of peptides to magnetic Luminex beads

The peptides were coupled to Luminex carboxylated magnetic beads (MagPlex-C microspheres, Luminex Corp, TX, USA) according to the manufacturer’s protocol ‘Sample protocol for two step carbodiimide coupling of protein to magnetic beads’ (http://cdn2.hubspot.net/hub/128032/file-213083097-pdf/Luminex-xMap_Cookbook.pdf).

The coupling was made with 200 µl of the stock microspheres, containing 1.25 × 107 beads per ml, or a smaller proportion of this volume. In the latter case, all volumes in the recipe were decreased proportionally. After the coupling, beads were incubated with 0.5 ml of PBS containing 0.05% (v/v) Tween 20 and 50 mM Tris in the dark for 15 min on a rocking mixer at room temperature to block unreacted carboxyl groups with primary amines. The microspheres were then washed once with 0.5 ml StabilGuard (SurModics, Eden Prairie, MN, USA, #SG01-1000) using the magnetic separator. The bead pellet was finally resuspended in 400 µl StabilGuard. This created a bead mixture consisting of 6250 beads/µl. The coupled beads were stored at 4°C in the dark.

Serum samples

Sera typed by focus reduction neutralization test (FRNT)

These samples (n = 11) were described in a previous study on hantavirus infections in the Netherlands [Citation30]. These sera were initially identified by screening by IFT and EIAs and subsequently confirmed (serotyped) by FRNT. This study was exempted from ethical review of human subject research by the Medical Ethical Review Committee of the Erasmus MC Medical Centre, University of Rotterdam. All data have been anonymized and are not attributable to individual patients.

Sera from Northern Sweden

Sera from Northern Sweden (n = 142) were obtained as part of a serological survey for PUUV-reactive antibodies [Citation31]. All samples were from healthy adults, and were obtained with ethical permission number Dnr2011/819-31/3). All samples were initially tested by a Puumala IgG EIA, see below.

Sera from Finland

Sera from Finland (n = 40) were sent to Helsinki University Central Hospital Laboratory Service (HUSLAB, Department of Virology and Immunology, Zoonosis Unit, Helsinki, Finland) for PUUV IgM and IgG analysis for suspicion of PUUV infection and had positive IgG but negative IgM result. The collection for the study was performed under HUSLAB’s research permission for OV (granted on 14 June 2013). The PUUV IgM and IgG reference results are accredited (SFS-EN ISO/IEC 17025 and SFS-EN ISO 15189; Finnish Standards Association). The IgM test (data not utilized in the present paper) was an IgM capture enzyme immunoassay based on recombinant PUUV-N produced in the baculovirus system [Citation32]. The IgG test was an in-house immunofluorescence assay (IFA) based on PUUV-infected acetone-fixed Vero E6 cells [Citation33,Citation34].

Blood donor samples

Blood donor samples (n = 89) from Uppsala, obtained anonymously with recorded donor consent according to the Swedish Biobank law and ethical consent UPS_01_367, from the Uppsala University Hospital blood bank. The Swedish Biobank law states that sera obtained with donor consent may be reused for developmental purposes if the purpose is related to the original purpose of taking the sample.

Serologies

Hantavirus suspension multiplex immunoassay (HSMIA)

Fifty µl of serum diluted 1/20 in PBS, pH 7.4, containing 0.05% (v/v) Tween 20, 50 mM Tris and 2% (v/v) Prionex (Sigma-Aldrich #81662) (PBSTP) was added to wells of a round bottom 96-well microtiter plate (Greiner #104650). Fifty µl of a vortexed and sonicated bead mixture consisting of 25 beads/µl was then added to each well. The plate was then incubated in the dark with gentle rotation for 1 h at 37°C. Wells were washed once with PBS using a magnetic plate separator (Life technologies #A14179). Beads were resuspended in 50 µl of PBSTP and 50 µl of biotinylated protein G (Pierce, article nr 29988) (4 µg/ml PBSTP) was added to each well. After 30 min of incubation in the dark with rotation at 37°C, the wells were washed once in PBS. Beads were resuspended again in 50 µl of PBSTP. Fifty µl of streptavidin-phycoerythrin (SA-PE) (InVitrogen-Thermo-Fisher, article nr S-866) (4 µg/ml in PBSTP) was added per well. The plate was incubated in the dark with gentle rotation for 15 min at 37°C. Beads were washed once in PBS before they were resuspended in 100 µl of PBS and analyzed in a Luminex-200 (Luminex corporation, Austin, Texas, USA) instrument following the manufacturer’s instructions. A minimum of 100 events per bead set were read and the median value obtained for each reaction event per bead set. A sample volume of 75 µl was analyzed for each sample. All samples were analyzed in duplicate and average readings calculated. To detect any antibody binding to the beads themselves, a naked non-peptide-containing (‘blank’) bead was included. A control His6 bead was also included in the SMIA reaction. The His6 tag allowed monitoring of coupling of the megapeptides to the magnetic beads using anti-hexahistidine antibodies (Antibodies on line, ABIN100493). We did not observe false positive reactions due to anti-His6 in the sera tested for this paper. One negative serum control, where PBSTP instead of serum was added, was also used in all experiments. A cutoff was first calculated as the PUUV average MFI of negative samples (OD < 0.1 in Puumala EIA) + 3 standard deviations, resulting in 108 MFI. In earlier SMIA work, we used a more rigorous cutoff of 200 MFI [Citation22,Citation23]. This higher value was subsequently used. The reproducibility of HSMIA was determined in 40 PUUV positive sera from Northern Sweden. They gave a standard deviation of 11% after retesting.

Comparative serological tests

A PUUV IgG EIA, based on Mab-captured baculovirus-expressed PUUV nucleocapsid protein (N), was performed as previously described [Citation35]. The Finnish sera were tested with a µ capture assay [Citation35] for PUUV IgM, and were positive for PUUV IgG according to an immunofluorescence (IFA) against PUUV and a TR-FRET assay [Citation36].

Focus reduction neutralization test (FRNT)

To confirm hantavirus-reactivity and -specificity, a focus reduction neutralization test (FRNT), a highly virus species-specific technique for hantavirus serology, was performed as previously described [Citation35,Citation37,Citation38]. Briefly, serially diluted serum samples were mixed with virus and incubated on confluent Vero E6 cell monolayers in six-well tissue culture plates for 7–12 days depending on the hantavirus. Monkey or rodent anti-hantavirus polyclonal sera were used as primary antibodies (35b, 35c) and peroxidase-labelled goat anti-human or anti-rabbit IgG (BioRad Laboratories, Hercules, CA, USA; and Jackson Laboratories, Bar Harbor, ME, USA) as conjugates, followed by 3, 3ʹ, 5, 5ʹ-tetramethylbenzidine substrate (Sigma, Stockholm, Sweden). An 80% reduction of the number of foci was used as cut-off for virus neutralization.

Statistical evaluation, computerized handling of results and bioinformatics

Principal Component Analysis (PCA) was performed using the multivariate program Unscrambler (CAMO AS, Bergen, Norway). Values from 282 samples, each comprising 16 observations, were subjected to the exact singular value decomposition PCA method. No missing data occurred in the dataset. All results were subtracted, averaged and collected in a .dbf table with columns in phylogenetic order, and categorized, using a program (‘SeroClass’) written by JB in Visual Foxpro. Categorization followed the two major principal components and the three minor patterns not detected by PCA. The PCA defined PUUV and DOBV defined patterns included weights (PCA loadings), for each megapeptide pattern component. The sum of MFI times weight for each megapeptide had to reach a threshold of 300. The inclusion criterion for the minor patterns was that the first MFI of the pattern-defining megapeptide was stronger than MFIs of other megapeptides for that serum, and higher than 300.

Results

Design of megapeptides and creation of a broad hantavirus serological array

The N terminus of the N protein was chosen as the antigenic target because of its earlier proven high antigenicity, high conservation and high cross-reactivity. Before megapeptide synthesis, the surface availability of epitopes of crystallized hantaviral N proteins was evaluated using the Discotope program at IEDB [Citation28,Citation39]. The N terminus up to about amino acid 120 forms a prominent, highly exposed, protrusion (not shown), predicted to contain conformational and linear epitopes, corroborated by experimentation [Citation32,Citation40,Citation41].

The absence of glycosylation in the N protein was also judged to be favorable for using synthetic antigen mimicks. We therefore synthesized NH2-PEG6-His6-PEG6 119–120mers from the N terminus of 16 hantaviral N proteins (), intending to cover all known hantavirus phylogroups. At least one representative of the major phylo- and serogroups was included (). The multiplexing capability and the high precision of quantification of the SMIA results allowed us to create a rational hantavirus serological grouping and detection system, HSMIA.

The HSMIA system allowed us to process the 282 samples in a short time, generating 16*282 = 4512 IgG reactivity values. Initial titrations indicated that a serum dilution of 1/20 gave strong reactions of up to about 5000 MFI and background values of around 40 MFI, which after background subtraction (subtraction of the value of the naked bead) became close to 0, yielding both high sensitivity and specificity.

Analysis of groups of sera

We first analyzed 11 sera that had earlier been typed by neutralization. Three were DOBV-specific and eight PUUV-specific, giving NT titers of >40 up to 2560. The HSMIA results followed the neutralization results, with some exceptions (). The three sera, which had a high DOBV NT titer also showed a DOBV HSMIA pattern, while five out of eight sera that were PUUV FRNT-specific also showed a PUUV HSMIA pattern. The remaining three PUUV FRNT-specific sera only reacted weakly (below 200 MFI) with PUUV, TULV, and THAIV in HSMIA. The FRNT sera had been heat inactivated, and gave higher reactions with the naked beads, for unknown reasons also with the Nova megapeptide. Results by the universal positive control – a rabbit anti-hexahistidine serum (not shown), assured us that the hantaviral megapeptides were antigenically active.

We then analyzed 104 sera from Northern Sweden which had previously been found to be negative (OD < 0.1) using the PUUV virus IgG EIA [Citation42] (). The bottom 28 rows in are from the serosurvey in Northern Sweden. One of these sera had an ‘Asama’ (ASAV) reactivity by HSMIA (727 MFI; see below). The remaining 103 became negative (with an MFI pattern sum cutoff set to 300 MFI). For example, the average MFI of the DOBV megapeptide was 6.3 (SD 15.8) after background subtraction and 41.0 (SD 45.3) before subtraction, and of the PUUV megapeptide 9.0 (SD 29.0) after subtraction, and 44.0 (SD 44.9) before subtraction. The concordance for the IgG reactivity of the two megapeptides with PUUV IgG EIA for negativity was 103/104, i.e. 99%.

Table 1. Overview of sera and results.

Subsequently we analyzed 31 samples from Northern Sweden by an EIA OD of over 0.2 (). A total of 30 reacted with the PUUV megapeptide with an MFI of over 200 (after background subtraction), giving a concordance of 30/31 = 97%). The correlation between the two methods is shown in . All 30 were classified as PUUV-specific (see below). The 31st serum (EIA absorbance 0.653) reacted only weakly, below 200 MFI, with PUUV (83MFI) and SWSV (54 MFI). Calculating sensitivity and specificity from gave a HSMIA PUUV component sensitivity versus PUUV IgG EIA EIA of 30/31 = 97%, and a specificity of 104/104 = 100%. MFI of the PUUV component correlated clearly (linear regression coefficient (r) = 0.88) with the OD of the PUUV IgG EIA ().

Figure 3. Correlation between HSMIA PUUV MFI and SIIDC PUUV EIA OD, for the 134 samples from Northern Sweden. X axis: SIIDC PUUV EIA OD, Y axis: HSMIA PUUV megapeptide MFI after background subtraction. A regression line (yielding a linear regression coefficient of 0.88) is also shown.

Figure 3. Correlation between HSMIA PUUV MFI and SIIDC PUUV EIA OD, for the 134 samples from Northern Sweden. X axis: SIIDC PUUV EIA OD, Y axis: HSMIA PUUV megapeptide MFI after background subtraction. A regression line (yielding a linear regression coefficient of 0.88) is also shown.

When the Finnish ‘PUUV Immune’ (IgG positive) sera were analyzed by HSMIA, 36 of 40 became positive against the PUUV component. The 36 also had a ‘PUUV-specific’ pattern as defined below. Of the four remaining sera, two were strong ‘bead-binders’ (2428 and 681 MFI, respectively, with naked beads; see Materials) leading to a probable oversubtraction of the naked bead background value. Thus, 36 of the 38 (95%) samples which could be adequately analyzed were detected by HSMIA.

Three of 89 Uppsala blood donors were found PUUV component positive. All three had a ‘Puumala’ pattern (see below). However, two other sera were ‘Andes’ (ANDV) positive, and a further two sera were ‘Seewis’ (SWSV) positive. These reactivity patterns are defined below.

Definition of major and minor hantaserological patterns

PCA was performed on all 282 results from the panels. The positive serological results fell into two major categories, principal component 1 (Pc1), and principal component 2 (Pc2) (). The two principal components were orthogonal with respect to each other, i.e. they varied independently. Pc1 corresponded to PUUV-specific reactions, corroborated by the PUUV-specific neutralization of some of these sera and a high correlation with the PUUV EIA. Pc2 corresponded to DOBV-specific reactions, corroborated by the DOBV-specific neutralization of some of these sera. A total of 77 sera had a Pc1 (PUUV-specific) pattern of reactivity. Three sera had a Pc2 (DOBV-specific) pattern of reactivity (, ).

Figure 4. Principal component analysis. The two first principal components are PUUV (PC-1) and DOBV (PC-2) specific. ‘Loadings’ depict the influence of megapeptides for each principal component (Pc1 and Pc2).The two principal components together explained 95% of the variation in the dataset (not shown). The figure was produced by the multivariate program ‘Unscrambler’.

Figure 4. Principal component analysis. The two first principal components are PUUV (PC-1) and DOBV (PC-2) specific. ‘Loadings’ depict the influence of megapeptides for each principal component (Pc1 and Pc2).The two principal components together explained 95% of the variation in the dataset (not shown). The figure was produced by the multivariate program ‘Unscrambler’.

The distribution of Pc1 megapeptide serological reactions largely followed the distribution of a subset of phylogroup IV [Citation2], serogroup B [Citation26], while Pc2 megapeptide reactions were confined to phylogroup III [Citation2], serogroup A [Citation26] ( and ).

Figure 5. Distribution of hantavirus N protein megapeptides (bold face and boxed) in a hantavirus S segment phylogenetic tree. The tree was redrawn from that of Holmes and Zhang [Citation2]. Viruses known to give disease in humans are shown with a human figurine. Phylogroups I–IV and serogroups A–D [Citation26,Citation27] are indicated. Laibin virus (Phylogroup II) was recently discovered [Citation42]. The distinction between serogroups D1 and D2 was introduced in this paper. PCA ‘loadings’ indicate the degree of influence for each antigen in the respective principal component. They are here plotted according to phylogenetic relationship. Judging from the loadings of megapeptide antigens, Pc1 covers serogroups B and C while Pc2 covers serogroup A. However, Oxbow virus (Serogroup E) gave a minor contribution to Pc1 (stippled arrow). Minor more solitary patterns are shown with arrows. Andes reactivity was both observed within the PUUV pattern and as a minor more solitary pattern.

Figure 5. Distribution of hantavirus N protein megapeptides (bold face and boxed) in a hantavirus S segment phylogenetic tree. The tree was redrawn from that of Holmes and Zhang [Citation2]. Viruses known to give disease in humans are shown with a human figurine. Phylogroups I–IV and serogroups A–D [Citation26,Citation27] are indicated. Laibin virus (Phylogroup II) was recently discovered [Citation42]. The distinction between serogroups D1 and D2 was introduced in this paper. PCA ‘loadings’ indicate the degree of influence for each antigen in the respective principal component. They are here plotted according to phylogenetic relationship. Judging from the loadings of megapeptide antigens, Pc1 covers serogroups B and C while Pc2 covers serogroup A. However, Oxbow virus (Serogroup E) gave a minor contribution to Pc1 (stippled arrow). Minor more solitary patterns are shown with arrows. Andes reactivity was both observed within the PUUV pattern and as a minor more solitary pattern.

Aside from the two major patterns there were minor patterns that occurred outside of them. The patterns can be studied in detail in , where MFI of the two major patterns were normalized to PUUV and DOBV megapeptide reactivities, respectively ().

Figure 6. Detailed presentation of major and minor HSMIA patterns. Roman numerals refer to phylogroups, alphabetical characters to serogroups of Megapeptide sequences. X axis = contiguous hantavirus megapeptides (order as in and ) belonging to a serological pattern. Y axis: either relative MFI, in % of the strongest reacting megapeptide or MFI, in linear or logarithmic representation. The PUUV pattern is shown with standard deviations (SD), where the red line shows the average relative MFI + 1 SD, and the green line shows the average − 1 SD.

Figure 6. Detailed presentation of major and minor HSMIA patterns. Roman numerals refer to phylogroups, alphabetical characters to serogroups of Megapeptide sequences. X axis = contiguous hantavirus megapeptides (order as in Figures 1 and 4) belonging to a serological pattern. Y axis: either relative MFI, in % of the strongest reacting megapeptide or MFI, in linear or logarithmic representation. The PUUV pattern is shown with standard deviations (SD), where the red line shows the average relative MFI + 1 SD, and the green line shows the average − 1 SD.

Discussion

Serology has advantages as compared to nucleic acid (PCR)-based systems. Our SMIA system is simple to perform, can easily be automated, and uses only a few microliters of serum. Desirable extensions of the HSMIA work are testing of sera from all major continents, and adaptation for IgM determination. PCR is an important adjunct to HSMIA but is more laborious, requires extraction of a larger volume of sample, and may require sequencing for typing. Moreover, many hantavirus infections give a rather short period of viremia, which may not reach high copy numbers, see e.g. [Citation43,Citation44], thus serology is required for the diagnosis of acute hantavirus infections.

The disadvantages of HSMIA relative to EIA are: (i) It requires a Luminex flow meter, which is not present in all laboratories; and (ii) the color-coded beads are relatively expensive, a drawback in the early period of development when many antigens need to be tested. Otherwise, HSMIA offers multiplexity and serotyping, is more economical in terms of antigen (10–50 fold less is used per sample) and sample volume (only 1–5 uL per assay) use.

Hantaviruses as a target for surveillance and diagnosis

Hantaviruses cause zoonotic infections and are a serious global medical problem. Appropriate diagnostic systems are not present in many parts of the world. The great number of newly discovered hantaviruses, which may or may not give human infections, requires appropriate surveillance and diagnostic systems. We took advantage of the high conservation and antigenicity of the N terminus of the N protein. This, coupled with the facile production of megapeptides, allows an efficient detection of past infections of diverse hantaviruses. We are now evaluating the generality of this serology in a broader geographic context than analyzed here.

HSMIA and its design

We used a long spacer, which included His6, to ensure optimal presentation of N protein epitopes in the megapeptides. The N protein contains more or less group-specific epitopes distributed over its entire length [Citation20]. Truncated N protein antigens were previously used [Citation20,Citation45,Citation46]. It was previously suggested to be a suitable antigen for serotyping [Citation26]. We strived for inclusion of the most cross-reactive epitopes [Citation25Citation27,Citation45] that would allow detection of antibodies to widely differing hantaviruses. At least one representative of all four known phylogroups [Citation2] was therefore synthesized. Coupling efficiency was monitored throughout the megapeptide synthesis. The SMIA system is useful for large scale serosurveys in humans and animals because of the use of protein G, which can bind to a large variety of mammalian IgGs. Calculation of median fluorescence intensity (MFI) of 100 determinations gave a great precision, facilitating distinction of serological patterns. The background was low. Two of the 282 sera gave an MFI exceeding 200 (‘bead binders’). Such occasional non-specfic reactions [Citation47,Citation48] were reduced by background subtraction.

Relevance for detection and seroepidemiological surveillance of emerging hantaviruses

The HSMIA patterns conformed to the serotyping of the same sera by FRNT with DOBV and PUUV. The patterns largely followed phylogenetic relatedness. Patterns occurred within or between related serogroups. The PCA results were clear-cut, with two principal components (DOBV and PUUV) that were orthogonal to each other, i.e. with minimal overlap. The two patterns conformed largely to the FRNT and EIA results. The ability of the system of 16 megapeptides to distinguish at least two major seroreactivity patterns, of DOBV and PUUV infections, respectively, makes it likely that infections with other hantaviruses can also be detected and broadly typed with this system. Once a pattern has been identified, it can easily be added in the automatic pattern recognition algorithm. A multiplex method like SMIA, where all antigens compete for antibody binding in the same solution, inherently has a high precision. This should allow the definition of antibody reactivity patterns for many different hantaviral species, also for hitherto unknown ones. The influence of host variation, including histocompatibility antigens, on the patterns is unknown, but the use of a broad vertebrate IgG detection (Protein G) allows this to be rationally tested.

We have currently no explanation for the weak atypical PUUV pattern in the nucleoprotein-based HSMIA in the three sera that had PUUV specificity in FRNT. FRNT depends on other epitopes than those present in the N terminus of the nucleoprotein. The chosen nucleoprotein portion is conserved and variation in it is an unlikely cause of the result. Reassortment between virus genetic segments [Citation49], time- or HLA-dependent epitope recognition could be involved.

Minor patterns were also observed. Two, ANDV and SWSV, occurred in blood donors. They could represent spurious immunological cross-reactions occurring in a few humans, but may indicate that a few Uppsala blood donors have been infected with hitherto unknown or underdiagnosed hantaviruses, during international travel or inside Sweden. HSMIA is tailored for hantavirus surveillance, allowing attempts to correlate serological patterns to newly detected, rare or emerging viruses and to disease.

For example, SEOV was recently found in rats in the Netherlands and Sweden [Citation12,Citation50]. Finnish sera reactive with a SWSV antigen in a singleplex EIA, ascribed to PUUV cross-reactions were recently described [Citation51]. Our SWSV reactive sera had a distinct pattern with little PUUV megapeptide reactivity, indicating a SWSV-like origin of this pattern. This could be further studied using serology with less cross-reactive antigens, like glycoproteins, of SWSV, ASAV and ANDV.

Conclusions

The simplicity as well as comprehensive and easily extensible hantaviral coverage of HSMIA and its ability to measure antibodies in a wide range of mammals makes it an attractive system for global hantaviral surveillance. Megapeptide-based arrays should also be useful in other virus families with a mixture of linear and conformational non-glycosylation dependent epitopes. It is likely that PCA and HSMIA will aid in the detection of new hantaviral serological patterns and perhaps also in the detection of novel hantaviruses. It should be followed up with nucleic acid (PCR- or next generation-sequencing) and other serological tests, such as neutralization assays.

Acknowledgments

JB was granted laboratory space at the laboratory of Clinical Microbiology, Uppsala University Hospital. He was also supported by a research grant from the Uppsala Academic Hospital.

Disclosure statement

JB is the sole inventor behind: (i) a patent application covering the use of megapeptides for serology; and (ii) a patent application covering multiplex assays with enhanced specificity.

Additional information

Funding

This work was supported by the Academic Hospital, Uppsala.

References

  • Vaheri A , Henttonen H , Voutilainen L , et al. Hantavirus infections in Europe and their impact on public health. Rev Med Virol. 2013;23:1–11.
  • Holmes EC , Zhang Y-Z. The evolution and emergence of hantaviruses. Curr Opin Virol. 2015;10:27–33.
  • Heyman P , Vervoort T , Colson P , et al. A major outbreak of hantavirus infection in Belgium in 1995 and 1996. Epidemiol Infect. 1999;122:447–453.
  • Heyman P , Plyusnina A , Berny P , et al. Seoul hantavirus in Europe: first demonstration of the virus genome in wild Rattus norvegicus captured in France. Eur J Clin Microbiol Infect Dis. 2004;23:711–717.
  • Heyman P , Baert K , Plyusnina A , et al. Serological and genetic evidence for the presence of seoul hantavirus in Rattus norvegicus in Flanders, Belgium. Scand J Infect Dis. 2009;41:51–56.
  • Sanfeliu I , Nogueras MM , Gegúndez MI , et al. Seroepidemiological survey of hantavirus infection in healthy people in Valles Occidental, Barcelona. Vector Borne Zoonotic Dis. 2011;11:697–700.
  • Lin X-D , Guo W-P , Wang W , et al. Migration of Norway rats resulted in the worldwide distribution of Seoul hantavirus today. J Virol. 2012;86:972–981.
  • Plyusnina A , Heyman P , Baert K , et al. Genetic characterization of Seoul hantavirus originated from Norway rats (Rattus norvegicus) captured in Belgium. J Med Virol. 2012;84:1298–1303.
  • Clement J , Maes P , Van Ranst M. Hemorrhagic fever with renal syndrome in the new, and hantavirus pulmonary syndrome in the old world: paradi(se)gm lost or regained? Virus Res. 2014;187:55–58.
  • Christova I , Plyusnina A , Gladnishka T , et al. Detection of Dobrava hantavirus RNA in apodemus mice in Bulgaria. J Med Virol. 2015;87:263–268.
  • Goeijenbier M , Verner-Carlsson J , van Gorp EC , et al. Seoul hantavirus in brown rats in the Netherlands: implications for physicians – epidemiology, clinical aspects, treatment and diagnostics. Neth J Med. 2015;73:155–160.
  • Verner-Carlsson J , Lohmus M , Sundstrom K , et al. First evidence of Seoul hantavirus in the wild rat population in the Netherlands. Infect Ecol Epidemiol. 2015;5:27215.
  • Avsic-Zupanc T , Petrovec M , Duh D , et al. Puumala hantavirus in Slovenia: analyses of S and M segment sequences recovered from patients and rodents. Virus Res. 2007;123:204–210.
  • Voutilainen L , Sironen T , Tonteri E , et al. Life-long shedding of Puumala hantavirus in wild bank voles (Myodes glareolus). J Gen Virol. 2015;96:1238–1247.
  • Watson DC , Sargianou M , Papa A , et al. Epidemiology of hantavirus infections in humans: a comprehensive, global overview. Crit Rev Microbiol. 2014;40:261–272.
  • Manigold T , Vial P . Human hantavirus infections: epidemiology, clinical features, pathogenesis and immunology. Swiss Med Wkly. 2014;144:w13937.
  • Reusken C , Heyman P . Factors driving hantavirus emergence in Europe. Curr Opin Virol. 2013;3:92–99.
  • Hartline J , Mierek C , Knutson T , et al. Hantavirus infection in North America: a clinical review. Am J Emerg Med. 2013;31:978–982.
  • Mattar S , Guzmán C , Figueiredo LT . Diagnosis of hantavirus infection in humans. Expert Rev Anti Infect Ther. 2015;13:939–946.
  • Lundkvist A , Meisel H , Koletzki D , et al. Mapping of B-cell epitopes in the nucleocapsid protein of Puumala hantavirus. Viral Immunol. 2002;15:177–192.
  • Heiskanen T , Lundkvist A , Soliymani R , et al. Phage-displayed peptides mimicking the discontinuous neutralization sites of Puumala hantavirus envelope glycoproteins. Virology. 1999;262:321–332.
  • Elfaitouri A , Herrmann B , Bolin-Wiener A , et al. Epitopes of microbial and human heat shock protein 60 and their recognition in myalgic encephalomyelitis. PLoS One. 2013;8:e81155.
  • Blomberg J , Blomberg F , Sjösten A , et al. No evidence for xenotropic murine leukemia-related virus infection in Sweden using internally controlled multiepitope suspension array serology. Clin Vaccine Immunol. 2012;19:1399–1410.
  • Rizwan M , Ronnberg B , Cistjakovs M , et al. Serology in the digital age: using long synthetic peptides created from nucleic acid sequences as antigens in microarrays. Microarrays (Basel). 2016;5.
  • Elgh F , Linderholm M , Wadell G , et al. Development of humoral cross-reactivity to the nucleocapsid protein of heterologous hantaviruses in nephropathia epidemica. FEMS Immunol Med Microbiol. 1998;22:309–315.
  • Yoshimatsu K , Arikawa J . Serological diagnosis with recombinant N antigen for hantavirus infection. Virus Res. 2014;187:77–83.
  • Yoshimatsu K , Arikawa J . Antigenic properties of N protein of hantavirus. Viruses. 2014;6:3097–3109.
  • Zhang Q , Wang P , Kim Y , et al. Immune epitope database analysis resource (IEDB-AR). Nucleic Acids Res. 2008;36:W513–W518.
  • Wang Y , Boudreaux DM , Estrada DF , et al. NMR structure of the N-terminal coiled coil domain of the Andes hantavirus nucleocapsid protein. J Biol Chem. 2008;283:28297–28304.
  • Goeijenbier M , Hartskeerl RA , Reimerink J , et al. The hanta hunting study: underdiagnosis of Puumala hantavirus infections in symptomatic non-travelling leptospirosis-suspected patients in the Netherlands, in 2010 and April to November 2011. Euro Surveill. 2014;19:20878.
  • Gherasim A , Hjertqvist M , Lundkvist A , et al. Risk factors and potential preventive measures for nephropatia epidemica in Sweden 2011–2012: a case-control study. Infect Ecol Epidemiol. 2015;5:27698.
  • Lundkvist A , Kallio-Kokko H , Sjölander KB , et al. Characterization of Puumala virus nucleocapsid protein: identification of B-cell epitopes and domains involved in protective immunity. Virology. 1996;216:397–406.
  • Kallio-Kokko H , Leveelahti R , Brummer-Korvenkontio M , et al. Human immune response to Puumala virus glycoproteins and nucleocapsid protein expressed in mammalian cells. J Med Virol. 2001;65:605–613.
  • Hedman K , Vaheri A , Brummer-Korvenkontio M . Rapid diagnosis of hantavirus disease with an igg-avidity assay. Lancet. 1991;338:1353–1356.
  • Sjölander KB , Elgh F , Kallio-Kokko H , et al. Evaluation of serological methods for diagnosis of Puumala hantavirus infection (nephropathia epidemica). J Clin Microbiol. 1997;35:3264–3268.
  • Hepojoki S , Rusanen J , Hepojoki J , et al. Competitive homogeneous immunoassay for rapid serodiagnosis of hantavirus disease. J Clin Microbiol. 2015;53:2292–2297.
  • Klingström J , Plyusnin A , Vaheri A , et al. Wild-type Puumala hantavirus infection induces cytokines, C-reactive protein, creatinine, and nitric oxide in cynomolgus macaques. J Virol. 2002;76:444–449.
  • Lohmus M , Verner-Carlsson J , Borg O , et al. Hantavirus in new geographic regions, Sweden. Infect Ecol Epidemiol. 2016;6:31465.
  • Haste Andersen P , Nielsen M , Lund O . Prediction of residues in discontinuous B-cell epitopes using protein 3D structures. Protein Sci. 2006;15:2558–2567.
  • Vapalahti O , Kallio-Kokko H , Närvänen A , et al. Human B-cell epitopes of Puumala virus nucleocapsid protein, the major antigen in early serological response. J Med Virol. 1995;46:293–303.
  • Dzagurova T , Tkachenko E , Slonova R , et al. Antigenic relationships of hantavirus strains analysed by monoclonal antibodies. Arch Virol. 1995;140:1763–1773.
  • Xu L , Wu J , He B , et al. Novel hantavirus identified in black-bearded tomb bats, China. Infect Genet Evol. 2015;31:158–160.
  • Noh JY , Cheong HJ , Song JY , et al. Clinical and molecular epidemiological features of hemorrhagic fever with renal syndrome in Korea over a 10-year period. J Clin Virol. 2013;58:11–17.
  • Terajima M , Hendershot JD 3rd , Kariwa H , et al. High levels of viremia in patients with the hantavirus pulmonary syndrome. J Infect Dis. 1999;180:2030–2034.
  • Elgh F , Lundkvist A , Alexeyev OA , et al. A major antigenic domain for the human humoral response to Puumala virus nucleocapsid protein is located at the amino-terminus. J Virol Methods. 1996;59:161–172.
  • de Carvalho Nicacio C , Sällberg M , Hultgren C , et al. T-helper and humoral responses to Puumala hantavirus nucleocapsid protein: identification of T-helper epitopes in a mouse model. J Gen Virol. 2001;82:129–138.
  • Todd DJ , Knowlton N , Amato M , et al. Erroneous augmentation of multiplex assay measurements in patients with rheumatoid arthritis due to heterophilic binding by serum rheumatoid factor. Arthritis Rheum. 2011;63:894–903.
  • Waterboer T , Sehr P , Pawlita M . Suppression of non-specific binding in serological luminex assays. J Immunol Methods. 2006;309:200–204.
  • Razzauti M , Plyusnina A , Sironen T , et al. Analysis of Puumala hantavirus in a bank vole population in northern Finland: evidence for co-circulation of two genetic lineages and frequent reassortment between strains. J Gen Virol. 2009;90:1923–1931.
  • Lundkvist A , Verner-Carlsson J , Plyusnina A , et al. Pet rat harbouring Seoul hantavirus in Sweden, June 2013. Euro Surveill. 2013;18:20521.
  • Ling J , Vaheri A , Hepojoki S , et al. Serological survey of Seewis virus antibodies in patients suspected for hantavirus infection in Finland; a cross-reaction between Puumala virus antiserum with Seewis virus N protein? J Gen Virol. 2015;96:1664–1675.