Next Article in Journal
Prevention of CD8 T Cell Deletion during Chronic Viral Infection
Next Article in Special Issue
Geographic Distribution and Phylogeny of Soricine Shrew-Borne Seewis Virus and Altai Virus in Russia
Previous Article in Journal
Progress in Anti-Mammarenavirus Drug Development
Previous Article in Special Issue
Continuing Orthohantavirus Circulation in Deer Mice in Western Montana
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 7
10.3390/v13071188
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Serological Cross-Detection of Bat-Borne Hantaviruses: A Valid Strategy or Taking Chances?

by
Renata Carvalho de Oliveira
1,
Jorlan Fernandes
1,*,
Elba Regina de Sampaio Lemos
1,
Fernando de Paiva Conte
2,* and
Rodrigo Nunes Rodrigues-da-Silva
2
1
Laboratory of Hantaviroses and Rickettsioses, Oswaldo Cruz Institute, FIOCRUZ, Rio de Janeiro 21040-900, Brazil
2
Laboratory of Monoclonal Antibodies Technology, Immunobiological Technology Institute, FIOCRUZ, Rio de Janeiro 21040-360, Brazil
*
Authors to whom correspondence should be addressed.
Viruses 2021, 13(7), 1188; https://0-doi-org.brum.beds.ac.uk/10.3390/v13071188
Submission received: 30 April 2021 / Revised: 21 May 2021 / Accepted: 24 May 2021 / Published: 22 June 2021
(This article belongs to the Special Issue Hantavirus)
Browse Figures
Review Reports Versions Notes

Abstract

:
Bats are hosts of a range of viruses, and their great diversity and unique characteristics that distinguish them from all other mammals have been related to the maintenance, evolution, and dissemination of these pathogens. Recently, very divergent hantaviruses have been discovered in distinct species of bats worldwide, but their association with human disease remains unclear. Considering the low success rates of detecting hantavirus RNA in bat tissues and that to date no hantaviruses have been isolated from bat samples, immunodiagnostic tools could be very helpful to understand pathogenesis, epidemiology, and geographic range of bat-borne hantaviruses. In this sense, we aimed to identify in silico immunogenic B-cell epitopes present on bat-borne hantaviruses nucleoprotein (NP) and verify if they are conserved among them and other selected members of Mammantavirinae, using a combination of (the three most used) different prediction algorithms, ELLIPRO, Discotope 2.0, and PEPITO server. To support our data, we in silico modeled 3D structures of NPs from representative members of bat-borne hantaviruses, using comparative and ab initio methods due to the absence of crystallographic structures of studied proteins or similar models in the Protein Data Bank. Our analysis demonstrated the antigenic complexity of the bat-borne hantaviruses group, showing a low sequence conservation of epitopes among members of its own group and a minor conservation degree in comparison to Orthohantavirus, with a recognized importance to public health. Our data suggest that the use of recombinant rodent-borne hantavirus NPs to cross-detect antibodies against bat- or shrew-borne viruses could underestimate the real impact of this virus in nature.

1. Introduction

Bats (order Chiroptera) are taxonomically diverse, representing approximately 21% of all classified mammal species worldwide, with nearly 1422 species recognized [1]. The diversity of bats and their broad ecological features create many potential viral niches. In recent years, severe infectious diseases have been emerging, many of them caused by viruses originating from bats, such as Ebola virus, Marburg virus, Nipah virus, Hendra virus, and SARS coronaviruses [2]. Therefore, there is increasing evidence that bats are reservoirs for a high number of viruses, including hantaviruses [3].
Hantaviruses are important (re)emerging zoonotic pathogens previously classified in an unassigned order and in the genus Hantavirus of the family Bunyaviridae, recently reclassified and reassigned to a new family, designated Hantaviridae, in the order Bunyavirales [4,5]. The family Hantaviridae comprises four subfamilies: Actantavirinae, Agantavirinae, Mammantavirinae, and Repantavirinae. The subfamily Mammantavirinae has mammals, particularly rodents, as the most frequent natural hosts and is the largest and most complex Hantaviridae group, composed by four genera: Loanvirus, Mobatvirus, Orthohantavirus, and Thottimvirus [5]. Novel genotypes of hantaviruses have been identified in moles, shrews, and bats, expanding our understanding of the nature and origins of the hantaviruses [5,6,7,8]. Among the members of the subfamily Mammantavirinae, only Orthohantavirus is associated with human diseases: hemorrhagic fever with renal syndrome (HFRS) and hantavirus pulmonary syndrome, also called hantavirus cardiopulmonary syndrome (HPS or HCPS), respectively.
The presence of newly described hantaviruses in bats has challenged the conventional view that hantaviruses originated from rodents [8]. Several bat-borne hantaviruses are presently known, which show large genetic diversities from currently known rodent- and insectivore-borne (shrews and moles) hantaviruses [8,9,10]. Since the first report of hantavirus in Nycteris hispida and Neoromicia nanus bats in Africa in 2012, twelve bat-borne hantaviruses have been described in different regions of the world [3,11,12]. Of these, only Laibin mobatvirus (LAIV), Longquan loanvirus (LQUV), Quezon mobatvirus (QZNV), Brno loanvirus (BRNV), Robina orthohantavirus (ROBV), and Xuan Son mobatvirus (XSV) are currently recognized as species within genus Mobatvirus and Loanvirus by the International Committee on Taxonomy of Viruses (ICTV) released in 2020 [5]. Phylogenetic analysis of bat-borne hantavirus has indicated that bats might be the natural original hantavirus hosts. However, due to a lack of sufficient bat-borne hantavirus genomic sequences, their evolutionary phylogeny and genetic diversity as well as biological features are not fully understood [3,8,9,10,13,14].
Although bat-borne hantaviruses have been identified worldwide, the virus detection rate is low and restricted to limited locations, and it is still unclear if they cause infection or disease in humans. So far, none have been detected in the New World. Taking this into consideration, and the fact that none of these newfound bat-borne viruses have been isolated in cell culture in addition to the low success rates of detecting hantavirus RNA in bat tissues, improved diagnostic methods should be developed for virus detection [3,8]. In this context, serological tests could be important to uncover the real situation of bat-borne hantaviruses prevalence, help to understand their geographic distribution, and estimate the potential risk of these viruses to human health.
Nucleoprotein (NP) is the major antigen that elicits early serological responses in infected animals and has been used as a biomarker to develop antibodies for diagnostic and surveillance for several hantavirus species and genotypes [15,16,17]. To date, a few serological studies on bat-borne hantaviruses were conducted, using recombinant NP (rNP) of LAIV, XSV, and Seoul orthohantavirus (SEOV) as antigen for hantavirus detection in bat samples from China, and rNP of Araraquara orthohantavirus (ARAV rodent-borne hantaviruses) in bat samples from Brazil [18,19,20]. Therefore, this study aimed to verify if NPs B-cell epitopes of bat-borne hantaviruses and other selected members of Mammantavirinae are conserved among and between bat-, rodent-, and shrew-borne hantaviruses, using a combination of the three most used prediction algorithms, and to validate the epitopes predicted as specific to bat-borne viruses or as conserved among other hantaviruses.

2. Materials and Methods

2.1. Phylogenetic Analyses

The representative of complete NP sequences of hantavirus species as well as some unclassified hantaviruses [5] were downloaded from the GenBank® database of NCBI (www.ncbi.nlm.nih.gov/nuccore; accessed on 28 October 2020). Multiple sequence alignment (MSA) of complete regions was performed using MAFFT version 7 (mafft.cbrc.jp/alignment/server; accessed on 28 October 2020). The best-fit evolutionary model was determined using MEGA version X, using the Bayesian Information Criterion [21]. We estimated phylogenetic relations of the NP sequences using maximum likelihood inference, as implemented in PhyML 3 [22] under the LG + G + I model of sequence evolution. Statistical support of the clades was measured by a heuristic search with 1000 bootstrap replicates in PhyML [23].

2.2. Amino Acid Composition Analysis

The full-length NP amino acid sequences from representative hantaviruses, such as, LAIV (YP_009507251), XSV (AHA83416), BRNV (APU53639), Dakrong virus (DKGV-AXQ03852), Lena River (LENV-QBH68029), LQUV (YP_009664869), Nova mobatvirus (NVAV-YP_009362031), and QZNV (YP_009361846) were obtained from the GenBank® database of NCBI (www.ncbi.nlm.nih.gov/protein; accessed on 9 November 2020), as mentioned previously. Moreover, sequences of orthohantaviruses associated to rodent-borne diseases, including SEOV (ABM67564), Hantaan orthohantavirus (HTNV-ALA20831), Andes orthohantavirus (ANDV-AAK14323), ARAV (ABQ45402), and Sin Nombre (SNV-NP_941975) were involved in the study.
To analyze the amino acid composition of studied hantavirus NPs, we used the Composition Profiler Server (www.cprofiler.org/index.html; accessed on 9 November 2020), a web-based tool that automates detection of enrichment or depletion patterns of individual amino acids [24]. Hantavirus NPs amino acid compositions were compared using the SwissProt database [25]. Statistical significance associated with a specific enrichment or depletion is estimated using the two-sample t-test between two sequences of binary indicator variables, one sequence for each of the samples. A particular enrichment or depletion is statistically significant when the p-value (the lowest value at which the null hypothesis that the same underlying Gaussian distribution generated for both samples can be rejected) is lower than or equal to a user-specified statistical significance (alpha) value.

2.3. Hantavirus NP Amino Acid Variability Analysis

Sequences of studied hantavirus NPs were aligned using DNAStar Lasergene–Muscle Align (DNASTAR Inc., Madison, USA). The variability of amino acids in hantavirus NP sequences were evaluated using the Protein Variability Server (imed.med.ucm.es/PVS accessed on 17 November 2020). This server calculates the sequence variability within a multiple sequence alignment using Shannon entropy analysis [26]. Values of Shannon variability range from 0 (only one residue present at that position) to 4.322 (all 20 residues are equally represented in that position). Here, highly conserved positions are those with “Shannon variability” < 1.0.

2.4. Analysis of Secondary Structure and Solvent Accessibility

The analysis of secondary structure and solvent accessibility of LAIV, XSV, BRNV, DKGV, LENV, LQUV, NVAV, and QZNV (representative hantaviruses from phylogroup II, see Results Section) and orthohantavirus associated to rodent-borne diseases (SEOV, HTNV, ANDV, ARAV, and SNV) were evaluated using the PredictProtein server (www.predictprotein.org; accessed on 9 December 2020) [27]. This server uses RePROFsec to predict secondary structure elements, such as helix, beta-strand, and loop. Furthermore, this server uses PROFacc to predict solvent accessibility of protein residues for 10 states of relative accessibility that are grouped into two main states: buried and exposed.

2.5. Three-Dimensional (3D) Structure Modeling

The 3D structures of selected NP from phylogroup II hantavirus were modeled using the Robetta server (new.robetta.org; accessed on 17 December 2021). This server is continually evaluated through CAMEO (Continuous Automated Model Evaluation) [28,29]. The quality of the protein structures generated was analyzed by MolProbity (molprobity.biochem.duke.edu; accessed on 6 January 2021) in order to choose the best predictive model generated for each NP. MolProbity is a widely used system of model validation for protein and nucleic acid structures [30]. Its unique feature of all-atom contact analysis (including hydrogens) was described in 1999 [31], followed by its complementary rotamer, Ramachandran, and Cβ deviation criteria [32,33], and the initial MolProbity web service [34].

2.6. Conformational B-Cell Epitope Prediction

After 3D protein structure modeling, structure-based discontinuous epitope prediction was initially performed using the ELLIPRO server (tools.iedb.org/ellipro; accessed on 12 January 2021), since it correlates epitope-antigenicity, solvent accessibility, and flexibility in antigen proteins for predicting conformational B-cell epitopes [4]. In addition to ELLIPRO analysis, we also ran the Discotope 2.0 server (tools.iedb.org/discotope; accessed on 15 January 2021), that uses a combination of amino acid statistics, spatial information, and surface exposure [5], and PEPITO server (pepito.proteomics.ics.uci.edu; accessed on 20 January 2021), that incorporates an amino acid propensity scale along with side-chain orientation and solvent accessibility information using half-sphere exposure values to predict conformational epitopes [6]. Potential conformational epitopes were considered only those that were more than 20% of Ellipro epitope predicted by all three algorithms, using their default thresholds.

2.7. Conservancy Analysis of the Selected Epitopes

To evaluate the conservancies of phylogroup II hantavirus epitopes among Orthohantavirus NPs used in serological tests, we used an approach based on the recent study of Ulah and collaborators [35]. Briefly, the IEDB server (http://tools.iedb.org/conservancy/; accessed on 22 January 2021) was used to determine the conservancy of the discontinuous epitopes setting the sequence identity threshold at ‘≥70’. Conservancy of epitopes can be defined as the fraction of protein sequences that contain the epitopes at/or above a certain level of identity [36,37]. Here, predicted epitopes were compared to Orthohantavirus commonly used to serologically cross-detect antibodies against Mammantavirinae: SEOV (ABM67564), HTNV (ALA20831), ANDVs (AAK14323), ARAV (ABQ45402), and SNV (NP_941975).

3. Results

3.1. Phylogenetic Analyses

To construct phylogenetic relationships, 70 complete amino acid sequences of hantaviral NP were used. Currently, there are only seven complete NP sequences of bat-borne hantavirus species available in GenBank®. As shown in Figure S1, bat-, rodent-, and insectivore-borne hantaviruses showed a similar topology to the current taxonomic classification of Hantavidae family, here identified as clade I, II, III, and IV and V. Bat-borne hantaviruses were all clustered together within clade II, with two insectivore-borne hantaviruses, NVAV and LENV, respectively identified in Talpa moles and Sorex shrews, while other insectivore-borne hantaviruses were divided into clades I, II, III, IV, and VI. Only clade V contains solely rodent-borne hantaviruses, detected in the American continent.

3.2. NPs Vary in Length and Amino Acid Composition Profiles among Hantaviruses

In order to compare hantavirus NPs, Table 1 summarizes basic information about length and molecular weight of studied proteins in addition to ecological and taxonomic information about hantaviruses and their reservoirs. The length of the majority of hantavirus NPs ranges from 427 to 429 mers, however in phylogroup II, BRNV and LQUV were highlighted as shorter NPs, composed by 423 mers, while LENV presented the longer protein sequence, composed by 448 mers. Besides the length variation, hantavirus NPs presented similar molecular weight, ranging from ARAV 47,878.46 Da to LENV 50,532.14 Da.
Regarding the amino acid composition of hantavirus NPs, we observed a diverse composition profile among studied proteins (Figure 1). Using the Composition Profiler server, we compared the amino acids of hantavirus NP with the SwissProt database and identified enriched and depleted amino acids in each studied protein. In our study, only ANDV, ARAV, HTNV, and LAIV did not present any enriched or depleted amino acids. On the other hand, methionine residues (M) were enriched in SEOV (p = 0.01), NVAV (p = 0.01), QZNV, and DKGV (p = 0.01). DKGV also presented enriched leucine residues (p = 0.02), such as XSV (p = 0.03), ROBV (p = 0.03), and QZNV (p = 0.01), that additionally presented enriched aspartate residues (p = 0.01). Moreover, BRNV presented enriched glutamate (p = 0.01) and arginine (p = 0.05) residues and depleted alanine residues (p = 0.02), while LQUV presented enriched glutamine residues (p = 0.02) and depleted asparagine residues (p = 0.04). This asparagine residue was also depleted in LENV NP (p = 0.02), which yet presented enriched aspartate (p = 0.04) and glutamine (p = 0.04) residues. Besides, glycine residues were depleted in SNV (p = 0.02) and ROBV (p = 0.02) NPs.

3.3. C-Terminal Region of NPs Were More Conserved among the Studied Hantavirus

As illustrated in Figure 2, hantavirus NPs presented about 44% of amino acid sequence conservancy with the consensus. Moreover, the C-terminal region of NPs, presenting 65% of highly conserved amino acids, was more conserved among the studied hantavirus than the N-terminal region. C-terminal and N-terminal regions presented 65% and 48% of amino acids highly conserved, respectively.

3.4. Hantavirus NPs Share Similar Secondary Structures and Solvent Accessibility

Based on the amino acid sequences of studied proteins, we analyzed the secondary structure and solvent accessibility of hantavirus NPs. Despite the low conservation degree, we observed a similar profile of secondary structures and proportion of exposed amino acids in solvent. In all studied hantavirus NPs, we identified about 50% of amino acids within helix, 5% of amino acids as β-strand, and 45% as loop secondary structure (Figure S2). Regarding the accessibility of amino acids in solvent, in all studied proteins, 53% to 56% of amino acids were predicted as exposed in solvent, while the proportion of buried amino acids in solvent ranged from 39% to 41% (Figure S2).

3.5. Validation of Phylogroup II Hantavirus NPs Models

The quality of the modeled proteins’ structures was evaluated and validated using the MolProbity server, based on the Ramachandran plot, that described the amino acid positions in the plot as well as the overall [38], and in the complementary rotamer, Ramachandran, and Cβ deviation criteria [32,33]. Here, in all selected models, the plot showed that more than 97.2% of amino acids were arranged in Ramachandran favored zones and presented more than 98.7% favored rotamers. Moreover, predicted models for LAIV, NVAV, QXNV, XSV, and DKGV NPs did not present outliers in the Ramachandran plot, only one was presented for BRNV and LQUV and two for ROBV and LENV. In the same way, among the used hantavirus NPs models, the LAIV model presented one poor rotamer, while all other models did not present poor rotamers (Table 2). These data showed that the generated structures of phylogroup II hantavirus NPs were good considering the overall geometry and can be used for further study.

3.6. Predicted Conformational B-Cell Epitopes in Hantavirus NPs Models

After validating and evaluating the quality of phylogroup II hantavirus NPs models and aiming to predict its main B-cell epitopes, we initially used Ellipro to predict discontinuous epitopes in combination with Discotope 2.0 and Pepito to refine the analysis. In this way, sixty epitopes were predicted by Ellipro, however, only 55% (33 epitopes) of these epitopes were also predicted by Discotope 2.0 and Pepito (Table S1). Conformational epitopes lengths ranged from 4 to 118 mers. As showed in Figure 3, LQUV, LAI, QZNV, and ROBV NPs presented three predicted epitopes, while four epitopes were predicted in BRNV, NVAV, XSV, and DKGV models, and five were predicted in the LENV virus NP model. Moreover, all predicted epitopes are located in exposed regions of hantavirus NPs and thus can be recognized by specific antibodies (Figure 3b,d,f,h,j,l,n,p,r).

3.7. Most of Conformational B-Cell Epitopes Are Not Conserved among Hantavirus

Aiming to verify the conservancy of predicted B-cell epitopes among hantaviruses, they were compared to phylogroup II hantavirus NPs (BRNV, LQUV, LAIV, NVAV, QZNV, XSV, LENV, DKGV, and ROBV). Besides, considering the use of Orthohantavirus NPs for serological detection of bat-borne hantavirus, these epitopes were also compared to ANDV, ARAV, SEOV, SNV, and HTNV NPs. Despite the high conservancy observed among the C-terminal region of hantavirus NPs, most of the predicted NPs B-cell epitopes were not conserved among phylogroup II hantavirus and Orthohantavirus enrolled in this study (Table 3). Focusing on the epitopes’ conservancy among phylogroup II hantavirus, only XSV and DKGV viruses presented 50% (2/4) of their epitopes (XS2, XS3, D1, and D4), showing a high average degree of conservation (≥70%) among the majority of representative phylogroup II hantaviruses, while BRNV, LAIV, QZNV, LENV, and ROBV NPs presented only one epitope highly conserved among this phylogroup. Moreover, BRNV, LAIV, XSV, and LENV viruses NPs presented one epitope (B4, LB3, XS3, and LS4, respectively) highly conserved (identity ≥ 70%) among 93% (13/15) of studied hantaviruses (Table 4). On the other hand, LQUV and NVAV did not present any highly conserved epitope among studied viruses.

4. Discussion

Several bat-borne hantaviruses have recently been identified, but to date, it is not clear whether they are pathogenic for humans. Thus far, a total of 12 hantaviruses were identified in 16 different species of Old-World bats in Africa, Asia, Europe, and Australia [3,39]. The discovery of highly divergent lineages of hantaviruses in bats of different species and their vast geographic distribution provides unlimited opportunities to search for other bat-associated hantaviruses, including in the New World. However, the low success rates of detecting hantavirus RNA in bat tissues observed in previous studies and the difficulty to isolate hantavirus from bat samples [3] reinforce the need to develop alternative and more sensitive diagnostic tools to help in understanding the role of hantaviruses harbored by bats on disease, epidemiology, and their geographic range.
In this sense, immunodiagnostic tools could be very helpful to uncover the real scenario of bat-borne hantaviruses, that is still poorly understood. Remarkably, NP is the major hantavirus antigen that elicits early humoral immune responses and has been used as a serological target and for antigen detection in animals and patients [40,41,42]. This protein contains both serotype-specific and common group epitopes [43,44,45] and has been used to investigate the presence of hantaviruses in bats [18] and shrews [46] and to screen bats that could harbor hantaviruses in the New World [19,20]. However, there are no studies comparing the similarity of B-cell epitopes among bat- and rodent-borne hantaviruses.
Here, our data suggest that the general tertiary structures of NP seem to be similar among members of the subfamily Mammantavirinae, since all studied proteins presented similar frequencies of exposed regions and secondary structures (Figure 2). Nevertheless, NPs of bat- and shrew-borne hantaviruses from phylogroup II were more variable in length and molecular weight than the ones belonging to Orthohantavirus, with lengths ranging from 423 to 448 mers and about 48,630 Da of mean molecular weight, while Orthohantavirus NP presented about 428 amino acids of length and 48,035 Da of mean molecular weight (Table 1). Moreover, we observed non-conserved profiles of amino acid composition among studied hantaviruses NPs, in which only ANDV, ARAV, HNTV, and LAIV did not present any enriched or depleted amino acid residue (Figure 1). These data are in agreement with the analysis of conservation degree among studied proteins that showed about 51% of amino acids’ sequence conservation. Epitopes can be distributed over the whole NP, but it is thought that the C-terminal half of NP contains rather conformation- and multimerization-dependent epitopes, which should be more serotype-specific [47,48]. In a recently published work from our group, we demonstrated that NP amino acid sequence conservation among Orthohantavirus associated to HFRS (SEOV, HNTV, Amur virus, Dobrava-Belgrade orthohantavirus, and Puumala orthohantavirus) and associated to HCPS (ANDV, SNV, and Laguna Negra orthohantavirus) ranged from 62% to 95% [45]. The N-terminal quarter of hantaviral NP bears linear and immunodominant epitopes, but as seen here, a possible and interesting antigenic site was found in the C-terminal half [47]. Thus, we believe that the observed low conservation among NPs from phylogroup II hantavirus in comparison with NPs from representative Orthohantavirus associated to rodent-borne disease corroborates the view that identification of common immunodominant epitopes in these proteins is critical for a better selection of targets to serological detection of antibodies.
Recently, Xu and collaborators [18], using recombinant NP from phylogroup II hantaviruses (LAIV and XSV) and SEOV, observed the cross-reaction of XSV-infected bat serum samples with SEOV rNP, while samples of bats infected by LAIV were not able to cross-react to SEOV rNP [18]. Indirectly, this finding makes us question whether the use of Orthohantavirus nucleoproteins cannot underestimate the presence of antibodies against bat-borne hantavirus, especially in the New World, where the knowledge about bat-borne hantavirus species and their reservoirs remains scarce. Considering this, our aim in this study was to predict immunodominant B-cell epitopes in phylogroup II hantavirus NPs and to compare them among phylogroup II and also among representative Orthohantavirus associated to rodent-borne diseases classically used in serological antibody detection.
B-cell epitope prediction aims to facilitate the identification of antigenic regions of proteins, that could be used to replace antigens for antibody production and to accelerate the development of novel vaccines and diagnostic/serologic tools [49]. In this context, conformational B-cell epitopes correspond to about 90% of protein surface-accessible clusters of amino acids able to be recognized by secreted antibodies or B-cell receptors eliciting humoral immune response [50]. Unfortunately, its prediction is limited by the knowledge of protein 3D structure, information only available for a few proteins [51], and absolutely scarce to hantavirus proteins. Despite the absence of crystallographic structures of studied proteins or similar models in the Protein Data Bank [52], we in silico modeled 3D structures of NPs from representative members of phylogroup II hantaviruses (BRNV, LQUV, LAIV, NVAV, QZNV, XSV, LENV, DKGV, and ROBV), using comparative and ab initio methods. All generated structures were well-evaluated (Table 2) and considered good models for conformational B-cell epitopes’ predictions [30,53,54,55].
Thus, in each studied NP, between three to five clusters of amino acids were predicted as exposed discontinuous B-cell epitopes, composed by 4 to 118 mers in each cluster (Figure 3). From our point of view, the similar number of predicted epitopes in studied proteins corroborates the similar profiles of exposed amino acids among NP hantavirus (data not shown). Moreover, the number of amino acids in predicted epitopes ranged from 30% to 44% of the protein length, and we believe that this variation was acceptable due to differences in protein length (Table 1) and 3D structure models. Besides, our data corroborates the study of Kalaiselvan and collaborators, who reported a similar number of B-cell linear epitopes among Orthohantavirus associated to HFRS [43]; however, this was the first study aiming to identify and to compare B-cell epitopes among phylogroup II hantaviruses.
Regarding the conservation of hantavirus NP-exposed regions, our observations showed that the majority of predicted epitopes are not conserved among hantavirus phylogroup II, since only 9 (27%) of 33 predicted epitopes share more than 70% of their amino acids with at least half of the selected representative members of this group. Interestingly, among studied representative members of phylogroup II hantavirus, LQUV and NVAV did not present conserved epitopes with other members of this group, while XSV and DKGV viruses presented two out of four highly conserved predicted epitopes. Remarkably, our data suggested that XSV NP is antigenically closer to other phylogroup II hantaviruses, since it shares three highly conserved epitopes with LQUV and LAIV, two epitopes with QZNV, DKGV, and ROBV viruses, and one epitope with BRNV and with hantaviruses harbored by shrews (NVAV and LENV). Considering the geospatial context, the geographic distribution of their reservoir hosts and their overlapping areas with other bat species (where hantavirus host-switching likely occurs), could also be taken into account in the serological surveillance of hantavirus regionally. XSV and DKGV were reported in bat species that belong to the family Hipposideridae, one of the most speciose of insectivorous bats. The reservoir hosts of XSV, the Pomona roundleaf bats (Hipposideros pomona) and Ashy roundleaf bats (Hipposideros cineraceus), have a higher geographic distribution area in comparison to the Stoliczka’s Asian trident bats (Aselliscus stoliczkanus), the reservoir host of DKGV. The vast geographic distribution of both of these Hipposideros species in Asia, that are sympatric and often roost in the same caves, and the distribution of 70 other species members of this large genus, provide opportunities to detect XSV-related hantaviruses or other genetically divergent mobatviruses throughout Africa, Europe, Asia, and Australia [3]. To date, only one study used two different recombinant proteins of phylogroup II hantavirus for serological detection of hantaviruses in bats, demonstrating a high reactivity among studied samples against XSV and LAIV recombinant NPs [18].
On the other side, recombinant Orthohantavirus NPs have been used to screen antibodies against bat-borne hantaviruses. Recently, recombinant SEOV NP was used to cross-detect antibodies against LAIV and XSV [18] and recombinant ARAV NP was used to investigate whether phyllostomid bats could harbor hantaviruses in Brazil [20]. Our data demonstrated that phylogroup II hantavirus NPs are antigenically low-conserved among Orthohantavirus (Table 4), presenting no more than one highly conserved B-cell epitope shared with the majority of compared Orthohantavirus (ANDV, ARAV, SNV, SEOV, and HTNV). Based on our findings, the use of recombinant Orthohantavirus NPs to cross-detect antibodies against bat-borne hantaviruses could lead to an underestimation of the real reactivity, resulting in low-sensibility strategies, since only seven (21%) of phylogroup II hantavirus NP predicted epitopes were highly conserved among Orthohantavirus NP. Corroborating with the findings of Xu and collaborators [18], who reported cross-reactivity among XSV and SEOV recombinant NP but observed no cross-reactivity against LAIV, another bat-borne hantavirus, and SEOV NPs. Here, we demonstrated that XSV NP shares two epitopes with SEOV, while LAIV shares only one short epitope, explaining the absence of cross-reactivity. Moreover, similar results were found by Wei at al. in their studies with shrew-borne hantavirus, which suggested low or no cross-reactivity among Seewis orthohantavirus, Altai orthohantavirus, Thottapalayam thottimvirus, Asama orthohantavirus, and rodent-borne hantaviruses [56]. Thereby, we believe that the knowledge about the conservation of natural epitopes is critical for serological diagnosis based on cross-reactions, especially considering the limited information about the spread of bat-borne hantaviruses around the world [3,39,57]. This suggestion is supported by Tischler’s study, that showed that ANDV NP presented different humoral key targets to humans and to rodents [58], reinforcing the critical role of conserved epitopes to higher chances of cross-reaction.
Efforts to develop novel and accurate diagnostic tools able to detect bat-borne hantaviruses are urgent to allow a better comprehension of its spread and its real impact on animal and human health. Here, we demonstrated the antigenic complexity of bat-borne hantaviruses, showing a low sequence conservation of epitopes among members of its own group and a minor conservation degree in comparison with rodent- and shrew-borne hantaviruses. Our data suggest that the use of recombinant rodent-borne hantavirus NPs to cross-detect antibodies against bat- or shrew-borne viruses could underestimate the real impact of this virus in nature. Diagnostic use of homologous virus antigen in the assay seems to be preferred to search for antibodies after infection by a particular hantavirus. For a broader detection, among studied proteins, XSV NP presented the higher sequence conservation among phylogroup II hantavirus members, and we believe that should be considered the better choice for serological cross-detection of bat-borne hantavirus. Besides, immunodominant linear B-cell epitopes of bat-borne hantavirus NPs can allow the development of novel and specific diagnostic approaches, but this strategy still needs to be explored by accurately identifying antibody targets.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/v13071188/s1, Figure S1: Phylogenetic tree-based complete hantavirus NP sequences segment; Figure S2: Analysis of hantavirus NPs structures; Table S1: Epitopes predicted by Ellipro, Discotope 2.0, and Pepito.

Author Contributions

Conceptualization, R.C.d.O., J.F. and E.R.d.S.L.; methodology and data curation, F.d.P.C. and R.N.R.-d.-S.; writing—original draft preparation, R.C.d.O., J.F., F.d.P.C. and R.N.R.-d.-S.; writing—review and editing, R.C.d.O., J.F., F.d.P.C., R.N.R.-d.-S. and E.R.d.S.L.; funding, E.R.d.S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazil (CAPES), grant number 001, and the APC was funded by Fundação Oswaldo Cruz (Fiocruz)—Brazil.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within article and supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ASM Mammal Diversity Database. Available online: Mammaldiversity.org (accessed on 27 July 2020).
  2. Wang, L.F.; Anderson, D.E. Viruses in bats and potential spillover to animals and humans. Curr. Opin. Virol. 2019, 34, 79–89. [Google Scholar] [CrossRef] [PubMed]
  3. Arai, S.; Yanagihara, R. Genetic Diversity and Geographic Distribution of Bat-borne Hantaviruses. Curr. Issues Mol. Biol. 2020, 39, 1–28. [Google Scholar] [CrossRef]
  4. Guterres, A.; Oliveira, R.C.; Fernandes, J.; de Lemos, E.R.; Schrago, C.G. New bunya-like viruses: Highlighting their relations. Infect. Genet. Evol. 2017, 49, 164–173. [Google Scholar] [CrossRef]
  5. Kuhn, J.H.; Adkins, S.; Alioto, D.; Alkhovsky, S.V.; Amarasinghe, G.K.; Anthony, S.J.; Avšič-Županc, T.; Ayllón, M.A.; Bahl, J.; Balkema-Buschmann, A.; et al. 2020 taxonomic update for phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales. Arch. Virol. 2020, 165, 3023–3072. [Google Scholar] [CrossRef]
  6. De Oliveira, R.C.; Guterres, A.; Fernandes, J.; D’Andrea, P.S.; Bonvicino, C.R.; de Lemos, E.R. Hantavirus reservoirs: Current status with an emphasis on data from Brazil. Viruses 2014, 6, 1929–1973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Kang, H.J.; Gu, S.H.; Yashina, L.N.; Cook, J.A.; Yanagihara, R. Highly Divergent Genetic Variants of Soricid-Borne Altai Virus (Hantaviridae) in Eurasia Suggest Ancient Host-Switching Events. Viruses 2019, 11, 857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Yanagihara, R.; Gu, S.H.; Arai, S.; Kang, H.J.; Song, J.W. Hantaviruses: Rediscovery and new beginnings. Virus Res. 2014, 187, 6–14. [Google Scholar] [CrossRef] [Green Version]
  9. Bennett, S.N.; Gu, S.H.; Kang, H.J.; Arai, S.; Yanagihara, R. Reconstructing the evolutionary origins and phylogeography of hantaviruses. Trends Microbiol. 2014, 22, 473–482. [Google Scholar] [CrossRef] [Green Version]
  10. Witkowski, P.T.; Drexler, J.F.; Kallies, R.; Ličková, M.; Bokorová, S.; Maganga, G.D.; Szemes, T.; Leroy, E.M.; Krüger, D.H.; Drosten, C.; et al. Phylogenetic analysis of a newfound bat-borne hantavirus supports a laurasiatherian host association for ancestral mammalian hantaviruses. Infect. Genet. Evol. 2016, 41, 113–119. [Google Scholar] [CrossRef]
  11. Sumibcay, L.; Kadjo, B.; Gu, S.H.; Kang, H.J.; Lim, B.K.; Cook, J.A.; Song, J.W.; Yanagihara, R. Divergent lineage of a novel hantavirus in the banana pipistrelle (Neoromicia nanus) in Cote d’Ivoire. Virol. J. 2012, 9, 34. [Google Scholar] [CrossRef] [Green Version]
  12. Weiss, S.; Witkowski, P.T.; Auste, B.; Nowak, K.; Weber, N.; Fahr, J.; Mombouli, J.V.; Wolfe, N.D.; Drexler, J.F.; Drosten, C.; et al. Hantavirus in Bat, Sierra Leone. Emerg. Infect. Dis. 2012, 18, 159–161. [Google Scholar] [CrossRef] [PubMed]
  13. Guo, W.P.; Lin, X.D.; Wang, W.; Tian, J.H.; Cong, M.L.; Zhang, H.L.; Wang, M.R.; Zhou, R.H.; Wang, J.B.; Li, M.H.; et al. Phylogeny and origins of hantaviruses harbored by bats, insectivores, and rodents. PLoS Pathog. 2013, 9, e1003159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Arai, S.; Kikuchi, F.; Bawm, S.; Sơn, N.T.; Lin, K.S.; Tú, V.T.; Aoki, K.; Tsuchiya, K.; Tanaka-Taya, K.; Morikawa, S.; et al. Molecular Phylogeny of Mobatviruses (Hantaviridae) in Myanmar and Vietnam. Viruses 2019, 11, 228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Vapalahti, O.; Kallio-Kokko, H.; Narvanen, A.; Julkunen, I.; Lundkvist, A.; Plyusnin, A.; Lehväslaiho, H.; Brummer-Korvenkontio, M.; Vaheri, A.; Lankinen, H. Human B-cell epitopes of Puumala virus nucleocapsid protein, the major antigen in early serological response. J. Med. Virol. 1995, 46, 293–303. [Google Scholar] [CrossRef]
  16. Padula, P.J.; Rossi, C.M.; Valle, M.O.D.; Martínez, P.V.; Colavecchia, S.B.; Edelstein, A.; Miguel, S.D.L.; Rabinovich, R.D.; Segura, E.L. Development and evaluation of a solid-phase enzyme immunoassay based on Andes hantavirus recombinant nucleoprotein. J. Med. Microbiol. 2000, 49, 149–155. [Google Scholar] [CrossRef] [Green Version]
  17. Olal, D.; Daumke, O. Structure of the Hantavirus Nucleoprotein Provides Insights into the Mechanism of RNA Encapsidation. Cell Rep. 2016, 14, 2092–2099. [Google Scholar] [CrossRef] [Green Version]
  18. Xu, L.; Wu, J.; Li, Q.; Wei, Y.; Tan, Z.; Cai, J.; Guo, H.; Yang, L.; Huang, X.; Chen, J.; et al. Seroprevalence, cross antigenicity and circulation sphere of bat-borne hantaviruses revealed by serological and antigenic analyses. PLoS Pathog. 2019, 15, e1007545. [Google Scholar] [CrossRef] [PubMed]
  19. Sabino-Santos, G., Jr.; Ferreira, F.F.; Da Silva, D.J.F.; Machado, D.M.; Da Silva, S.G.; Sao Bernardo, C.S.; Dos Santos Filho, M.; Levi, T.; Figueiredo, L.T.M.; Peres, C.A.; et al. Hantavirus antibodies among phyllostomid bats from the arc of deforestation in Southern Amazonia, Brazil. Transbound. Emerg. Dis. 2020, 67, 1045–1051. [Google Scholar] [CrossRef] [PubMed]
  20. Sabino-Santos, G., Jr.; Maia, F.G.; Vieira, T.M.; De Lara Muylaert, R.; Lima, S.M.; Goncalves, C.B.; Barroso, P.D.; Melo, M.N.; Jonsson, C.B.; Goodin, D.; et al. Evidence of Hantavirus Infection Among Bats in Brazil. Am. J. Trop. Med. Hyg. 2015, 93, 404–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  22. Guindon, S.; Gascuel, O. A Simple, Fast, and Accurate Algorithm to Estimate Large Phylogenies by Maximum Likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef] [Green Version]
  23. Anisimova, M.; Gascuel, O. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst. Biol. 2006, 55, 539–552. [Google Scholar] [CrossRef] [PubMed]
  24. Vacic, V.; Uversky, V.N.; Dunker, A.K.; Lonardi, S. Composition Profiler: A tool for discovery and visualization of amino acid composition differences. BMC Bioinform. 2017, 8, 211. [Google Scholar] [CrossRef] [Green Version]
  25. Boeckmann, B.; Blatter, M.C.; Famiglietti, L.; Hinz, U.; Lane, L.; Roechert, B.; Bairoch, A. Protein variety and functional diversity: Swiss-Prot annotation in its biological context. C. R. Biol. 2005, 328, 882–899. [Google Scholar] [CrossRef] [PubMed]
  26. Shannon, C.E. The mathematical theory of communication. 1963. MD Comput. 1997, 14, 306–317. [Google Scholar]
  27. Rost, B.; Yachdav, G.; Liu, J. The PredictProtein server. Nucleic Acids Res. 2004, 32, W321–W326. [Google Scholar]
  28. Kim, D.E.; Chivian, D.; Baker, D. Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004, 32, W526–W531. [Google Scholar]
  29. Song, Y.; Dimaio, F.; Wang, R.Y.; Kim, D.; Miles, C.; Brunette, T.; Thompson, J.; Baker, D. High-resolution comparative modeling with Rosetta CM. Structure 2013, 21, 1735–1742. [Google Scholar]
  30. Williams, C.J.; Headd, J.J.; Moriarty, N.W.; Prisant, M.G.; Videau, L.L.; Deis, L.N.; Verma, V.; Keedy, D.A.; Hintze, B.J.; Chen, V.B.; et al. MolProbity: More and better reference data for improved all-atom structure validation. Protein Sci. 2018, 27, 293–315. [Google Scholar] [CrossRef] [PubMed]
  31. Word, J.M.; Lovell, S.C.; Labean, T.H.; Taylor, H.C.; Zalis, M.E.; Presley, B.K.; Richardson, J.S.; Richardson, D.C. Visualizing and quantifying molecular goodness-of-fit: Small-probe contact dots with explicit hydrogen atoms. J. Mol. Biol. 1999, 285, 1711–1733. [Google Scholar] [CrossRef] [PubMed]
  32. Lovell, S.C.; Word, J.M.; Richardson, J.S.; Richardson, D.C. The penultimate rotamer library. Proteins 2000, 40, 389–408. [Google Scholar] [CrossRef]
  33. Prisant, M.G.; Richardson, J.S.; Richardson, D.C. Structure validation by Calpha geometry: Phi, psi and Cbeta deviation. Proteins 2003, 50, 437–450. [Google Scholar]
  34. Davis, I.W.; Murray, L.W.; Richardson, J.S.; Richardson, D.C. MOLPROBITY: Structure validation and all-atom contact analysis for nucleic acids and their complexes. Nucleic Acids Res. 2004, 32, W615–W619. [Google Scholar] [CrossRef] [PubMed]
  35. Ullah, M.A.; Sarkar, B.; Islam, S.S. Exploiting the reverse vaccinology approach to design novel subunit vaccines against Ebola virus. Immunobiology 2020, 225, 151949. [Google Scholar] [CrossRef] [PubMed]
  36. Bui, H.H.; Sidney, J.; Li, W.; Fusseder, N.; Sette, A. Development of an epitope conservancy analysis tool to facilitate the design of epitope-based diagnostics and vaccines. BMC Bioinform. 2007, 8, 361. [Google Scholar] [CrossRef] [Green Version]
  37. Vita, R.; Mahajan, S.; Overton, J.A.; Dhanda, S.K.; Martini, S.; Cantrell, J.R.; Wheeler, D.K.; Sette, A.; Peters, B. The Immune Epitope Database (IEDB): 2018 update. Nucleic Acids Res. 2019, 47, D339–D343. [Google Scholar] [CrossRef] [Green Version]
  38. Laskowski, R.A.; Moss, D.S.; Thornton, J.M. Main-chain bond lengths and bond angles in protein structures. J. Mol. Biol. 1993, 231, 1049–1067. [Google Scholar] [CrossRef]
  39. Zana, B.; Kemenesi, G.; Buzas, D.; Csorba, G.; Gorfol, T.; Khan, F.A.; Tahir, N.; Zeghbib, S.; Madai, M.; Papp, H.; et al. Molecular Identification of a Novel Hantavirus in Malaysian Bronze Tube-Nosed Bats (Murina aenea). Viruses 2019, 11, 887. [Google Scholar] [CrossRef] [Green Version]
  40. Raboni, S.M.; Levis, S.; Rosa, E.S.; Bisordi, I.; Delfraro, A.; Lemos, E.; Correia, D.C.; Duarte Dos Santos, C.N. Hantavirus infection in Brazil: Development and evaluation of an enzyme immunoassay and immunoblotting based on N recombinant protein. Diagn. Microbiol. Infect. Dis. 2007, 58, 89–97. [Google Scholar] [CrossRef]
  41. Schlegel, M.; Tegshduuren, E.; Yoshimatsu, K.; Petraityte, R.; Sasnauskas, K.; Hammerschmidt, B.; Friedrich, R.; Mertens, M.; Groschup, M.H.; Arai, S.; et al. Novel serological tools for detection of Thottapalayam virus, a Soricomorpha-borne hantavirus. Arch. Virol. 2012, 157, 2179–2187. [Google Scholar] [CrossRef]
  42. Amada, T.; Yoshimatsu, K.; Yasuda, S.P.; Shimizu, K.; Koma, T.; Hayashimoto, N.; Gamage, C.D.; Nishio, S.; Takakura, A.; Arikawa, J. Rapid, whole blood diagnostic test for detecting anti-hantavirus antibody in rats. J. Virol. Methods 2013, 193, 42–49. [Google Scholar] [CrossRef]
  43. Kalaiselvan, S.; Sankar, S.; Ramamurthy, M.; Ghosh, A.R.; Nandagopal, B.; Sridharan, G. Prediction of B Cell Epitopes Among Hantavirus Strains Causing Hemorragic Fever with Renal Syndrome. J. Cell Biochem. 2017, 118, 1182–1188. [Google Scholar] [CrossRef]
  44. Kalaiselvan, S.; Sankar, S.; Ramamurthy, M.; Ghosh, A.R.; Nandagopal, B.; Sridharan, G. Prediction of Pan-Specific B-Cell Epitopes from Nucleocapsid Protein of Hantaviruses Causing Hantavirus Cardiopulmonary Syndrome. J. Cell Biochem. 2017, 118, 2320–2324. [Google Scholar] [CrossRef]
  45. Conte, F.P.; Tinoco, B.C.; Santos Chaves, T.; Oliveira, R.C.; Figueira Mansur, J.; Mohana-Borges, R.; Lemos, E.R.S.; Neves, P.; Rodrigues-Da-Silva, R.N. Identification and validation of specific B-cell epitopes of hantaviruses associated to hemorrhagic fever and renal syndrome. PLoS Negl. Trop. Dis. 2019, 13, e0007915. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Okumura, M.; Yoshimatsu, K.; Kumperasart, S.; Nakamura, I.; Ogino, M.; Taruishi, M.; Sungdee, A.; Pattamadilok, S.; Ibrahim, I.N.; Erlina, S.; et al. Development of serological assays for Thottapalayam virus, an insectivore-borne Hantavirus. Clin. Vaccine Immunol. 2007, 14, 173–181. [Google Scholar] [CrossRef] [Green Version]
  47. Yoshimatsu, K.; Arikawa, J. Serological diagnosis with recombinant N antigen for hantavirus infection. Virus Res. 2014, 187, 77–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Reuter, M.; Krüger, D.H. The nucleocapsid protein of hantaviruses: Much more than a genome-wrapping protein. Virus Genes 2018, 54, 5–16. [Google Scholar] [CrossRef] [PubMed]
  49. Sanchez-Trincado, J.L.; Gomez-Perosanz, M.; Reche, P.A. Fundamentals and Methods for T- and B-Cell Epitope Prediction. J. Immunol. Res. 2017, 2680160. [Google Scholar] [CrossRef] [Green Version]
  50. Van Regenmortel, M.H. What is a B-cell epitope? Methods Mol. Biol. 2009, 524, 3–20. [Google Scholar]
  51. Levitt, M. Nature of the protein universe. Proc. Natl. Acad. Sci. USA 2009, 106, 11079–11084. [Google Scholar] [CrossRef] [Green Version]
  52. Sussman, J.L.; Lin, D.; Jiang, J.; Manning, N.O.; Prilusky, J.; Ritter, O.; Abola, E.E. Protein Data Bank (PDB): Database of three-dimensional structural information of biological macromolecules. Acta Crystallogr. D Biol. Crystallogr. 1998, 54, 1078–1084. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Eisenberg, D.; Luthy, R.; Bowie, J.U. VERIFY3D: Assessment of protein models with three-dimensional profiles. Methods Enzymol. 1997, 277, 396–404. [Google Scholar] [PubMed]
  54. Von Grotthuss, M.; Pas, J.; Wyrwicz, L.; Ginalski, K.; Rychlewski, L. Application of 3D-Jury, GRDB, and Verify3D in fold recognition. Proteins 2003, 53 (Suppl. S6), 418–423. [Google Scholar] [CrossRef] [PubMed]
  55. Prisant, M.G.; Williams, C.J.; Chen, V.B.; Richardson, J.S.; Richardson, D.C. New tools in MolProbity validation: CaBLAM for CryoEM backbone, UnDowser to rethink "waters," and NGL Viewer to recapture online 3D graphics. Protein Sci. 2020, 29, 315–329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Wei, Z.; Shimizu, K.; Nishigami, K.; Tsuda, Y.; Sarathukumara, Y.; Muthusinghe, D.S.; Gamage, C.D.; Granathne, L.; Lokupathirage, S.M.W.; Nanayakkara, N.; et al. Serological methods for detection of infection with shrew-borne hantaviruses: Thottapalayam, Seewis, Altai, and Asama viruses. Arch. Virol. 2021, 166, 275–280. [Google Scholar] [CrossRef] [PubMed]
  57. Yon, L.; Duff, J.P.; Agren, E.O.; Erdelyi, K.; Ferroglio, E.; Godfroid, J.; Hars, J.; Hestvik, G.; Horton, D.; Kuiken, T.; et al. Recent Changes in Infectious Diseases in European Wildlife. J. Wildl. Dis. 2019, 55, 3–43. [Google Scholar] [CrossRef] [PubMed]
  58. Tischler, N.D.; Galeno, H.; Rosemblatt, M.; Valenzuela, P.D. Human and rodent humoral immune responses to Andes virus structural proteins. Virology 2005, 334, 319–326. [Google Scholar] [CrossRef] [Green Version]
Viruses 13 01188 g001 550
Figure 1. Composition profiles of hantavirus NPs. Amino acids are represented on the x-axis: alanine (A; dark grey bar), cysteine (C; yellow bar), aspartate (D; bright red bar), glutamate (E; bright red bar), phenylalanine (F; mid blue bar), glycine (G; light grey bar), histidine (H; pale blue bar), isoleucine (I; green bar), lysine (K; blue bar), leucine (L; green bar), methionine (M; yellow bar), asparagine (N; cyan bar), proline (P; flesh bar), glutamine (Q; cyan bar), arginine (R; blue bar), serine (S; orange bar), threonine (T; orange bar), valine (V; green bar), tryptophan (W; pink bar), tyrosine (Y; mid blue bar). Depleted and/or enriched amino acids were indicated by *. Statistical significance associated with a specific enrichment or depletion is estimated using the two-sample t-test between two sequences of binary indicator variables, one sequence for each of the samples. A particular enrichment or depletion is statistically significant when the p-value (the lowest value at which the null hypothesis that the same underlying Gaussian distribution generated for both samples can be rejected) is lower than or equal to a defined statistical significance value (p = 0.05).
Figure 1. Composition profiles of hantavirus NPs. Amino acids are represented on the x-axis: alanine (A; dark grey bar), cysteine (C; yellow bar), aspartate (D; bright red bar), glutamate (E; bright red bar), phenylalanine (F; mid blue bar), glycine (G; light grey bar), histidine (H; pale blue bar), isoleucine (I; green bar), lysine (K; blue bar), leucine (L; green bar), methionine (M; yellow bar), asparagine (N; cyan bar), proline (P; flesh bar), glutamine (Q; cyan bar), arginine (R; blue bar), serine (S; orange bar), threonine (T; orange bar), valine (V; green bar), tryptophan (W; pink bar), tyrosine (Y; mid blue bar). Depleted and/or enriched amino acids were indicated by *. Statistical significance associated with a specific enrichment or depletion is estimated using the two-sample t-test between two sequences of binary indicator variables, one sequence for each of the samples. A particular enrichment or depletion is statistically significant when the p-value (the lowest value at which the null hypothesis that the same underlying Gaussian distribution generated for both samples can be rejected) is lower than or equal to a defined statistical significance value (p = 0.05).
Viruses 13 01188 g001
Viruses 13 01188 g002 550
Figure 2. Analysis of hantavirus NPs variability. The variability analyses were evaluated comparing Laibin mobatvirus (YP_009507251), Xuan Son virus (AHA83416), Brno loanvirus (APU53639), Dakrong virus (AXQ03852), Lena River virus (QBH68029), Longquan loanvirus (YP_009664869), Nova mobatvirus (YP_009362031), Quezon mobatvirus (YP_009361846), Seoul orthohantavirus (ABM67564), Hantaan orthohantavirus (ALA20831), Andes orthohantavirus (AAK14323), Araraquara virus (ABQ45402), and Sin Nombre orthohantavirus (NP_941975) NPs with consensus sequence. The graph y-axis represents the Shannon entropy. The traced line indicates the threshold value of conservation, and values ≤ 1.0 indicate highly conserved positions.
Figure 2. Analysis of hantavirus NPs variability. The variability analyses were evaluated comparing Laibin mobatvirus (YP_009507251), Xuan Son virus (AHA83416), Brno loanvirus (APU53639), Dakrong virus (AXQ03852), Lena River virus (QBH68029), Longquan loanvirus (YP_009664869), Nova mobatvirus (YP_009362031), Quezon mobatvirus (YP_009361846), Seoul orthohantavirus (ABM67564), Hantaan orthohantavirus (ALA20831), Andes orthohantavirus (AAK14323), Araraquara virus (ABQ45402), and Sin Nombre orthohantavirus (NP_941975) NPs with consensus sequence. The graph y-axis represents the Shannon entropy. The traced line indicates the threshold value of conservation, and values ≤ 1.0 indicate highly conserved positions.
Viruses 13 01188 g002
Viruses 13 01188 g003 550
Figure 3. Mapping of predicted B-cell conformational epitopes is highlighted on the hantavirus NPs models. Mapping of: (a,c,e,g,i,k,m,o,q) secondary structural elements and (b,d,f,h,j,l,n,p,r) surface accessibility.
Figure 3. Mapping of predicted B-cell conformational epitopes is highlighted on the hantavirus NPs models. Mapping of: (a,c,e,g,i,k,m,o,q) secondary structural elements and (b,d,f,h,j,l,n,p,r) surface accessibility.
Viruses 13 01188 g003
Table
Table 1. General information of studied hantaviruses and their NPs.
Table 1. General information of studied hantaviruses and their NPs.
Taxonomic DataEcological DataNucleocapsid
Protein Data
HantavirusPhylogroupGenusReservoir
Order		Reservoir
Families		Reservoir SpeciesVirus DistributionLengthMolecular Weight (Da)
Andes virus (ANDV)VOrthohantavirusRodentiaCricetidaeOligoryzomys longicaudatusArgentina, Chile42848,039.69
Araraquara virus (ARAV)VOrthohantavirusRodentiaCricetidaeNecromys lasiurusBrazil42847,878.46
Sin Nombre virus (SNV)VOrthohantavirusRodentiaCricetidaePeromyscus maniculatusUSA, Canada, and Mexico42848,174.57
Seoul virus (SEOV)IVOrthohantavirusRodentiaMuridaeRattus rattus and R. norvegicusWorldwide42947,930.31
Hantaan virus (HTNV)IVOrthohantavirusRodentiaMuridaeApodemus agrariusChina, Korea, Japan, and Russia42948,152.61
Brno virus (BRNV)IILoanvirusChiropteraVespertilionidaeNyctalus noctulaCzech Republic42348,699.39
Longquan virus (LQUV)IILoanvirusChiropteraRhinolophidaeRhinolophus sinicus; R. affinis
and R. monoceros		China42348,193.69
Nova virus (NVAV)IIMobatvirusSoricomorphaTalpidaeTalpa europaeaBelgium, France, Hungary,
and Poland		42848,418.26
Laibin virus (LAIV)IIMobatvirusChiropteraEmballonuridaeTaphozous melanopogonChina, Myanmar42748,108.92
Quezon virus (QZNV)IIMobatvirusChiropteraPteropodidaeRousettus amplexicaudatusPhilippines42948,441.9
Lena River virus (LENV)IIMobatvirusSoricomorphaSoricidaeSorex caecutiensRussia44850,532.14
Xuan Son virus (XSV)IIMobatvirusChiropteraHipposideridaeHipposideros pomon;
H. cineraceus		Vietnam42748,136.93
Dakrong virus (DKGV)IIUnclassifiedChiropteraHipposideridaeAselliscus stoliczkanusVietnam42748,612.78
Robina virus (ROBV)IIOrthohantavirusChiropteraPteropodidaePteropus alectoAustralia42948,525.26
Table
Table 2. Geometry analysis of hantavirus NPs models.
Table 2. Geometry analysis of hantavirus NPs models.
HantavirusPoor
Rotamers		Favored
Rotamers		Ramachandran
Outliers		Ramachandran
Favored		Z-Score
Brno037714121.15 ± 0.39
0.00%99.21%0.24%97.86%
Longquan036914140.66 ± 0.40
0.00%99.73%0.24%98.34%
Nova036504180.88 ± 0.38
0.00%99.46%0.00%98.12%
Laibin136404130.59 ± 0.38
0.27%99.18%0.00%97.18%
Quezon037504201.04 ± 0.38
0.00%99.47%0.00%98.36%
Lena River03782438−0.08 ± 0.37
0.00% 98.69%0.45%98.21%
Xuan Son036704180.20 ± 0.37
0.00%99.73%0.00%98.35%
Dakrong037004180.43 ± 0.38
0.00%98.93%0.00%98.35%
Robina036924200.79 ± 0.39
0.00%100%0.47%98.36%
Goals<0.3%>98%<0.05%>98%<2
Table
Table 3. Analysis of epitopes’ conservation of hantavirus phylogroup II NPs.
Table 3. Analysis of epitopes’ conservation of hantavirus phylogroup II NPs.
VirusEpitopeLengthConservation among OrthohantavirusConservation among Phylogroup II
High Similar NPs (≥70%)
(%—N/5)		Mean Conservation
% (Min–Max)		High Similar NPs (≥70%)
(%—N/9)		Mean Conservation
% (Min–Max)
BRNVB1350.00 (0/5)21.1 (20–22.9)11.11 (1/9)35.9 (20–100)
B2280.00 (0/5)40 (35.7–42.9)11.11 (1/9)50.4 (21.4–100)
B31030.00 (0/5)22.1 (20.4–24.3)11.11 (1/9)35.5 (21.4–100)
B47100.00 (5/5)80 (71.4–85.7)88.89 (8/9)92.1 (57.1–100)
LQUVLQ1250.00 (0/5)54.4 (48–60)11.11 (1/9)52.4 (44–100)
LQ2370.00 (0/5)22.7 (18.9−29.7)11.11 (1/9)36.6 (18.9–100)
LQ3880.00 (0/5)44.5 (42–46.6)11.11 (1/9)55.6 (40.9–100)
NVAVN1440.00 (0/5)35.5 (22.7–43.2)11.11 (1/9)31.8 (18.2–100)
N2530.00 (0/5)38.9 (32.1–45.3)11.11 (1/9)42.6 (7.6–100)
N3570.00 (0/5)47 (43.9–49.1)11.11 (1/9)57.3 (31.6–100)
N490.00 (0/5)53.3 (44.4–55.6)11.11 (1/9)66.7 (55.6–100)
LAIVLB1720.00 (0/5)27.5 (26.4–29.2)22.22 (2/9)53.1 (33.3–100)
LB2650.00 (0/5)39.7 (35.4–44.6)33.33 (3/9)60.9 (40–100)
LB38100.00 (5/5)80 (75–87.5)88.89 (8/9)93.1 (62.5–100)
QZNVQ1950.00 (0/5)14.1 (10.5–18.9)22.22 (2/9)29.5 (8.4–100)
Q2320.00 (0/5)55 (50–56.3)22.22 (2/9)62.5 (37.5–100)
Q3640.00 (2/5)73.3 (66.7–83.3)88.89 (8/9)90.7 (50–100)
LENVLR1100.00 (0/5)50 (50–50)11.11 (1/9)61.9 (57.1–100)
LR2850.00 (0/5)10.4 (9.4–11.8)11.11 (1/9)20.7 (7.1–100)
LR3360.00 (0/5)25.6 (22.2–30.6)11.11 (1/9)30.6 (19.4–100)
LR4480.00 (4/5)80 (50–100)100.00 (9/9)77.8 (75–100)
LR5120.00 (0/5)41.7 (41.7–41.7)11.11 (1/9)42.6 (33.3–100)
XSVXS1990.00 (0/5)20.8 (18.9–24.2)22.22 (2/9)48.4 (19.2–100)
XS21540.00 (2/5)65.3 (60–73.3)77.78 (7/9)75.6 (46.7–100)
XS311100.00 (5/5)80 (72.7–81.82)88.89 (8/9)79.8 (54.6–100)
XS440.00 (0/5)50 (50–50)22.22 (2/9)58.3 (50–100)
DKGVD150.00 (0/5)52 (40–60)55.56 (5/9)77.8 (60–100)
D21180.00 (0/5)28.6 (26.3–30.5)33.33 (3/9)49.4 (17.8–100)
D3590.00 (0/5)16.3 (15.3–18.6)11.11 (1/9)40.1 (15.3–100)
D4640.00 (2/5)73.3 (66.7–83.3)88.89 (8/9)90.7 (50–100)
ROBVR11020.00 (0/5)15.4 (11.1–20.2)22.22 (2/9)31.6 (13.1–100)
R2630.00 (0/5)45.1 (42.9–47.6)22.22 (2/9)59.4 (31.8–100)
R3640.00 (2/5)73.4 (66.7–83.3)88.89 (8/9)90.7 (50–100)
Table
Table 4. Conservation degree of hantavirus phylogroup II NPs conformational epitopes.
Table 4. Conservation degree of hantavirus phylogroup II NPs conformational epitopes.
VirusEpitopeLengthConservation Degree (%)
Rodent-BorneBat-BorneShrew-Borne
ANDVARAVSNVSEOVHTNVBRNVLQUVLAIVQZNVXSVDKGVROBVLENVNVAV
BrnoB13520202320231005420262329262620
B22843434336361005743433954572139
B310320212024241004322252127243026
B478686718671100100100100861001008657
LongquanLQ12548485660605610044524448602444
LQ23719191927302210041244335241922
LQ38847474243446910047504148515242
NovaN14443434325233227182318202720100
N2534543423232813495147424528100
N35749494447464447566060615632100
N4956565644565667676756676756100
LaibinLB17228262628292128100577869603233
LB26537383543454045100578282622952
LB388888757575100100100100881001008863
QuezonQ195111112191998151001414712213
Q23250505656635356531006659813856
Q368383676767100100100100831001008350
Lena RiverLR11050505050505757575757575710057
LR2859911121171111139121310011
LR33622222528312522192219222210022
LR447575100501007575757575757510075
LR51242424242423342333333334210033
Xuan SonXS19920191822241929784010070432233
XS21560606073736773808010087734773
XS1182827382827382918210082827355
XS4450505050505075505010050505050
DakrongD156060604040606010080100100606080
D211829282630311820755170100482537
D35919171515152425463447100471522
D468383676767100100100100831001008350
RobinaR110211121320201313177317161002213
R26343444844465148607359631003249
R368383676767100100100100831001008350
The epitopes’ degree of conservation among studied hantaviruses were indicated in the table (%) and classified as low conservancy (light-blue cells; conservancy < 40%), intermediate conservancy (yellow cells, ≥40% conservancy < 70%), and high conservancy (red cells, conservancy ≥ 70%). The columns represent hantavirus species: Andes orthohantvirus (ANDV); Araraquara virus (ARAV); Sin Nombre orthohantavirus (SNV); Seoul orthohantavirus (SEOV); Hantaan orthohantavirus (HTNV); Brno loanvirus (BRNV); Longquan loanvirus (LQUV); Nova mobatvirus (NVAV); Laibin mobatvirus (LAIV); Quezon mobatvirus (QZNV); Lena River mobatvirus (LENV); Xuan Son virus (XSV); Dakrong virus (DKGV); Robina orthohantavirus (ROBV).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

de Oliveira, R.C.; Fernandes, J.; de Sampaio Lemos, E.R.; de Paiva Conte, F.; Rodrigues-da-Silva, R.N. The Serological Cross-Detection of Bat-Borne Hantaviruses: A Valid Strategy or Taking Chances? Viruses 2021, 13, 1188. https://0-doi-org.brum.beds.ac.uk/10.3390/v13071188

AMA Style

de Oliveira RC, Fernandes J, de Sampaio Lemos ER, de Paiva Conte F, Rodrigues-da-Silva RN. The Serological Cross-Detection of Bat-Borne Hantaviruses: A Valid Strategy or Taking Chances? Viruses. 2021; 13(7):1188. https://0-doi-org.brum.beds.ac.uk/10.3390/v13071188

Chicago/Turabian Style

de Oliveira, Renata Carvalho, Jorlan Fernandes, Elba Regina de Sampaio Lemos, Fernando de Paiva Conte, and Rodrigo Nunes Rodrigues-da-Silva. 2021. "The Serological Cross-Detection of Bat-Borne Hantaviruses: A Valid Strategy or Taking Chances?" Viruses 13, no. 7: 1188. https://0-doi-org.brum.beds.ac.uk/10.3390/v13071188

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop