Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Fine mapping of the antigenic epitopes of the Gc protein of Guertu virus

  • Meilipaiti Yusufu ,

    Contributed equally to this work with: Meilipaiti Yusufu, Ayipairi Abula

    Roles Conceptualization, Data curation, Methodology, Software, Validation, Writing – original draft

    Affiliation Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China

  • Ayipairi Abula ,

    Contributed equally to this work with: Meilipaiti Yusufu, Ayipairi Abula

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Writing – original draft

    Affiliation Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China

  • Boyong Jiang,

    Roles Data curation, Formal analysis, Methodology, Software, Visualization

    Affiliation Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China

  • Jiayinaguli Zhumabai,

    Roles Formal analysis, Investigation, Supervision, Visualization

    Affiliation Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China

  • Fei Deng,

    Roles Investigation, Project administration, Supervision, Writing – review & editing

    Affiliation State Key Laboratory of Virology, Wuhan Institute of Virology, Chinese Academy of Sciences, Wuhan, China

  • Yijie Li,

    Roles Formal analysis, Investigation, Software, Supervision

    Affiliation Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China

  • Yujiang Zhang ,

    Roles Investigation, Software, Supervision, Writing – review & editing

    sr_sun2005@163.com (SS); dingjuntao2004@126.com (JD); xjsyzhang@163.com (YZ)

    Affiliation Center for Disease Control and Prevention of Xinjiang Uygur Autonomous Region, Urumqi, China

  • Juntao Ding ,

    Roles Conceptualization, Investigation, Supervision, Visualization, Writing – original draft

    sr_sun2005@163.com (SS); dingjuntao2004@126.com (JD); xjsyzhang@163.com (YZ)

    Affiliation Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China

  • Surong Sun

    Roles Investigation, Software, Supervision, Writing – review & editing

    sr_sun2005@163.com (SS); dingjuntao2004@126.com (JD); xjsyzhang@163.com (YZ)

    Affiliation Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China

Abstract

Guertu virus (GTV), a newly discovered member of the genus Banyangvirus in the family Phenuiviridae, poses a potential health threat to humans and animals. The viral glycoprotein (GP) binds to host cell receptors to induce a neutralizing immune response in the host. Therefore, identification of the B-cell epitopes (BCEs) in the immunodominant region of the GTV Gc protein is important for the elucidation of the virus–host cell interactions and the development of GTV epitope assays and vaccines. In this study, an improved overlapping biosynthetic peptide method and rabbit anti-GTV Gc polyclonal antibodies were used for fine mapping of the minimal motifs of linear BCEs of the GTV Gc protein. Thirteen BCE motifs were identified from eleven positive 16mer-peptides, namely EGc1 (19KVCATTGRA27), EGc2 (58KKINLKCKK66), EGc3 (68SSYYVPDA75), EGc4 (75ARSRCTSVRR84), EGc5 (79CTSVRRCRWA88), EGc6 (90DCQSGCPS97), EGc7 (96PSHFTSNS103), EGc8 (115AGLGFSG121), EGc9 (148ENPHGVI154), EGc10 (179KVFHPMS185), EGc11 (230QAGMGVVG237), EGc12 (303RSHDSQGKIS312), and EGc13 (430DIPRFV435). Of these, 7 could be recognized by GTV IgG-positive sheep sera. Three-dimensional structural analysis revealed that all 13 BCEs were present on the surface of the Gc protein. Sequence alignment of the 13 BCEs against homologous proteins from 10 closely related strains of severe fever with thrombocytopenia syndrome virus from different geographical regions revealed that the amino acid sequences of EGc4, EGc5, EGc8, EGc11, and EGc12 were highly conserved, with 100% similarity. The remaining 8 epitopes (EGc1, EGc2, EGc3, EGc6, EGc7, EGc9, EGc10, and EGc13) showed high sequence similarity in the range of 71.43%–87.50%. These 13 BCEs of the GTV Gc protein provide a molecular foundation for future studies of the immunological properties of GTV glycoproteins and the development of GTV multi-epitope assays and vaccines.

Introduction

Ticks, the species of which are either generalist or specialist parasites, are responsible for transmission of a variety of pathogens to animals and humans [1]. The incidence of tick-borne viral diseases has increased over the past few decades, thereby compromising the health of humans significantly. In recent years, tick-borne severe fever with thrombocytopenia syndrome virus (SFTSV) and heartland virus (HRTV) infections have been shown to be associated with serious human diseases and fatalities in the USA and East Asian countries [24]. Studies on SFTSV and HRTV indicate that novel tick-borne viruses may pose a significant threat to public health. Therefore, the discovery and investigation of novel tick-borne viruses are of great importance for ensuring the control and prevention of potential disease outbreaks.

Guertu virus (GTV), which is a member of the newly defined genus Banyangvirus belonging to the family Phenuiviridae of Order Bunyavirales [5], was first isolated in 2014 from Dermacentor nuttalli, a tick found in the Guertu mountain region of Wusu, Xinjiang, China [6]. Researchers reported the morphological characteristics of GTV for the first time, and phylogenetic analysis showed that it was closely related to SFTSV [710]. The genome of GTV consists of the following three single-stranded RNA segments: an S segment that contains open reading frames that encode a nonstructural (NS) protein and a nucleoprotein (NP); an L segment that encodes an RNA-dependent RNA polymerase; and M segment that encodes a glycoprotein precursor that is cleaved into two mature envelope proteins (Gn and Gc) during co-translational modification [6]. Studies have shown that viral expression of the Gn protein is directly correlated with the Gc protein-facilitated binding of the viral particles to the host cell surface and that the Gc protein is efficiently translated from the Gn/Gc-encoding mRNA, thus driving the fusion between the virus and the host cell membrane [11]. At present, presumable N-glycosylation sites is located on SFTSV Gc glycoproteins [10, 11]. DC-SIGN has been identified as the receptor factor on the host cell membrane [10]. Some studies have also shown that HS and NMMHC-IIA may also be the entry factors related to cell entry of SFTSV [12]. Additionally, Gn and Gc facilitate viral entry into the target cells and are the primary targets of neutralizing antibodies [13].

B-cell epitopes (BCEs), which are recognized by B-cell receptors, induce cellular and humoral immune responses in the host [14]. Therefore, the identification of viral BCEs would deepen our understanding of the host immune response, provide candidate epitope peptides for the development of epitope vaccines or the establishment of disease diagnostic methods, and help to reveal the mechanisms of therapeutic antibodies [15, 16]. Specifically, the identification of BCEs of the conserved domain or immunodominant region of the Gc protein of GTV is crucial for the elucidation of the virus–host cell interactions and development of GTV assays and vaccines [17]. However, up to date, there are no reports on epitope mapping of GTV Gc segment and their fine localization.

The currently available epitope mapping methods include recombinant DNA techniques, biosynthetic peptide (BSP) methods, chemically synthesized peptides or peptide microarrays, and phage-display libraries [1822]. The BSP method developed by Xu et al. in 2009 was used for the identification of linear BCEs and minimal epitope motifs, and they successfully mapped linear of BCEs of three structural proteins of human papillomavirus type 58, using anti-recombinant protein pAbs [23, 24]. In our previous study, this method was used to identify the linear BCEs of structural proteins from a strain of Crimean-Congo hemorrhagic fever virus that was isolated in Xinjiang, where its feasibility for epitope identification was demonstrated [2527]. Additionally, Zhang et al. [17] used pAbs and the BSP method to identify nine BCEs on the Gn protein of GTV. However, the identification of BCEs on the Gc protein of GTV has not been reported thus far.

Therefore, this study was performed to identify the minimal BCEs on the Gc protein of GTV, using biosynthetic peptide (BSP) methods [14], and rabbit anti-GTV-Gc pAbs. Subsequently, 13 minimal epitope motifs were obtained and their antigenicity was confirmed using positive sera from sheep naturally infected with GTV. Python 2.7 was used to analyze the three-dimensional (3D) conformation of the Gc protein, thus revealing that the identified BCEs were all located on the protein surface. This study provides fundamental data that can be used in future studies aimed at elucidating the immunological properties of the Gc protein of GTV and for developing diagnostic reagents and vaccines for infection prevention and control.

Materials and methods

Ethical statement

This study was approved by the Animal Ethics Committee of the Key Laboratory of Bioresources and Genetic Engineering, Xinjiang University (Approval number: BRGE-AE001). The animal serum samples were collected using a random sampling method that did not involve sacrifice of the animals, and serum collection was performed according to this approved animal protocol [14].

Antibodies

The rabbit pAbs raised against GTV Gc1 (aa 1–291) and Gc2 (aa 252–549) were prepared by the Wuhan Institute of Virology, Chinese Academy of Sciences. His-tag mouse monoclonal antibody, goat anti-mouse horseradish peroxidase-conjugated IgG (IgG-HRP), and goat anti-rabbit IgG-HRP were purchased from Beijing Solarbio Science & Technology Co., Ltd. The serum samples from sheep infected or not infected with GTV were previously identified by indirect immunofluorescence assay and reverse transcription-polymerase chain reaction (RT-PCR) assay [6].

Construction of the truncated Gc fragments for recombinant plasmid expression

The prokaryotic expression plasmid pET-32a (+) containing the full-length Gc gene fragment of GTV strain DXM (GenBank accession no. KT328592.1) was stored at our laboratory. The prokaryotic expression vector pXXGST-3 was provided by Professor Xu Wanxiang of the Shanghai Institute for Biomedical Pharmaceutical Technologies, China. The Gc protein of GTV was divided into two truncated fragments Gc1 and Gc2 (Fig 1). These fragments were then PCR-amplified using the previously constructed pET-32a-GTV-Gc plasmid as a template [28] for the design of the following primers: GcI F: CGGGATCCATGAAACTGGTTCGCTTAACC; GcI R: AATGCGGCCGCTTATGATTTCAGTTCTTCAATTTCAC; GcII F: CGGGATCCATGAGTGGTGTTCCGACCCAG; and GcII R: AATGCGGCCGCTTAGGTCAGCAGCACGCCT.

thumbnail
Fig 1. Schematic of epitope mapping strategy.

(A) The black band indicates the full-length sequence of GTV-Gc, and the gray bands indicate two of truncated GcⅠ and GcⅡ segments. (B) Schematic of epitope mapping strategy involves 57 overlapping 16mer-peptides spanning Gc sequences.

https://doi.org/10.1371/journal.pone.0271878.g001

Strategy for designing GTV Gc overlapping peptides

To map the linear BCEs of the GTV Gc protein, we designed 16mer-peptides with eight overlapping amino acid residues for the first round of antigenic peptide mapping (Fig 1). Additionally, a series of 8mer-peptides with seven overlapping amino acid residues that spanned every positive 16mer-peptides sequence were also designed for the second round of fine mapping of the BCE motif [29]. The GTV Gc protein was designed into fifty-seven 16mer-peptides using modified overlapping peptide biosynthesis. These 16mer-peptides were further shortened into a total of seventy-nine 8mer-peptides overlapping 7mer-peptides. The BSP method had the following features: (1) By selecting the truncated GST188 protein as the fusion tag for the recombinant expression of short peptides, the fragment size of the expressed fusion peptide was maintained within 21.5–22.5 kDa, which allowed its easy identification by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and subsequent western blot analysis. This also facilitated the easy distinction of the short peptides from the proteins expressed by the host (Escherichia coli); (2) The antigenic proteins were first divided into 16mer-peptides with eight overlapping amino acid residues, following which basic and complete screening was performed for those containing BCEs; (3) Further truncation of the positive 16mer-peptides into 8mer-peptides with seven overlapping amino acid residues allowed the precise determination of the minimal epitope motif for each linear BCE [14].

Construction and prokaryotic expression of the GTV Gc overlapping peptides

A suspension of bacterial cells harboring the empty prokaryotic expression plasmid pXXGST-3 was added to 10 mL of liquid medium containing 0.1% ampicillin and incubated overnight. The plasmid was then extracted using the Plasmid Mini-Prep Kit. The DNA fragments intended for incorporation into the plasmid were first digested with BamH I and Sal I at a constant temperature of 37°C for 4–5 h. According to the biosynthesis report of the biotechnology company, we considered the final concentration of each synthesized DNA fragment as 5 μM per tube. These fragments were mixed and heated to 95°C for 5 min in a constant temperature bath and then annealed by natural cooling to ambient temperature. The annealed oligonucleotides were mixed with the digested pXXGST-3 plasmid and then ligated at 16°C for 12–15 h. The vector was transfected into E. coli BL21(DE3) competent cells and cultured overnight, following which single clones were added to ampicillin-containing Luria broth and incubated overnight at 37°C under shaking conditions. On the next day, 1% of the overnight culture was inoculated into new Luria broth and cultured at 30°C under shaking conditions until the OD600 reached 0.6–0.8. Heat induction at 42°C was then conducted for another 4 h [29].

SDS-PAGE and western blot analysis of the GTV Gc overlapping peptides

The bacterial culture was centrifuged at 12000 rpm for 1 min, and the supernatant was discarded. After resuspending the cells in 160 μL of 1×phosphate-buffered saline, 40 μL of 5×loading buffer was added and the mixture was boiled at 100°C for 15 min and then centrifuged at 12000 rpm for 15 min. The protein samples were then separated by 15% SDS-PAGE. The plasmids expressing the correct protein size were send to Shanghai Sangon for sequencing. Western blot analysis was performed on the recombinant proteins, using rabbit anti-Gc pAbs to analyze their antigenicity. The proteins from strains with the correct sequencing results were separated by 15% SDS-PAGE. Polyvinylidene difluoride (PVDF) membranes were purchased from Whatman International Limited. The separation gel was then removed with a transfer clamp and placed in an electrophoresis tank for electrotransfer of the proteins onto a PVDF membrane at 90 V for 90 min. Thereafter, the PVDF membrane was immersed in blocking buffer at 37°C for 2 h, washed with Tris-buffered saline-Tween 20 (TBS-T), and then incubated with rabbit serum (1:300) at 4°C overnight. Subsequently, the PVDF membrane was first rinsed with the washing buffer for 1 h and then incubated with HRP-labeled secondary antibody (1:2000) at 37°C for 1 h. Finally, after rinsing the membrane with the washing buffer for 1 h, the ECL (Wuhan Boster Biological Technology Ltd.) chromogenic reagent was added and the membrane was imaged using the LAS-4000 chemiluminescence imaging system.

Analyses of the similarities of the conformational antigenic epitopes and their localization on Gc by three-dimensional modeling

To analyze the similarity of each BCE with homologous proteins, Gc sequences of SFTSV strains from different countries and genetic lineages were downloaded from the GenBank database on the basis of the phylogenetic tree of SFTSV strains [9]. The BCE locations were experimentally identified from the 3D structure of the Gc protein using the PyMOL software (https://pymol.org/2/). Predictions of the secondary structure were conducted according to the Garnier-Robson [30] and Chou-Fasman [31] methods. The hydrophilicity, flexibility, surface accessibility, and antigenicity indices were also analyzed and predicted.

Results

Induced expression of the overlapping 16mer-peptides and western blot analysis of their antigenicity

The SDS-PAGE results showed that the 57 overlapping 16mer-peptides spanning the GTV Gc1-549 fragment were successfully expressed under the 42°C heat induction conditions (Fig 2A and 2C). The corresponding western blot analysis of their antigenicity showed that P3 (Gc17-32), P8 (Gc57-72), P9 (Gc65-80), P10 (Gc73-89), P12 (Gc89-104), P15 (Gc113-128), P19 (Gc145-160), P23 (Gc177-192), P29 (Gc225-240), P38 (Gc297-312), and P54 (Gc425-440) could bind specifically to the rabbit pAbs (Fig 2B and 2D) and were therefore potential linear antigenic epitope motifs.

thumbnail
Fig 2. SDS-PAGE and Western blot analysis of expressed 16mer-peptides.

(A and C) SDS-PAGE analysis of expressed 16mer-peptides. The numbers of P1-P57 indicate each 16mer-peptides in cell total proteins. The cell proteins of each recombinant clone were resolved by 12% SDS-PAGE and stained with Coomassie brilliant blue. M, the protein molecular marker; PC, Positive control of expressed 16mer peptide in GTV-NP; NC, Negative control of GST188 carrier protein expressed by pXXGST-2. (B and D) Western blot analysis for mapping reactive 16mer-peptides in P1-P57. The rabbit antiserum (1:300 dilution) against GTV-Gc was used in Western blot analysis. The reactive bands in Western blot analysis were visualized by enhanced chemiluminescence.

https://doi.org/10.1371/journal.pone.0271878.g002

Induced expression and western blot analysis of the overlapping 8mer-peptides

To identify the minimal antigenic epitopes that could bind specifically to the rabbit pAbs, the eleven positive 16mer peptide fragments were further designed as 8mer-peptides with seven overlapping amino acids. The SDS-PAGE and sequencing results showed that all overlapping 8mer-peptides with a GST188 tag were correctly expressed in E. coli. Moreover, the western blot analysis confirmed that P9-3 (Gc68-75) in the P9 epitope was specifically recognized by the rabbit pAbs, indicating that the minimal motif of the P9 epitope was “SSYYVPDA” (named EGc3). Similarly, in the P12, P15, P19, P23, and P54 epitopes, the minimal motifs recognized by the rabbit pAbs were found to be P12-1 (Gc90-97), P12-7 (Gc96-103), P15-1 (Gc114-121), P15-2 (Gc115-122), P19-2 (Gc147-154), P19-3 (Gc148-155), P23-1 (Gc178-185), P23-2 (Gc179-186), P54-3 (Gc428-435), P54-4 (Gc429-436), and P54-5 (Gc430-437). Interestingly, none of the overlapping 8mer-peptides generated from the other five positive 16mer-peptides (i.e., P3, P8, P10, P29, and P38) showed a positive reaction, suggesting that their minimal epitope motifs might contain longer antigenic fragments. Finally, as shown in Fig 3, seven BCEs were identified, namely EGc3 (68SSYYVPDA75), EGc6 (90DCQSGCPS97), EGc7 (96PSHFTSNS103), EGc8 (115AGLGFSG121), EGc9 (148ENPHGVI154), EGc10 (179KVFHPMS185), and EGc13 (430DIPRFV435).

thumbnail
Fig 3. SDS-PAGE and Western blot analysis of expressed 8mer-peptides.

(A, C, E and G) SDS-PAGE analysis of expressed 8mer-peptides. It indicates each 8mer peptide of P3, P8, P9, P10, P12, P15, P19, P23, P29, P38 and P54. The cell proteins of each recombinant clone were resolved by 12% SDS-PAGE and stained with Coomassie brilliant blue. M, the protein molecular marker. NC, Negative control of GST188 protein. PC, Positive control of mapped reactive P3. (B, D, F and H) Western blot analysis for mapping fine epitopes in each reactive 8-mer-peptides. The rabbit antiserum against GTV-Gc (1:300 dilution) was used in Western blot analysis. The reactive bands in Western blot analysis were visualized by enhanced chemiluminescence.

https://doi.org/10.1371/journal.pone.0271878.g003

Induced expression and western blot analysis of the overlapping 10mer-peptides

To identify the minimal antigenic epitope motifs of the positive 16mer-peptides P3, P8, P10, P29, and P38 that could bind specifically to rabbit pAbs, we designed them as 10mer-peptides with seven overlapping amino acid residues for validation. The SDS-PAGE and sequencing results showed that all overlapping 10mer-peptides with a GST188 tag were correctly expressed in E. coli, and western blot analysis showed that P3-9 (Gc18-27), P3-10 (Gc19-28), P8-8 (Gc57-66), P8-9 (Gc58-67), P10-14 (Gc76-85), P10-17 (Gc79-88), P29-11 (Gc228-237), P29-12 (Gc229-238), P29-13 (Gc230-239), and P38-16 (Gc303-312) were specifically recognized by the rabbit pAbs. The minimal motifs in these epitopes were EGc1 (19KVCATTGRA27), EGc2 (58KKINLKCKK66), EGc4 (75ARSRCTSVRR84), EGc5 (79CTSVRRCRWA88), EGc11 (230QAGMGVVG237), and EGc12 (303RSHDSQGKIS312) (Fig 4). Thus, a total of 13 minimal antigenic epitope motifs were identified in GTV Gc1-549. The overlapping 16/10/8mer peptide sequences and their corresponding locations on the Gc protein are shown in Fig 5.

thumbnail
Fig 4. SDS-PAGE and Western blot analysis of expressed 10mer-peptides.

(A and C) SDS-PAGE analysis of expressed 10mer-peptides for P3, P8, P10, P29 and P38. The cell proteins of each recombinant clone were resolved by 12% SDS-PAGE and stained with Coomassie brilliant blue. M, the protein molecular marker; NC, Negative control of GST188 carrier protein. PC, Positive control of P3. (B and D) Western blot analysis for mapping fine epitopes in each reactive 10mer peptide. The rabbit antiserum against GTV-Gc (1:300 dilution) was used in Western blot analysis. The reactive bands in Western blot analysis were visualized by enhanced chemiluminescence.

https://doi.org/10.1371/journal.pone.0271878.g004

thumbnail
Fig 5. The bio-synthetic overlapping peptides derived from IgG-reactive peptides of Gc1-459.

The yellow highlighting represents the common sequences among immunodominant peptides that react with pAbs according to western blot analysis.

https://doi.org/10.1371/journal.pone.0271878.g005

Validation of the antigenicity of the peptides using GTV-positive sheep sera

To investigate whether BCE screening with a rabbit anti-GTV Gc polyclonal antiserum could detect the IgG antibodies present naturally in livestock, the 13 short peptides of BCEs identified in this experiment were analyzed by western blot analysis using the positive serum of a sheep naturally infected with GTV and the negative serum of an uninfected sheep. The results showed that all 13 short peptides containing BCEs were correctly expressed (Fig 6A). Among the 13 minimal BCE motifs identified, only EGc2, EGc3, EGc4, EGc9, EGc10, EGc11, and EGc12 were specifically recognized by the positive sheep serum (Fig 6B). None of the BCEs identified could be recognized by the negative sheep serum (Fig 6C). The Western blot results showed that there was a good consistency between the immune response of rabbit and sheep to GTV GC protein.

thumbnail
Fig 6. Western blot analysis of mapped BCEs using sheep sera.

(A) SDS-PAGE analysis of expressed BCEs. (B) Using a positive serum from a sheep confirmed GTV-infection. (C) Using a serum from healthy sheep with no history of GTV infection as a negative control. PC, Positive control of GTV-Gc1.

https://doi.org/10.1371/journal.pone.0271878.g006

Sequence conservation analysis and three-dimensional modeling

Using the Gendoc software, the amino acid sequences of the 13 BCEs were aligned with those of homologous proteins from 10 genetically closely related SFTSV strains from different regions (Japan, Accession no. AB817987.1; Japan, Accession no. AB985297.1; China, Accession no. HQ141605.1; China, Accession no. HQ830167.1; China, Accession no. HQ830167.1; China, Accession no. JQ670931.1; China, Accession no. KC505127.1; South Korea, Accession no. KF358692.1; China, Accession no. KP202164.1; China, Accession no. KX302598.1; and South Korea, Accession no. KY789438.1). The results revealed that EGc4 (75ARSRCTSVRR84), EGc5 (79CTSVRRCRWA88), EGc8 (115AGLGFSG121), EGc11 (230QAGMGVVG237), and EGc12 (303RSHDSQGKIS312) were identical among the representative SFTSV strains and showed 100% conserved sequences (Fig 7). Thus, they can be used as candidates for broad-spectrum multi-epitope vaccine design.

thumbnail
Fig 7. Sequence comparison between GTV-Gc1-549 and 10 SFTSV strains.

The GenBank codes and sources are shown at left and the sequence analysis was based on the ClustalW program. The nine of fine BCEs and APs recognized by pAbs are highlighted, and the variable aa residue within the BCE motif are highlighted in red. Dots (.) indicate identical aa residue in ten SFTSV strains (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article).

https://doi.org/10.1371/journal.pone.0271878.g007

As shown in Fig 7, there was variability between the other epitopes and the SFTSV sequences. The sequence similarities of EGc1, EGc2, EGc3, EGc6, EGc7, EGc9, EGc10, and EGc13 were 87.50%, 87.50%, 87.50%, 87.50%, 87.50%, 85.71%, 71.43%, and 71.43%, respectively. In EGc1 (19KVCATTGRA27), “V” was changed to “E” at position Gc20; in EGc2 (58KKINLKCKK66), “V” was changed to “R” at position Gc59; in EGc3 (68SSYYVPDA75), “Y” was changed to “F” at position Gc71; in EGc6 (90DCQSGCPS97) and EGc7 (96PSHFTSNS103), “S” was changed to “P” at position Gc97; in EGc9 (148ENPHGVI154), “V” was changed to “I” at position Gc153; in EGc10 (179KVFHPMS185), “K” was changed to “V” and “R” was changed to “T” at positions Gc179 and Gc180; and in EGc13 (430DIPRFV435), “I” was changed to “R” and “V” was changed to “K” at positions Gc431 and Gc433.

The 3D structure of the Gc protein of GTV strain DXM was modeled using the Python 2.7 software, where different color markers were used to locate the BCEs in the model. As shown in Fig 8B, all the identified antigenic epitopes were located on the surface of the Gc protein. Moreover, the 13 BCEs were located at positions that were easily accessible for antibody chimerization and thus could facilitate the binding of specific antibodies.

thumbnail
Fig 8. Prediction of Gc secondary structure and 3D localization of each BCEs.

(A) Epitope prediction for GTV-Gc using DNAStar-Protean software. The secondary structure, flexibility plot, hydrophilicity, surface probability, and antigenicity index for GTV-Gc were taken into consideration. (B) Location distribution on 3D structure of mapped BCEs and APs on molecular surface were shown in different colors. EGc1 (red), EGc2 (salmon), EGc3 (blue), EGc4 (slate), EGc5 (yellow), EGc6 (wheat), EGc7 (hot pink), EGc8 (light magenta), EGc9 (brown),EGc10 (firebrick), EGc11 (gray), EGc12 (orange), EGc13 (white). The figures were generated using the PyMOL™ molecular graphics system (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article).

https://doi.org/10.1371/journal.pone.0271878.g008

The secondary structure of Gc1-549 was predicted using the DNASTAR Protean software (DNASTAR Inc., Madison, WI, USA), which showed that the identified epitopes formed part of the β-sheet and α-helix regions (Fig 8B). Moreover, the identified epitopes were predicted to be located on the surface of the GTV Gc protein with a high antigenic index and hydrophilicity, suggesting that they might be important BCEs of this glycoprotein (Fig 8A). The 3D modeling and protein secondary structure analysis indicated that the predictions were consistent with the identified epitopes and that the identified BCEs were all located in the hydrophilic region.

Discussion

Bunyavirales are the largest order of RNA viruses with the ability to infect a wide range of hosts, including humans, arthropods, and plants [32]. Several emerging bunyaviruses, such as the Rift Valley fever virus [33], Crimean-Congo hemorrhagic fever virus [34], HRTV [7], and SFTSV [35], have been shown to pose a serious threat to human health. SFTSV can be transmitted from ticks to humans, with human-to-human transmission occurring in rare instances [3638], and can lead to the manifestation of diseases characterized by fever, gastrointestinal symptoms, leukopenia, and thrombocytopenia. SFTSV and its related viruses are important emerging pathogens for which there are currently no specific antiviral drugs or vaccines [32]. As a novel bunyavirus, GTV belongs to the genus Banyangvirus together with HRTV and SFTSV [5]. Shen et al. [6]. performed in vitro and in vivo experiments to characterize the infection properties of GTV in human- and animal-derived cells. They found that infection of C57BL/6 mice with GTV could cause viremia or organopathological changes, indicating that this virus was a potential pathogen capable of infecting humans and animals. Serological testing of local residents in the Guertu area of Wusu City, Xinjiang revealed that 19.8% were positive for anti-GTV antibodies, and three of them had neutralizing antibody activity against the virus, which further suggested the potential risk of this pathogen to human health.

Studies have shown that viral glycoproteins are the targets of neutralizing antibodies and that the envelope glycoproteins are responsible for the binding of the virus particles to cellular receptors and subsequent viral fusion [39]. In viral infections, almost all neutralizing polyclonal or monoclonal antibodies are specific to the Gn and Gc glycoproteins [40, 41]. Therefore, the identification of BCE domains in GTV glycoproteins will provide new targets for the design and development of vaccines and diagnostic reagents for this virus. Although several studies on the identification of BCEs of the GTV Gn protein have been reported [17], there are no published studies on epitope mapping of the GTV Gc protein or of its IgG motif recognition sites.

In this study, an improved BSP method was used to map linear BCEs spanning the full length of the GTV Gc protein. Using the overlapping BSP method, the Gc protein was divided into 57 overlapping 16mer-peptides. Of these, 11 positive 16mer-peptides were selected through western blot analysis using rabbit anti-GTV-Gc pAbs and were further divided into 79 overlapping 8mer-peptides and 35 overlapping 10mer-peptides. Using immunoblot analysis, the following 13 antigenic epitope motifs were finally identified: EGc1 (19KVCATTGRA27), EGc2 (58KKINLKCKK66), EGc3 (68SSYYVPDA75), EGc4 (75ARSRCTSVRR84), EGc5 (79CTSVRRCRWA88), EGc6 (90DCQSGCPS97), EGc7 (96PSHFTSNS103), EGc8 (115AGLGFSG121), EGc9 (148ENPHGVI154), EGc10 (179KVFHPMS185), EGc11 (230QAGMGVVG237), EGc12 (303RSHDSQGKIS312), and EGc13 (430DIPRFV435). Amino acid sequence alignment of the BCEs with homologous proteins from 10 closely related SFTSV strains showed that EGc4, EGc5, EGc8, EGc11, and EGc12 were the most conserved epitopes (i.e., with 100% sequence identity) and could be used as candidates for the future design of broad-spectrum multi-epitope vaccines against GTV and SFTSV. EGc1, EGc2, EGc3, EGc6, EGc7, EGc9, and EGc13 differed more significantly from the SFTSV Gc sequences, showing differences of 1–2 amino acids, and thus could be used for the serological verification of GTV infection in patients in the future. Since no GTV-infected patients have been identified thus far, it is not possible to confirm whether these epitopes are specific to this virus. In the future, detection of serum antibodies to SFTSV and GTV simultaneously in patients will remain necessary to compare the specificity and sensitivity of these epitopes and to establish a theoretical basis for the development of broad-spectrum and specific detection methods and vaccines for bunyaviruses.

Additionally, the 13 BCEs identified in this study were all confirmed by secondary and 3D structural analyses to be located on the surface of the protein structure and showed good antigenicity. The prediction results included most of the experimentally identified epitope motifs, indicating that the epitope prediction tool combined with the BSP method was a reliable approach for epitope identification and mapping and could reduce the experimental workload and costs of epitope mapping and immunodiagnosis [14, 26, 27]. Moreover, none of the 13 BCEs required purification, mainly because the recombinant fusion peptide expressed by the pXXGST-3 vector was larger than the protein in the empty vector, and the expression products were located in the weak antigenic region of E. coli. Therefore, the expression products could be used directly for western blot analysis without a purification step, thereby reducing the time needed for screening of the recombinant bacteria and the workload [23, 24].

Among study limitations, this study only identified linear epitopes of the Gc protein. Also, there are currently no B cell epitopes predicted to be affected by different MHC haplotypes. The mechanism of virus neutralization by anti-Gc protein antibodies remains unknown, and neutralization through identified linear epitopes of Gc protein has not been characterized. Notably, improved BSP method adopted is only suitable for screening simple linear epitopes. In the future, peptide microarray high-throughput screening, mass spectrometry, and other techniques may be necessary to identify conformational epitopes in order to avoid missing the detection of specific epitopes of SFTSV and GTV that could be used to distinguish the two viruses. In the future, serological validation of GTV should be conducted to investigate whether these epitopes are specific to this virus and whether they can be used to detect SFTSV and GTV antibodies or to design broad-spectrum vaccines.The epitopes of the GTV Gc protein obtained in this study can serve as fundamental data for future research works aimed at elucidating the immunological properties of the GTV glycoproteins and for developing multi-epitope diagnostic reagents and vaccines.

References

  1. 1. Estrada-Pena , Agustín , de la Fuente . The ecology of ticks and epidemiology of tick-borne viral diseases. Antiviral Res. 2014; 108: 104–128. pmid:24925264
  2. 2. Kurihara S, Satoh A, Yu F, Hayasaka D,Shimojima M, Tashiro M et al. The world first two cases of severe fever with thrombocytopenia syndrome: An epidemiological study in Nagasaki, Japan. J Infect Chemother. 2016; 22(7): 461–465. pmid:27142979
  3. 3. Takahashi T, Maeda K, Suzuki T, lshido A, Shigeoka T, Tominaga T, et al. The first identification and retrospective study of severe fever with thrombocytopenia syndrome in Japan. J Infect Dis. 2014; 209(6): 816–827. pmid:24231186
  4. 4. Fill MA, Compton ML, Mcdonald EC, Moncayo AC, Dunn JR, Schaffner W, et al. Novel clinical and pathologic findings in a Heartland virus-associated death. Clin Infect Dis. 2017; 64(4): 510–512. pmid:27927857
  5. 5. Abudurexiti A, Adkins S, Alioto D, Alkhovsky SV, Avšič-Županc T, Ballinger MJ, et al. Taxonomy of the order Bunyavirales: update 2019. Arch Virol. 2019; 164(7): 1949–1965. pmid:31065850
  6. 6. Shen S, Duan X, Wang B, Zhu L, Zhang Y, Zhang J, et al. A novel tick-borne phlebovirus, closely related to severe fever with thrombocytopenia syndrome virus and Heartland virus, might be a potential pathogen. Emerg Microbes Infect. 2018; 7(1): 95–108. pmid:29802259
  7. 7. Yoo JR, Heo ST, Park D, et al. Family cluster analysis of severe fever with thrombocytopenia syndrome virus infection in Korea. Am J Trop Med Hyg. 2016; 95(6): 1351–1357. pmid:27928083
  8. 8. Yoshikawa T, Shimojima M, Fukushi S, Tani H, Fukuma A, Taniguchi S, et al. Phylogenetic and geographic relationships of severe fever with thrombocytopenia syndrome virus in China, South Korea, and Japan. J Infect Dis. 2015; 212(6): 889. pmid:25762790
  9. 9. Shi J, Hu S, Liu X, Yang J, Liu D, Wu L, et al. Migration, recombination, and reassortment are involved in the evolution of severe fever with thrombocytopenia syndrome bunyavirus. Infect Genet Evol. 2017; 47: 109–117. pmid:27884653
  10. 10. Hofmann H, Li X, Zhang X, Liu W, Kühl A, Kaup F, et al. Severe fever with thrombocytopenia virus glycoproteins are targeted by neutralizing antibodies and can use DC-SIGN as a receptor for pH-dependent entry into human and animal cell lines. J Virol. 2013; 87(8): 4384–4394. pmid:23388721
  11. 11. Halldorsson S, Behrens AJ, Harlos K, T.Huiskonen J, M.Elliott R, Crispin M, et al. Structure of a phleboviral envelope glycoprotein reveals a consolidated model of membrane fusion. Proc Natl Acad Sci U S A. 2016; 113(26): 7154–7159. pmid:27325770
  12. 12. Yuan F, Zheng AH. Entry of severe fever with thrombocytopenia syndrome virus. Virologica Sinica, 2017; 32(001): 44–50. pmid:27995422
  13. 13. Elliott RM, Brennan B. Emerging Phleboviruses. Curr Opin Virol. 2014; 5(100): 50–57. pmid:24607799
  14. 14. Zhu R, Zhen L, Dong SJ, Rui-Song YU, Zhu YM. Progresses on B cell epitope mapping, microbio china. 2017; 44(9): 2172–2183. https://doi.org/10.13344/j.microbiol:china.170234
  15. 15. Potocnakova L, Bhide M, Pulzova LB. An introduction to B-cell epitope mapping and In Silico epitope prediction. J Immunol Res. 2016; 2016: 1–11. pmid:28127568
  16. 16. Abbott WM, Damschroder MM, Lowe DC. Current approaches to fine mapping of antigen-antibody interactions. Immunology. 2014; 142(4): 526–535. pmid:24635566
  17. 17. Zhang J, Moming A, Yue X, Shen S, Liu D, Xu WX, et al. Fine epitope mapping of glycoprotein Gn in Guertu virus. PLoS One. 2019; 14(10): e0223978. pmid:31618247
  18. 18. Morrow JF, Cohen SN, Chang AC, Boyer HW, Goodman HM, Helling RB. Replication and transcription of eukaryotic DNA in Escherichia coli. Proc Natl Acad Sci U S A. 1974; 71(5): 1743–1747. pmid:4600264
  19. 19. Houen G. Peptide antibodies: Past, present, and future. Methods Mol Biol. 2015; 1348: 1–6. pmid:26424257
  20. 20. Ladner RC. Mapping the epitopes of antibodies. Biotechnol Genet Eng Rev. 2007; 24: 1–30. pmid:18059626
  21. 21. Han Z, Zhao F, Shao Y, Liu X, Kong X, Song Y, et al. Fine level epitope mapping and conservation analysis of two novel linear B-cell epitopes of the avian infectious bronchitis coronavirus nucleocapsid protein. Virus Res. 2013; 171(1): 54–64. pmid:23123213
  22. 22. Davidson E, Doranz BJ. A high-throughput shotgun mutagenesis approach to mapping B-cell antibody epitopes. Immunology. 2014; 143(1): 13–20. pmid:24854488
  23. 23. Xu WX, He YP, Tang HP, Jia XF, Ji CN, Gu SH, et al. Minimal motif mapping of a known epitope on human zona pellucida protein-4 using a peptide biosynthesis strategy. J Reprod Immunol. 2009; 81(1): 9–16. pmid:19539378
  24. 24. Xu WX, Wang J, Tang HP, He YP, Zhu QX, Gupta SK, et al. Epitomics: IgG-epitome decoding of E6, E7 and L1 proteins from oncogenic human papillomavirus type 58. Sci Rep. 2016; 6: 34686. pmid:27708433
  25. 25. Moming A, Tuoken D, Yue X, Xu WX, Guo R, Liu D, et al. Mapping of B-cell epitopes on the N-terminal and C-terminal segment of nucleocapsid protein from Crimean-Congo hemorrhagic fever virus. PLoS One. 2018; 13(9): e0204264. pmid:30235312
  26. 26. Shalitanati A, Yu H, Liu D, Xu WX, Yue X, Guo R, et al. Fine mapping epitope on glycoprotein-Gn from Crimean-Congo hemorrhagic fever virus. Comp Immunol Microbiol Infect Dis. 2018; S0147957118300572. pmid:30290884
  27. 27. Zhang J, Simayi A, Wang M, Moming A, Xu WX, Wang C, et al. Fine mapping epitope on glycoprotein Gc from Crimean-Congo hemorrhagic fever virus. Comp Immunol Microbiol Infect Dis. 2019; 67: 101371. pmid:31627038
  28. 28. Abula A, Shen S, Zhang J, Liu X, Li Y, Deng F, et al. Expression and antibody preparation of recombinant truncated glycoprotein of Guertu virus. Chin J Microbiol Immunol. 2020; 40(03): 178–184. https://doi.org/10.3760/cma.j.cn112309-20190727-00231
  29. 29. Xu WX, Wang J, Tang HP, Chen LH, Lian WB, Zhan JM, et al. A simpler and more cost-effective peptide biosynthetic method using the truncated GST as carrier for epitope mapping. PLoS One. 2017; 12(10): e0186097. pmid:29023483
  30. 30. Garnier J, Gibrat JF, Robson B. GOR method for predicting protein secondary structure from amino acid sequence. Methods Enzymol. 1996; 266: 540–53. pmid:8743705
  31. 31. Chou PY, Fasman GD. Prediction of the secondary structure of proteins from their amino acid sequence. Adv Enzymol Relat Areas Mol Biol. 1978; 47(6): 145–148. pmid:364941
  32. 32. Plegge T, Hofmann-Winkler H, Spiegel M, Pöhlmann S. Evidence that processing of the severe fever with thrombocytopenia syndrome virus Gn/Gc polyprotein is critical for viral infectivity and requires an internal Gc signal peptide. PLoS One. 2016; 11(11): e0166013. pmid:27855227
  33. 33. Mansfield KL, Banyard AC, McElhinney L, Johnson N, Horton DL, Hernández-Triana LM, et al. Rift Valley fever virus: A review of diagnosis and vaccination, and implications for emergence in Europe. Vaccine. 2015; 33(42): 5520–5531. pmid:26296499
  34. 34. Bente DA, Forrester NL, Watts DM, McAuley AJ, Whitehouse CA, Bray M. Crimean-Congo hemorrhagic fever: history, epidemiology, pathogenesis, clinical syndrome and genetic diversity. Antiviral Res. 2013; 100(1): 159–89. pmid:23906741
  35. 35. Xu B, Liu L, Huang X, Ma H, Zhang Y, Du Y, et al. Metagenomic analysis of fever, thrombocytopenia and leukopenia syndrome (FTLS) in Henan Province, China: discovery of a new bunyavirus. PLoS Pathog. 2011; 7(11): e1002369. pmid:22114553
  36. 36. Bao CJ, Guo XL, Qi X, Hu JL, Zhou MH, Varma JK, et al. A family cluster of infections by a newly recognized bunyavirus in eastern China, 2007: further evidence of person-to-person transmission. Clin Infect Dis. 2011; 53(12): 1208–14. pmid:22028437
  37. 37. Liu Y, Li Q, Hu W, Wu J, Wang Y, Mei L, et al. Person-to-person transmission of severe fever with thrombocytopenia syndrome virus. Vector Borne Zoonotic Dis. 2012; 12(2): 156–60. pmid:21955213
  38. 38. Tang X, Wu W, Wang H, Du Y, Liu L, Kai K, et al. Human-to-human transmission of severe fever with thrombocytopenia syndrome bunyavirus through contact with infectious blood. J Infect Dis. 2013; 207(5): 736–9. pmid:23225899
  39. 39. Hinkula J, Devignot S, Åkerström S, Karlberg H, Wattrang E, Bereczky S, et al. Immunization with DNA plasmids coding for Crimean-Congo hemorrhagic fever virus capsid and envelope proteins and/or virus-like particles induces protection and survival in challenged mice. J Virol. 2017; 91(10): e02076–16. pmid:28250124
  40. 40. Yu L, Bai W, Wu X, Zhang L, Zhang L, Li P, et al. A recombinant pseudotyped lentivirus expressing the envelope glycoprotein of Hantaan virus induced protective immunity in mice. Virol J. 2013; 10(1): 301. pmid:24093752
  41. 41. Arikawa J, Yao JS, Yoshimatsu K, Takashima I, Hashimoto N. Protective role of antigenic sites on the envelope protein of Hantaan virus defined by monoclonal antibodies. Arch Virol. 1992; 126(1–4): 271–81. pmid:1381911