Entao Li1†
Qiuxue Han1,2†
Jinhao Bi1,3†Shimeng Wei1,3Shen Wang1Ying Zhang1,4Jun Liu1Na Feng1Tiecheng Wang1Jun Wu5
Songtao Yang1Yongkun Zhao1*Bo Liu5*
Feihu Yan1*Xianzhu Xia1*- 1Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun, China
- 2College of Veterinary Medicine, Jilin Agriculture University, Changchun, China
- 3Institute of Laboratory Animal Science, Chinese Academy of Medical Science and Comparative Medicine Center, Peking Union Medical College, Beijing, China
- 4College of Wildlife and Protected Area, Northeast Forestry University, Harbin, Heilongjiang, China
- 5Department of Microorganism Engineering, Beijing Institute of Biotechnology, Beijing, China
The emergence of SARS-CoV-2 variants stresses the continued need for broad-spectrum therapeutic antibodies. Several therapeutic monoclonal antibodies or cocktails have been introduced for clinical use. However, unremitting emerging SARS-CoV-2 variants showed reduced neutralizing efficacy by vaccine induced polyclonal antibodies or therapeutic monoclonal antibodies. In our study, polyclonal antibodies and F(ab’)2 fragments with strong affinity produced after equine immunization with RBD proteins produced strong affinity. Notably, specific equine IgG and F(ab’)2 have broad and high neutralizing activity against parental virus, all SARS-CoV-2 variants of concern (VOCs), including B.1.1,7, B.1.351, B.1.617.2, P.1, B.1.1.529 and BA.2, and all variants of interest (VOIs) including B.1.429, P.2, B.1.525, P.3, B.1.526, B.1.617.1, C.37 and B.1.621. Although some variants weaken the neutralizing ability of equine IgG and F(ab’)2 fragments, they still exhibited superior neutralization ability against mutants compared to some reported monoclonal antibodies. Furthermore, we tested the pre-exposure and post-exposure protective efficacy of the equine immunoglobulin IgG and F(ab’)2 fragments in lethal mouse and susceptible golden hamster models. Equine immunoglobulin IgG and F(ab’)2 fragments effectively neutralized SARS-CoV-2 in vitro, fully protected BALB/c mice from the lethal challenge, and reduced golden hamster’s lung pathological change. Therefore, equine pAbs are an adequate, broad coverage, affordable and scalable potential clinical immunotherapy for COVID-19, particularly for SARS-CoV-2 VOCs or VOIs.
Introduction
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a causative agent of coronavirus disease 2019 (COVID-19) discovered in December 2019 (1). SARS-CoV-2 is a single-stranded enveloped, positive-sense RNA virus. The virion is mainly composed of membrane protein (M), and nucleocapsid protein (N), spike protein (S), and envelope protein (E) (2). The S protein is a highly glycosylated type I membrane protein composed of S1 and S2 subunits on the surface of the virus particle. It is displayed on the surface of the virus particle in the form of a homotrimer, forming a unique coronavirus S structure that mediates virus entry (3). The receptor binding domain (RBD) is located on the S1 subunit, which mediates SARS-CoV-2 binding to the angiotensin-converting enzyme 2 (ACE2) receptor of the host cell (4) and acts as an antigenic epitope inducing neutralizing antibodies (5, 6). The function of the transmembrane unit S2 is to promote membrane fusion between the virus and the host cell (7, 8). The above characteristics of the S protein make it an important target for the development of vaccines and antibodies against SARS-CoV-2 (9). In the course of the COVID-19 pandemic, the host spectrum of the virus has been expanding and the cellular microenvironment has been different. This can lead the virus to make a series of adaptive mutations to adapt to their hosts, which influence its virulence, transmissibility, or pathogenesis (10). Numerous SARS-CoV-2 variants have already been documented globally during the COVID-19 pandemic, classified as VOC and VOI, and high consequence (11). Five variants have rapidly become dominant in several countries, raising concerns as they possess mutations of interest and providing evidence of international spreading: B.1.1,7 (Alpha), B.1.351 (Beta), B.1.617.2 (Delta), P.1 (Gamma), and B.1.1.529 (Omicron) (12). The strike of the increasing number of confirmed cases and deaths globally is unimaginable. Although a variety of therapeutic antibody drugs, such as homologous or heterologous monoclonal antibodies (13–15) or their cocktails (16), nanobodies (17–20), humanized antibodies (21) and chimeric antibodies (22), have been licensed or approved to enter clinical trials, the current epidemic situation and the emergence of SARS-CoV-2 variants have prompted the acceleration of the development of therapeutic antibody drugs with broad-spectrum neutralizing activity (23, 24).
At the early stage of the COVID-19 global pandemic, plasma transfusion from recovered individuals was identified as a potential therapeutic method for COVID-19 in clinical practice (25–29). Obviously, the quantity and quality of plasma were not guaranteed to satisfy the standard of clinical use, limiting the widespread use of this treatment method (30–32). Immunoglobulin purified from COVID-19 convalescent plasma showed higher activity than the original plasma (33). However, treatments relying on human plasma have obvious and neglected shortcomings. Blood sources are highly dependent on human volunteers, and human pathogens have to be screened. It is also worth noting that the nucleic acid of SARS-CoV-2 can be detected in the plasma, and the risk of transmission through blood transfusion still exists. Fortunately, animal-derived immunoglobulins can solve the problem of a lack of blood sources. Antiserum produced after active immunization of large animals such as horses has long been approved for clinical application to treat complex and difficult diseases, especially venomous snake bites and various highly pathogenic infectious diseases (34–37).
With no licensed antibodies with broad-spectrum neutralizing activity against various SARS-CoV-2 variants, there is an urgent, cross-reactive and affordable therapeutic approach to improve the prognosis of COVID-19 patients and reduce person-to-person spread. Aninsworth S et al. proposed that animal-derived antibodies should be considered alongside convalescent human plasma to deliver treatments for COVID-19 (38). Therefore, another promising therapy for COVID-19 patients is the intravenous administration of heterologous pAbs, which are purified from the plasma of horse hyperimmunized with SARS-CoV-2-related antigens. In this study, we described an application of equine anti-SARS-CoV-2 purified polyclonal IgG and F(ab’)2 with demonstration of their cross-reactive neutralizing activity against all SARS-CoV-2 VOCs and VOIs and evaluated their protective efficacy in a BALB/c and hamster model of COVID-19, respectively.
Results
Generation of the recombinant RBD protein of SARS-CoV-2
To effectively induce a robust nAb response against SARS-CoV-2 in vivo, the RBD of the S protein was selected as the immunogen (Supplementary Data, Figure S1A). Based on our previous study (39), honeybee melittin (HBM) signal peptide was added to the N-terminal of the RBD region to promote the expression of SARS-CoV-2 RBD protein in secretory forms. Baculovirus expression system and 6×His tags were applied for mass expression and purification of the recombinant RBD protein. As shown in Supplemental Figure 1, Western blot results suggested that the recombinant RBD protein was represented by a 35 kDa band (Figure S1C). Meanwhile, the purity of the recombinant RBD preparation, as determined by SDS‒PAGE followed by thin layer chromatography, was determined to be 90% (Figure S1D). The above results show that the RBD protein of SARS-CoV-2 was well prepared for further research.
Equine immunization and production of equine antisera products
Three horses (No. 15, No. 16, and No. 17) were immunized with adjuvanted RBD protein purified from Sf9 cells. The results showed that the neutralizing antibody titers of the sera from the three horses all showed a trend of an unremitting increase in the number following a series of immunizations. The neutralizing serum titer was detected with the SARS-CoV-2 Wuhan01 strain in the BSL-3 laboratory. The neutralizing antibody titers of the three horses were 1:34,000, 1:51,200, and 1:25,600 after five immunizations, respectively (Figure S2B). In view of the above results, the sera of horses No. 15 and No. 16 after five immunizations were selected for subsequent IgG and F(ab’)2 purification. The saturated ammonium sulfate precipitation method was applied to purify IgG from equine serum, and pepsin was used to digest the Fc region to obtain purified F(ab’)2 fragments. Finally, the sample was run through a Protein A column to remove the residual immunoglobulin G, and the purified F(ab’)2 fragments were collected by running through a Protein-G column. The integrity of IgG and F(ab’)2 fragments was determined by SDS‒PAGE (Figures S2C, D). After purification, the purities of IgG and F(ab’)2 were determined to be 92.6% (Figure S1E) and 93.7%, respectively (Figure S2F). These results collectively demonstrated that the purity of IgG and F(ab’)2 products were obtained from equines.
In vitro characterization of equine IgG and F(ab’)2
Neutralizing activity against SARS-CoV-2 is an indicator of functional antibodies that confer a protective immune response. The binding activity to RBD and neutralizing activity of equine IgG and F(ab’)2 were investigated. The RBD protein produced by the prokaryotic expression system was used to detect the binding ability of purified IgG. The results showed that the ELISA titer of IgG was 1:12,280 (Figure 1B). Hyperimmune serum, purified IgG and F(ab’)2 derived from the No. 15 and No. 16 equines were selected to evaluate whether the purification operation affected the neutralizing ability. As shown in Figure 1A, purified IgG and F(ab’)2 of equine No. 15 were found to possess similar neutralizing activity levels, with a neutralizing titer of over 50,000, while the titer of equine No. 15 serum products was over 40,000. The neutralizing titer of serum and purified IgG derived from equine No. 16 was over 30,000, but the purified F(ab’)2 has a slightly decreased neutralizing titer of over 20,000. No significant differences in neutralizing activity between hyperimmune serum and purified equine immunoglobulins were found after purification. These results indicated that the preparation and purification processes do not adversely affect the neutralizing capacities of these antibodies.
Figure 1 In vitro characterization of purified equine immunoglobulin against SARS-CoV-2. (A) The neutralizing titers of hyperimmune serum, purified IgG, and F(ab’)2 derived from equine No. 15 and No. 16 were tested with wild type SARS-CoV-2 Wuhan 01. The serum neutralizing antibody titer was defined as the reciprocal of the highest dilution showing a 100% CPE reduction compared to the virus control. (B) The titers of purified SARS-CoV-2-specific IgG in equine sera were examined via RBD-capture ELISA. Two repeated tests were performed on each sample.
Broad-spectrum neutralizing activity test
To evaluate the broad-spectrum neutralizing activity of equine immunoglobulin, VSV-vectored pseudoviruses of the five SARS-CoV-2 mutants recommended by the WHO were applied. Table S1 shows the source of each variant of VOCs and VOIs and the specific information of mutation sites on the S protein. It is not difficult to find that most variants of VOCs and VOIs shared D614G mutation sites since the D614G mutation occurred in the B.1.1.7 variant. Additionally, the N501Y mutation was identified in most SARS-CoV-2 VOCs except for the B.1.617.2 variant. N501Y and D614G mutations were stable in the Omicron variant (BA.1 and BA.2), which has recently continued to circulate around the world (Table S1). The results of the neutralization assay suggested that IgG or F(ab’)2 derived from equine No. 15 and No. 16 have similar neutralizing activity against SARS-CoV-2 Wuhan01 (Figures 2A, B). The mean neutralizing titer of IgG from the two equines was 1:51,200, while the mean neutralizing titers of F(ab’)2 were 51,200 and 42,000, respectively. However, the SARS-CoV-2 VOCs defined by the WHO, including Alpha B.1.1,7 (Alpha), B.1.351 (Beta), B.1.617.2 (Delta), P.1 (Gamma), and B.1.1.529 (Omicron), weakened the neutralizing effect of IgG and F(ab’)2 to varying degrees. Compared to SARS-CoV-2 Wuhan01, the VOCs weakened the neutralizing titer by 1.5~4.0 folds of IgG and 2.4~8.0 folds of F(ab’)2 for No.15 equine (Figure 2A), 2.4~9.6 folds of IgG and 2.0~13.3 folds of F(ab’)2 for No.16 equine (Figure 2C). VOIs weakened the neutralizing titer by 1.2~3.0 folds of IgG and 1.5~4.8 folds of F(ab’)2 for No.15 equine (Figure 2B), 1.2~3.0 folds of IgG and 1.7~6.7 folds of F(ab’)2 for No.16 equine (Figure 2D). Among the VOCs, regardless of whether they were IgG or F(ab’)2, Omicron had the most significant neutralization weakening ability, with a mean neutralizing titer of 12,800 for No. 15 (Figure 3A) and 4,200 for No. 16 (Figure 3C). IgG and F(ab’)2 products still maintained a strong neutralizing effect. These results indicated that equine polyclonal antibodies against SARS-CoV-2 have potential broad-spectrum neutralizing activity against SARS-CoV-2 variants. For the VOIs, IgG and F(ab’)2 of No. 15 equine have a neutralizing capability against pseudoviruses of B1,1,42 (Epsilon) and P.2 (Zeta) similar to that of the parental virus SARS-CoV-2 Wuhan01. Meanwhile, other variants of VOIs to a lesser extent weakened the neutralizing ability (Figure 2C). In comparison to equine No. 15, the neutralizing ability of equine No. 16 IgG or F(ab’)2 was slightly weaker than that of equine No. 15 (Figure 2D). In the following study, we selected the equine serum product of horse No. 15 to evaluate the protective efficacy in rodent models.
Figure 2 Broad-spectrum neutralizing activity test against SARS-CoV-2 VOC and VOI. The neutralizing antibody titers were calculated as the highest dilution of sera that completely inhibited virus-caused CPE. The serum neutralizing antibody titer was defined as the reciprocal of the highest dilution showing a 100% CPE reduction compared to the virus control. (A) Neutralizing antibody titers of purified IgG and F(ab’)2 of equine No.15 against SARS-CoV-2 VOC; (B) Neutralizing antibody titers of purified IgG and F(ab’)2 of equine No.16 against SARS-CoV-2 VOC; (C) Neutralizing antibody titers of purified equine immunoglobulin of equine No.15 against SARS-CoV-2 VOI; (D) Neutralizing antibody titers of purified equine immunoglobulin of equine No.16 against SARS-CoV-2 VOI. Comparison to the wild type SARS-CoV-2 Wuhan01, the number above the column represented the fold by which the neutralizing titer of the IgG or F(ab’)2 was weakened by the SARS-CoV-2 VOC and VOI. Samples were processed in triplicate, and error bars indicate standard error. Data are presented as the mean ± SEM. (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001).
Figure 3 Evaluation of the protective efficacy of purified equine immunoglobulin in a mouse model. Groups of 13 BALB/c mice were administered with IgG or F(ab’)2 at 1 day before mouse-adapted SARS-CoV-2 (BMA8) infection or 1 dpi with BMA8. Each mouse was given 250 µg of antibody at a dose of 10 mg/kg. BALB/c mice were challenged intranasally with a lethal dose 50 LD50 of BMA8 before treatment or after administration. The survival rate, weight change, body temperature and clinical scores of BALB/c mice were monitored daily after SARS-CoV-2 BMA8 infection. (A) Schematic diagram of the administration of equine immunoglobulin drugs and virus challenge procedure; (B) Survival rate. (C) Percent weight change. (D) Body temperature change. Body weight change of mice in a with comparison to isotype control was measured by repeated measurements two-way analysis of variance (ANOVA) with Tukey’s post hoc test. Data are mean ± s.e.m. of each experimental group. (****P < 0.0001).
Protection efficacy test of equine immunoglobulin IgG and F(ab’)2 in a mouse model
The half-lives of IgG and F(ab’)2 were assessed in guinea pigs and found to be 5~6 days and 24~48 h, respectively (Table 1). To verify the protective effect of equine antibodies in vivo, purified IgG and F(ab’)2 were given to wild-type (WT) BALB/c mice by intraperitoneal injection with 250 µg per dose (10 mg/kg) at 1 day post-infection (dpi) or 1 day preinfection with a lethal dose mouse-adapted SARS-CoV-2 (Figure 3A). The results indicated that all BALB/c mice injected with IgG or F(ab’)2 survived after SARS-CoV-2 infection, while the prevention control and treatment control succumbed to the infection at 4~7 dpi and 5~7 dpi, respectively (Figure 3B). Throughout the challenge experiment, all mice in control group had obvious weight loss. At the endpoint of the study, the weight of the surviving mice recovered to their original weight, and the weight of the IgG treatment group increased by 1.50 ± 3.04% (mean ± SD). The body weight of the two control groups dropped to 80% of the initial weight within 2~5 dpi (Figure 3C). The mice were euthanized if their body weight was reduced by 20% based on the requirements of the Animal Experimentation Committee. The results of temperature changes showed that the body temperature of mice decreased continuously after SARS-CoV-2 infection, and the temperature changes in the control group were 5 ± 1.55°C in the prevention control group and 5 ± 1.14°C in the therapeutic control group. The body temperature of the prevention group and therapeutic group fluctuated within a normal range. The temperature changes were 0.14 ± 1.09°C in the IgG prevention group and 0.44 ± 0.62°C in the F(ab’)2 prevention group; the temperature changes were 2.28 ± 0.82°C in the IgG treatment group and 1.56 ± 0.67°C in the F(ab’)2 treatment group (Figure 3D).
Table 1 Determination of half-life of purified equine IgG and F(ab’)2 in guinea pigs.
Hematology test
Hematological analysis revealed that BMA8-infected BALB/c mice in the control group showed a moderate increase in white blood cell (WBC) count (Figure 4A), with significant decreases in the percentage of lymphocytes (LYM%) (Figure 4C) and monocytes (Mon%) (Figure 4E) and a decrease in platelet count (PLT) (Figure 4D) at 3~5 dpi in comparison to all treatment groups. A marked increase in the percentage of neutrophils (Neu%) (Figure 4B) was also observed in the control group rather than the treatment groups injected with IgG or F(ab’)2, which is a typical severe characteristic of SARS-CoV-2 infection observed in clinical COVID-19 patients.
Figure 4 Blood counts in SARS-CoV-2-infected mice. The hematological values of BALB/c mice were analysed, including lymphocyte (LYM), neutrophil percentage (Neu%), monocytes (Mon), platelet count (PLT) and white blood cell count (WBC), at 3 dpi after SARS-CoV-2 BMA8 infection. Four infected mice were sacrificed at 3 dpi to collect the whole blood for blood counts test. (A) White blood cell (WBC) count; (B) neutrophil (Neu) percentage; (C) lymphocyte (LYM) percentage; (D) platelet (PLT) (E) Monocyte(Mno). Data are presented as the mean ± SEM (n=4). (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001).
Viral load
To assess the bio-distribution of SARS-CoV-2 infection in the upper and lower respiratory tracts, the turbinates and lungs were harvested at selected times (3 dpi and 5 dpi) after infection, and viral RNA (vRNA) was quantified by RT-qPCR and TCID50. As shown in Figure 5, there were much higher levels of vRNA loads in the lung and turbinate at 3 and 5 dpi in both the prevention control group and the treatment control group than in all administration groups. The vRNA loads reached 105~1011 copies/g at 3 dpi (Figure 5A) and 107~1011 copies/g at 5 dpi (Figure 5B) in the prevention control group, while the viral loads reached 108~1012 copies/g at 3 dpi (Figure 5C) and 107~1011 copies/g at 5 dpi (Figure 5D) in the treatment control group. Regardless of the prevention group or the treatment group, no significant differences in viral copies between IgG and F(ab’)2 were found at 3 dpi or 5 dpi. The TCID50 results tended to coincide exactly with the RT-qPCR results, and significantly higher viral titers were confirmed in all control groups (Figures 5E–H). Although the treated mice were positive for vRNA in the turbinates and lungs, infectious viruses were not recovered from those tissues from most of the surviving animals at 5 dpi (Figures 5F, H). Taken together, the above results illustrated that equine IgG or F(ab’)2 treatment significantly inhibited viral replication in the upper and lower tracts, which may block shedding of infectious SARS-CoV-2.
Figure 5 Viral loads in the upper and lower respiratory tracts of mice. Four infected mice were sacrificed at 3 dpi and 5 dpi, respectively, and the turbinate and lungs were harvested to analyze the viral RNA loads by RT‒qPCR and TCID50. (A) The viral load quantified by RT‒qPCR at 3 dpi in the prevention group. (B) The viral load quantified by RT‒qPCR at 5 dpi in the prevention group. (C) The viral load quantified by RT‒qPCR at 3 dpi in the treatment group. (D) The viral load quantified by RT‒qPCR of the turbinate and lungs at 5 dpi in the treatment group. (E) The viral load confirmed by TCID50 at 3 dpi in the prevention group; (F) The viral load confirmed by TCID50 at 5 dpi in the prevention group; (G) The viral load confirmed by TCID50 at 3 dpi in the treatment group; (H) The viral load confirmed by TCID50 at 5 dpi in the treatment group. Data are presented as the mean ± SEM (n=4). (**P < 0.01, ***P < 0.001, ****P < 0.0001).
Histopathology and immunohistochemistry findings in the internal organs of mice
The gross pathology and histological assays were evaluated to check the pathological changes in the lungs at 3 dpi. The pathological section of the lung from controlled animals showed that the tissue airway epithelial cells were tightly arranged, and there was no obvious epithelial cell degeneration and necrosis; more necrotic epithelial cells (blue arrows) can be found in individual airway cavities in the local tissues (Figure 6E), and the number of airway epithelial cells were reduced; a large range of alveolar cavities were narrowed, and a small amount of alveolar wall capillary congestion (black arrow), and a small amount of neutrophil infiltration was obvious and perivascular edema with a handful of inflammatory cell infiltration (yellow arrow) in the local alveolar cavity (Figures 6A, B, F); spleen result prompted that the red pulp and white pulp were clearly demarcated and spotted apoptosis of lymphocytes, nuclear pyknosis and deep staining or fragmentation (black arrows) were shown in the spleen nodules, and the expansion of germinal centers (yellow arrows) could be seen locally (Figures 6C, G); scattered neutrals were mostly seen in the red pulp granulocyte infiltration (red arrow) (Figure 6D); more brown‒yellow particles (blue arrow) were found in the red pulp (Figures 6G, H).
Figure 6 Histopathological and immunohistochemistry findings in SARS-CoV-2-infected mice. The lungs and spleens were collected from the control mice infected with SARS-CoV-2 without equine immunoglobulin drug injection at 3dpi, and the lungs, spleens, livers and kidneys were harvested from recovered mice. After each tissue was embedded in paraffin, the sections were sectioned for HE staining. (A, B, E, F) Lung tissue changes of control mice were characterized by more necrotic epithelial cells (blue arrow), a small amount of neutrophil infiltration, and perivascular edema with a small amount of inflammatory cell infiltration in the local alveolar cavity (yellow arrow). (C, D, G, H) Spleen tissue changes of control mice were characterized with spotted apoptosis of lymphocytes, nuclear pyknosis and deep staining or fragmentation in the spleen nodules (black arrows), and the expansion of germinal centers (yellow arrow), scattered neutrophils mostly seen in the red pulp granulocyte infiltration (red arrow), and more brown‒yellow particles in the red pulp (blue arrow). (I-L) The basically normal structure of the lung, spleen liver, and kidney tissues were found in administration groups given equine IgG or F(ab’)2. The figure showed immunohistochemistry (IHC) labeling against SARS-CoV-2 N. (M) Viral antigen was not detectable in prevention group given purified IgG; (N) Viral antigen was not detectable in prevention group given purified F(ab’)2; (O) Viral antigen was detected for positive in prevention control group; (P) Viral antigen was not detectable in treatment group given purified IgG; (Q) Viral antigen was not detectable in treatment group given purified IgG F(ab’)2; (R) Viral antigen was detected for positive in treatment control group. (scale bar = 100 μm).
The lungs of the surviving mice showed that the bronchial epithelial structure was intact, the epithelial cells were normal and tightly arranged, and the alveolar structure was clear and not significantly thickened without any obvious inflammation (Figure 6I). The splenic tissue capsule was composed of dense connective tissue with uniform thickness and rich in elastic fibers and smooth muscle fibers. The capsule connective tissue extended into the spleen to form trabeculae, with no obvious abnormalities. Hepatic cords were arranged neatly, with no obvious expansion or squeezing of liver sinusoids (Figure 6J). The liver lobule was clearly demarcated, and the staining nuclei in the liver cells were blue. The nuclei were large, round, and centered, with little heterochromatin and light staining, and the nucleoli could be clearly seen (Figure 6K). The renal tissue was evenly stained, the renal cortex and medulla were clearly demarcated, the glomerulus was normal in shape and structure, and the epithelial cells of renal tubules were closely arranged without obvious inflammation (Figure 6L).
To further identify the effectiveness of purified equine immunoglobulin IgG and F(ab’)2, after combining the observation of pathological sections, we continued to perform a detailed histopathological assessment of the lungs of the mice at 5 dpi after challenge using immunohistochemistry by detecting SARS-CoV-2 nucleocapsid. Regardless of the prevention group or the treatment group, after the application of IgG or F(ab’)2, the mouse lungs were normal and no virus was detectable (Figures 6M, N, P, Q). SARS-CoV-2 nucleocapsid was detected in the lung tissue by immunohistochemistry (IHC) staining, especially in conducting airway epithelia and in the alveoli (Figures 6O, R). Based on the above histopathological results, equine immunoglobulin IgG and F(ab’)2 could effectively clear SARS-CoV-2 in the lungs of mice.
Protection efficacy test of equine immunoglobulin IgG and F(ab’)2 in the golden hamster model
The immunocompetent Syrian golden hamster has also been used to evaluate equine pAb activity against SARS-CoV-2 infection in the upper and lower respiratory tracts (40, 41). We used this animal model to independently assess the inhibitory activity of equine IgG or F(ab’)2. Weights were followed for 6 days, and then tissues were harvested for virological, histopathology and immunohistochemistry analysis. The treatment strategy is shown in Figure 7A. At 500 μg of equine immunoglobulins, hamsters treated with equine immunoglobulin IgG or F(ab’)2 and infected with wild-type SARS-CoV-2 Wuhan01 showed a 100% survival rate (Figure 7B), protection against weight loss (Figure 7C), and reduced viral burden levels in the lungs and nasal swabs compared to the independent control (Figures 7D, E). Correspondingly, RT‒qPCR and TCID50 analysis showed reduced SARS-CoV-2 RNA replication in the lungs of hamsters treated with equine immunoglobulin IgG or F(ab’)2. Although treated animals were positive for vRNA in turbinates and lungs by RT-qPCR (Figures 7D, E), infectious viruses were not recovered from most of these animals via TCID50 analysis (Figures 7F, G). Overall, when a 500 μg dose of IgG or F(ab’)2 was administered, equine pAbs established a trend toward protection against weight loss in hamsters infected with SARS-CoV-2 Wuhan01.
Figure 7 Evaluation of the protective efficacy of purified equine immunoglobulin in the golden hamster model. Each golden hamster was given 500 µg of antibody at a dose of 10 mg/kg. Groups of golden hamsters were infected intranasally with 1,000 TCID50 of wild-type SARS-CoV-2 Wuhan 01 before treatment and or after administration. The survival rate and weight change of BALB/c mice were monitored daily after SARS-CoV-2 Wuhan01 infection. Four infected golden hamsters in each group were sacrificed at 3 dpi, and the turbinate and lung samples were collected to analyze the viral RNA loads by RT‒qPCR and TCID50, respectively. (A) Schematic diagram of the administration of equine immunoglobulin drugs and virus challenge procedure. (B) Survival rate. (C) Percent weight change; Body weight change of mice in a with comparison to isotype control was measured by repeated measurements two-way analysis of variance (ANOVA) with Tukey’s post hoc test. Data are mean ± s.e.m. of each experimental group. (D) The viral loads of turbinate were quantified by RT‒qPCR at 3 dpi in each group; (E) The viral loads of lung were quantified by RT‒qPCR at 3 dpi in each group; (F) The viral loads of turbinate were determined by TCID50 at 3 dpi in each group; (G) The viral loads of lung were determined by TCID50 at 3 dpi in each group. Data are presented as the mean ± SEM (n=5). (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001).
For the surviving hamster, the structure of the lung tissue was basically normal, and the alveolar wall in the visual field was uniform without thickening. The alveolar outline was clear with normal size. The alveolar epithelial cell nucleus was round with no necrosis; no inflammatory cell infiltration was found in the lung parenchyma (Figures S3A, E). For the control group, the lung tissue structure was severely abnormal (Figures S3B, F), and proliferation of alveolar epithelial cells is shown by the yellow arrow (Figure S3B). A small amount of inflammatory cell infiltration and hemorrhage was found in the lung parenchyma, as shown by the red arrow (Figure S3F). Immunohistochemistry results indicated that SARS-CoV-2 was negative in all treatment groups (Figures S3C, G) and was detectable in all control groups (Figures S3D, H).
Discussion
The emergence and dissemination of SARS-CoV-2 and the resulting COVID-19 pandemic triggered a global public health crisis. Although several COVID-19 vaccines have been approved for clinical use, demand far exceeds supply, and access to them is inequitable. Moreover, the appearance of SARS-CoV-2 variants and reports of reinfections associated with immune escape highlight the urgent need for effective and broad coverage COVID-19 therapeutics. Of note, the emergence of antibody-resistant SARS-CoV-2 variants that might limit the therapeutic usefulness of monoclonal antibodies can be mitigated by the use of polyclonal antibodies or monoclonal antibody combinations that simultaneously target distinct neutralizing epitopes. A previous study confirmed that the delta VOC has resistance to therapeutic monoclonal antibodies due to binding impairment to the spike protein, indicating increased immune escape (42). Furthermore, innovative therapies such as mAb cocktails of COVID-19 will probably be too expensive for public use in the face of such large numbers of hospitalized infected people in desperate need of drug treatment. Polyclonal antibodies (pAbs) and convalescent plasma are limited to wide application because it is difficult to simultaneously satisfy two strict donor rigorous health conditions with free blood-borne pathogens and high levels of neutralization anti-SARS-CoV-2 antibodies (43). Additionally, the therapeutic strategy of plasma transfusion from convalescent patients cannot meet the large demand for plasma from severe COVID-19 patients. The new generation of processed and purified equine pAbs containing highly purified IgG or F(ab’)2 fragments was demonstrated to be safe and well tolerated. Equine pAbs are easy to manufacture, allowing fast development and scaling up for treatment (44).
Equine polyclonal immunoglobulins with Fc cleavage have been approved for the treatment of snake envenomation or exposure to rabies virus (45, 46). Polyclonal equine antibodies consisting of a wider range of antibody types can be induced because intact RBDs are applied as immunogens, including all neutralizing and nonneutralizing epitopes. Moreover, RBD without the NTD domain was selected as an immunogen rather than the complete spike protein or the SARS-CoV-2 whole virus in our study to avoid inducing NTD-specific ADE antibodies and the occurrence of antibody-dependent enhancement (ADE). Additionally, optimized hyperimmunization protocols using repeated doses of SARS-CoV-2 RBD protein in the presence of adjuvants likely contributed to the much higher neutralizing potency in comparison to the mAb. It is also worth mentioning that equine pAbs against the SARS-CoV-2 RBD recognize multiple epitopes, which reduce the risk of viral escape. A previous has indicated that a polyclonal approach of equine anti-SARS-CoV-2 F(ab’)2 antibodies that achieved a high level of neutralizing potency against all SARS-CoV-2 variants of concern tested including Omicron BA.1, BA.2, BA.2.12 and BA.4/5 (47). Another study also reported that equine polyclonal antibodies (pAbs) have been assessed in clinical trials in Costa Rica against five globally circulating variants of concern: alpha, beta, epsilon, gamma and delta (48). Although the above studies all evaluated the broad-spectrum neutralizing activity of equine antibodies against SARS-CoV-2 variants in vitro, the protective effect of antibodies in vivo has not been thoroughly explored. In our study, we not only tested neutralizing activity of all SARS-CoV-2 VOCs and VOIs using VSV vectored pseudoviruses, but also evaluated the protective effect of antibodies in vivo using two animal models.
Replication-competent pseudovirus-based vesicular stomatitis virus (VSV) vectors have been used as a useful tool to evaluate the neutralizing activity of various kinds of SARS-CoV-2 antibodies, and the results were similar to those of wild-type SARS-CoV-2 (49, 50). Our neutralizing assay of VOCs based on VSV-vectored pseudovirus also supported the above standpoint, in which the purified equine IgG or F(ab’)2 fragment still has high neutralizing potency against all SARS-CoV-2 variants of concern, although to a lesser extent, the neutralizing potency of those equine antibody products was weakened by the variants. Mutants B.1.351 and P.1 have been shown to have potential threats to reduce the efficacy of certain vaccines by weakening the neutralizing capability (51–53). The Omicron variant, as the currently circulating strain and the most divergent variant, has been reported to reduce the sensitivity of SARS-CoV-2 Omicron to antibody neutralization elicited by booster vaccination (54). Notably, in our study, the evaluation results of the neutralizing activity of pseudoviruses against IgG antibodies and F(ab’)2 showed that IgG has better neutralizing activity against each mutant, especially for SARS-CoV-2 Omicron. However, plasma derived from vaccinated or convalescent individuals or therapeutic mAbs with neutralizing activity has reduced neutralization of multiple SARS-CoV-2 variants efficiently (55–58). A previous study confirmed that the delta VOC has shown resistance to therapeutic monoclonal antibodies due to binding impairment to the spike protein, indicating increased immune escape (42).
In our in vitro study, treatment with purified equine IgG or the F(ab’)2 fragment was shown to be 100% effective before challenge and 100% effective postexposure in a BALB/c mouse model. IgG and F(ab’)2 have a strong capability to inhibit SARS-CoV-2 replication in the upper and lower respiratory tracts in the BALB/c mouse model and hamster model. Many treated animals were positive for vRNA in turbinate and lung samples; however, infectious virus was not recovered from most animals. Detection of vRNA without recovery of infectious virus particles has also been observed with Ebola virus (EBOV) and Nipah virus (59, 60). Therefore, future studies should be performed to confirm how long vRNA may persist despite a lack of viremia and virus shedding. It is interesting to illustrate whether vRNAs remain in circulation for extended periods of time and are generated from ongoing viral replication. It is worth mentioning that SARS-CoV-2-specific non-neutralizing (nnAbs) antibodies in the purified equine IgG-mediated Fc-dependent antibody mechanism may also contribute to survival and virus clearance. Fc-Fcγ receptors are expressed on the surface of various immune effector cells and helper cells, such as monocytes, macrophages and dendritic cells. Activating FcγRs are expressed in various effector cells, such as NK cells and DC cells, which regulate lymphocyte proliferation and differentiation and promote the secretion of cytokines and other inflammatory factors and other immunological functions (61). The protective mechanisms employed by nnAbs depend on engagement with effector cells expressing Fc-Fcγ receptors and consist of antibody-dependent cellular cytotoxicity (ADCC), complement-dependent cytotoxicity (CDC), and antibody-dependent cell-mediated phagocytosis (ADCP) (62). When human neutralizing antibodies were given after infection, intact mAbs reduced SARS-CoV-2 burden and lung disease in mice and hamsters better than loss-of-function Fc variant mAbs. Fc engagement of neutralizing antibodies mitigates inflammation and improves respiratory mechanics (63). A previous study demonstrated that non-neutralizing antibodies against filoviruses, including EBOV and Marburg virus (MARV), mediated protection by Fc-factor function and enhanced the efficacy of other antibodies (64, 65). ADCC-mediated antibodies (ADCC-Abs) can be induced from monovalent inactivated influenza A (H1N1) virus vaccine, and high titers of ADCC-Abs can contribute to lower virus replication and reduced clinical symptoms (62). A similar phenomenon was also found in which nnAbs antibodies elicited by recombinant rabies virus vectored vaccine of Lassa fever virus (LASV) are critical for protection against Lassa fever via the ADCC and ADCP pathways (66). In summary, SARS-CoV-2-specific antibodies could limit infection by directing virion neutralization and/or by targeting infected cells for viral clearance through complement or antibody-mediated cytotoxicity/phagocytosis (67). However, beyond neutralization, whether non-neutralizing antibodies in equine polyclonal antibodies mediate the protective effect of COVID-19 still needs to be further verified, and the broad-spectrum antiviral mechanism of a cocktail of antibody treatment drugs, such as equine polyclonal antibodies, which target multiple epitopes of non-neutralizing antibodies and neutralizing antibodies, also needs to be further explored.
Interestingly, the existence of the Fc fragment of IgG was a “two-edged sword”. On the one hand, we will benefit from the Fcγ receptor pathways in driving antibody-mediated antiviral immunity and exclude the possibility of pathogenic or disease-enhancing effects of Fcγ receptor engagement of anti-SARS-CoV-2 antibodies upon infection (68). On the other hand, equine immunoglobulin is a heterologous protein of the human immune system, and the removal of the Fc fragment also avoids the risk of adverse effects of equine immunoglobulin in the human body, eliminating the major concern of ADE. Our study showed that the removal of the Fc fragment did not impair the neutralizing capabilities and cross-reactivity. The ability of F(ab’)2 to neutralize SARS-CoV-2 variants and the curative ratio in animal models was similar to that of purified IgG. Nevertheless, it is still unknown whether the protection mechanism mediated by the Fc fragment will be impaired following the absence of the Fc fragment of equine IgG.
In conclusion, the devastation caused by the COVID-19 outbreak and emerging causative variants highlights the deficiencies in available and affordable treatment options. Our study highlighted that purified equine IgG and F(ab’)2 fragments were not only a potential therapeutic countermeasure for COVID-19 but also a broad-coverage therapy for low- and middle-income countries to rapidly establish a large-scale manufacturing capacity to prepare antiviral drugs against SARS-CoV-2 variants.
Materials and methods
Animals, viruses and cells
Female BALB/c mice aged 9 months, female hartly guinea pigs aged 6~8 weeks, and female golden hamsters aged 4~6 weeks were purchased from Changchun Yisi Experimental Animal Technology Center. Horses aged 5~6 years were kept in Red Hill Horse Farm and given proper animal welfare. BALB/c mouse-adapted SARS-CoV-2 (BMA8), VSV-vectored pseudoviruses of SARS-CoV-2 (VSV-eGFP-SARS-CoV-2) VOCs and VOIs expressing the reporter gene of enhanced green fluorescent protein (eGFP) were generated and stored in the Biosafety Level 3 Laboratory (BSL-3). Our previous study has reported that BMA8 was obtained by passaging the SARS-CoV-2 Wuhan01 in aged mice for 8 generations, which a mouse-adapted virus cause 100% death in aged BALB/c mice. The mouse model was characterized by high viral loads in the upper and lower respiratory tract, accompanied by lymphopenia, and neutrophilia. Q498H (S), A22D (E) and A36V (E) emerged in structural protein, while five mutations were identified in non-structural protein (69). Wild type SARS-CoV-2 Wuhan01 was gifted from Pro.Chenfeng Qin, an isolation derived from a COVID-19 patient of Wuhan (BetaCov/Wuhan/AMMS01/2020). Spodoptera frugiperda Sf9 (ATCC® CRL-1711™) insect ovary cells were grown in Sf-900™ II serum-free medium (Life Technologies, San Diego, CA, USA) at 27°C at 200 rpm in suspension culture for baculovirus infection and protein expression. Vero E6 cells (ATCC® CRL-1587™) were cultured at 37°C in Dulbecco’s modified Eagle medium (DMEM) (DMEM; GIBCO, Grand Island, NY) containing 10% fetal bovine serum (Gibco, ThermoFisher, MA, USA) with 5% CO2 for BMA8 propagation and titration.
Protein expression and purification
The complete genome sequence of SARS-CoV-2 was obtained from GenBank (ID: NC_045512.2), and codons in the RBD (aa319~541) were optimized for protein expression in Sf9 cells (Sangon Biotech, Shanghai, CN). The RBD region gene was cloned into the pFastBac1 vector to generate the recombinant plasmid pFastBac1-RBD. pFastBac1-RBD was transformed into E. coli DH10Bac competent cells (Life Technologies, San Diego, CA, USA) to produce recombinant bacmid. According to the manufacturer’s instructions, ExpiFectamine Sf transfection reagent (Gibco, Thermo Fisher, MA, USA) was used to transfect recombinant bacmid into Sf9 cells. The recombinant baculovirus rBV-RBD was harvested 4 days after transfection. To generate recombinant RBD protein, Sf9 cells were infected with rBVs in a volume of 300 mL with a cell density of 1~2×107/mL. At 120 h post-infection, cells and debris were removed from the supernatant by centrifugation at 5,000 rpm for 20 min in a JA-10 rotor (Beckman J-20XP, Fullerton, USA). The cell culture supernatant was harvested and subjected to affinity purification through HisPur™ Ni-NTA Spin Columns (Thermo Fisher, MA, USA) to obtain high-purity recombinant RBD protein.
Detection of RBD protein by western blot analyses
For western blot analysis, 10 ug high-purity recombinant RBD protein samples separated through SDS-PAGE were transferred onto Polyvinylidene fluoride (PVDF) membranes (Immobilin-p, Millipore, USA). The membranes were then blocked with 5% milk blocking buffer in Tris-buffered saline with Tween 20 (TBST) for 1 hour. A 1,000-fold dilution of the mouse polyclonal antibody (SinoBiological, CN) were used to detect the presence of RBD protein. A 10,000-fold dilution of the HRP-labeled goat anti-mouse IgG was utilized to identify and combine the primary antibody (Abcam, USA). The membranes were visualized using Western Blue stabilized substrate (Promega, USA).
Indirect immunofluorescence test
Briefly, recombinant baculoviruses-infected Sf9 cells were fixed with 80% cold acetone. Next, the cells were subjected to IFA analysis by using 1,000-fold dilution of mouse monoclonal antibodies against SARS-CoV-2 RBD (SinoBiological, Peking, CN). After incubation with primary antibody, the cells were thoroughly washed with PBS and blocked with PBST. Then, FITC-labeled goat anti-mouse IgG (Millipore, Bedford, USA) was added with 0.3% Evans blue for 1 h as a secondary antibody. The plate was placed under a fluorescence microscope (ZEISS Axio Vert. A1, Germany) to observe and capture the fluorescence.
Equine immunization
Three 5~6-year-old healthy male horses were fed under standard conditions. After emulsifying 1.5 mg of RBD protein in 1 mL PBS with an equal volume of Freund’s complete adjuvant (Thermo Fisher, USA), horses were vaccinated intramuscularly via a subcutaneous multipoint injection route. The subsequent immunization dosage was, respectively given 2.5, 2.5, 3.0, and 4.0 mg of protein in 1 mL PBS mixed with equal volume Freund’s incomplete adjuvant (Thermo Fisher, USA) for each horse via subcutaneous multipoint injection route. Each horse was immunized 5 times in total at 14-day intervals (Figure S2A). One week after each immunization, and the serum was harvested for further use.
SARS-CoV-2-specific antibody detection by ELISA
The RBD protein expressed in E. coli was coated into a 96-well microtiter plate at a final concentration of 1 μg/ml and incubated overnight at 4°C. Subsequently, the plate were blocked with 2% BSA solution for 2 hours at 37°C, 100 μL of serum samples after serial twofold dilutions were added to each well and incubated at 37°C for 2 h. After the cells were thoroughly washed with PBS containing 0.5% Tween-20, 100 μL of 5,000-fold diluted HRP-labeled goat anti-horse antibody was added to each well, and the plates were incubated at 37°C for 1 h. After three washes, 100 μL TMB (3,3’,3,5’-tetramethyl benzidine) was added to each well and incubated for 15 min. Then, 50 μL 0.5 M H2SO4 solution was added to the well to stop the reaction in each well. The optical density (OD) was measured at 450 nm. The titer result was defined as the inverse of the maximum dilution of serum.
Preparation of immunoglobulin and F(ab’)2 fragments
A 1/2 volume of saturated ammonium sulfate solution was added to an equal volume of serum diluted with saline. The mixture was gently stirred for 30 min at room temperature, centrifuged at 5,000 rpm for 20 min, and 1/3 volume of saturated ammonium sulfate solution was added. The mixture was incubated at room temperature for 30 min and separated at 5,000 rpm for 20 min. The precipitate was dialyzed in saline overnight at 4°C to remove salt.
The pH of equine antiserum was adjusted to 3.3 with 1 mol/L HCl. The equine antiserum was digested at 30°C for 2.5 h with pepsin activated by NaAc solution (pH=3.3) at a final concentration of 10,000 IU/mL. The digestion reaction was terminated by adding 1 mol/L NaOH to the phase solution to adjust the pH to 7.2. Subsquently, the above solution was purified by protein a column and then by protein g column. The purified protein was confirmed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). The purity of F(ab’)2 was detected by a thin slice scan.
Purification of SARS-CoV-2 specific equine IgG and F(ab’)2
Immunoaffinity chromatography was used to purify RBD-specific IgG and F(ab’)2. Sepharose 4B gel was utilized to prepare the RBD immunoaffinity column. The gel was resuspended in the coupling solution in which the RBD protein was dissolved and was placed on an overturning shaker overnight at 4°C. Subsequently, the gel was transferred to an empty purification column. After washing with coupling buffer, the gel was resuspended in 2% sodium borohydride (V/V) and incubated at room temperature for 30 min. After washing again, 0.1 mol/L Tris-HCl (pH=8.0) solution was used to resuspend and block at room temperature for 2 h. After washing, the gel was washed alternately with 0.1 mol/L glycine-HCl (pH=2.7) solution and washing buffer. The residual liquid was removed, and 0.02 mol/L PBS was used for equilibration. Total IgG prepared by the precipitation method or F(ab’)2 digested by pepsin was mixed and combined with RBD-4B gel overnight at 4°C. RBD-specific equine immunoglobulin IgG and F(ab’)2 were eluted with 0.1 mol/L glycine-HCl (pH=2.7). Residual glycine-HCl in the eluent was terminated with 0.1 mol/L Tris-HCl (pH=9.0) to prevent antibody inactivation. After purification, RBD-specific equine immunoglobulin IgG and F(ab’)2 were obtained.
SDS-PAGE and thin layer chromatography
The purified immunoglobulin and F(ab’)2 fragment prepared were analyzed by SDS-PAGE. Samples were mixed with SDS-PAGE loading buffer and boiled for 10 min. The prepared samples were electrophoresed through a 10% precast SDS‒PAGE gel. Finally, the SDS-PAGE gel was stained with Coomassie Brilliant Blue staining solution and incubated at room temperature for 30 min. A Tanon 3,500 automatic gel image analysis system (Bio-Equip, Pecking, CN) was used to capture images after decolorization.
Wild type SARS-CoV-2 and pseudovirus neutralization assay variants
We measured the cross-neutralizing activity of equine antibodies by wild type SARS-CoV-2 Wuhan01 and pseudovirus-based neutralizing assay. After the inactivated serum samples were diluted from 1:40 with a serial twofold dilution, 100 μL of diluted serum samples were added to a 96-well plate, and 100 TCID50 wild type SARS-CoV-2 or VSV-eGFP-SARS-CoV-2 was added to each well. 1×104 Vero E6 cells were added to each well after the virus-serum mixture was incubated at 37°C for 1 hour. Cytopathic effect (CPE) was detected under a microscope after the cells were cultured at 37°C for 48~72 hours. The neutralizing antibody titers were defined as the inverse of the highest dilution, which could completely inhibit the production of the cytopathic effect (CPE) compared with viral controls. Virus control and cell control were included in each test.
Half-life test in guinea pigs
Two groups of 5 inbreed hartly guinea pigs were administered an intraperitoneal (IP) injection of l mL purified equine antisera (neutralization titer 1:20,480) or F(ab’)2 (neutralization titer 1:20,480). Blood samples were collected at 0, 1, 4, 8, 12, 24, 48, 72, 96, 120, and 144 h post-injection, and sera were used for determination of neutralization titers against wild-type SARS-CoV-2 Wuhan01. The half-life was defined as the time range post-infection in which during which the serum neutralizing antibody titer decreased by 50%.
Mouse challenge study
Mice were randomized into six groups of 13. Each mouse was given 250 µg of antibody at a dose of 10 mg/kg. In the prevention group, BALB/c mice were first intraperitoneally given 250 µg purified equine IgG or F(ab’)2 and challenged by intranasal administration of 100 TCID50 mouse-adapted SARS-CoV-2 in 50 µL DMEM after 24 h. In the treatment group, BALB/c mice were first challenged with 100 TCID50 mouse-adapted SARS-CoV-2 in 50 µL DMEM and given 250 µg purified equine IgG or F(ab’)2 via IP injection 24 h post-infection. At 3 dpi and 5 dpi, four mice in each group were sacrificed. Lungs and turbinates were harvested to test the viral RNA loads, histopathology and immunohistochemistry test. The challenge strain of mouse-adapted SARS-CoV-2 named BMA8 was previously described (69). All mice were monitored for clinical signs of disease, survival rate, weight change and temperature changes.
Viral load determination by RT-qPCR and TCID50
Briefly, viral RNA in each tissue homogenates sample was extracted with an RNeasy Mini kit (QIAGEN, USA) following the manufacturer’s protocol. The viral RNA quantification was performed by RT-qPCR targeting the N gene of SARS-CoV-2. RT-qPCR was performed using One Step PrimeScript RT-PCR Kit (Takara, Dalian, China) with the following primers and probes: NF (5’-GGGGAACTTCTCCTGCTAGAAT-3’); NR (5’-CAGACATTTTGCTCTCAAGCTG-3’); and NP (5’-FAM-TTGCTGCTGCTTGACAGATT-TAMRA-3’). The details can be referred to our previous study (69).
The supernatant of tissue homogenate was serially diluted at a tenfold ratio, and 100μL were added into 96-well cell culture plates that had been plated with Vero E6 cells. The plates were incubated at 37°C for 1 h and then washed 3 times with sterile PBS. 100μL of incomplete medium containing 2% FBS and 1% penicillin‒streptomycin was added to each well of 96-well plate. The median tissue culture infective dose (TCID50) was calculated by observing the number of wells with cytopathic effect (CPE) based on the Reed and Muench method after incubating for 72 h.
Hematological tests
To determine the complete blood cell counts in mice, Blood was collected through the orbital venous plexus, blood was collected through the orbital venous plexus into anticoagulant tubes containing EDTAK2. Blood cell count was carried out under the automatic blood analyzer (BC-5000vet, Mindray, China), and the operation procedures were carried out according to the instrument operation instructions.
Histopathology and immunohistochemistry
Mice were euthanized at the appointed time. Lungs, spleens, kidneys and livers were collected. Tissue samples from mice were fixed in 10% buffered formalin and embedded in paraffin wax. Lung tissue was stained with hematoxylin and eosin (H&E) for histopathological assessment. For immunohistochemistry (IHC), quenching of paraffin-embedded lung tissues was performed with 3% hydrogen peroxide in methanol for 10 min. The tissue was added with 500-fold dilution of mouse monoclonal anti-SARS-CoV-2 N antibody at 4°C for overnight. After incubation with primary antibody, sections were washed with PBS for 3 times and species-matched secondary antibodies were applied for 2 h at 4°C. The sections were observed under ZEISS Axio Vert. A1 microscope. The details of the histopathology and immunohistochemistry assays can be found in our previous study (69).
Hamster challenge study
Hamsters were randomized into six groups of 9. Each Syrian hamster was given 250 µg of antibody at a dose of 10 mg/kg. In the prevention group, hamsters were administered purified equine IgG or F(ab’)2 by IP injection and challenged with 1,000 TCID50 wild-type SARS-CoV-2 Wuhan01 by the intranasal route in a final volume of 100 μL after 24 h. In the treatment group, hamsters were first challenged with 1,000 TCID50 SARS-CoV-2 Wuhan01 and given purified equine IgG or F(ab’)2 by IP injection 24 hours post-infection. All hamsters were monitored for body weight loss. At 3 dpi, four hamsters in each group were euthanized. Lungs and turbinates were harvested to detect the viral loads, histopathology and immunohistochemistry test.
Ethical statement
All procedures involving cells and animals were conducted in a biosafety level 3 laboratory (BSL-3) and approved by the animal experimental committee of the Laboratory Animal Center, Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences (approval number: IACUC of AMMS-11-2020-020). Animals are given enough drinking water, adequate food and comfortable living conditions in accordance with animal welfare principles.
Statistical analysis
Figures were generated using GraphPad Prism 8.0 software. Differences between means were evaluated using one-way ANOVA or two-way ANOVA and were deemed significant at P values of 0.05 or less.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The animal study was reviewed and approved by Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences.
Author contributions
Conceptualization, YoZ, BL and XX. Methodology, FY. Software, SWe, YiZ. Formal analysis, SWa, QH, BL. Investigation, NF, TW, JL. Resources, JH. Data curation, FY, JB, QH and EL. Writing original draft preparation, FY, JB and EL, SWa. Writing review and editing, FY. Supervision, YoZ, XX and SY. Project administration, JW, XX. Funding acquisition, YoZ. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Science and Technology Project (2020YFC0846100).
Acknowledgments
All authors would like to thank Chengfeng Qin for the gift of SARS-CoV-2 Wuhan01 and thank Dr. Jun Liu for assistance with BSL-3 animal experiments.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1066730/full#supplementary-material
References
1. Hu B, Guo H, Zhou P, Shi ZL. Characteristics of SARS-CoV-2 and COVID-19. Nat Rev Microbiol (2021) 19:141–54. doi: 10.1038/s41579-020-00459-7
2. Jiang S, Zhang X, Du L. Therapeutic antibodies and fusion inhibitors targeting the spike protein of SARS-CoV-2. Expert Opin Ther Targets (2021) 25:415–21. doi: 10.1080/14728222.2020.1820482
3. Walls AC, Park YJ, Tortorici MA, Wall A, McGuire AT, Veesler D. Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell (2020) 181:281–292.e6. doi: 10.1016/j.cell.2020.02.058
4. Yan R, Zhang Y, Li Y, Ye F, Guo Y, Xia L, et al. Structural basis for the different states of the spike protein of SARS-CoV-2 in complex with ACE2. Cell Res (2021) 31:717–9. doi: 10.1038/s41422-021-00490-0
5. Jackson CB, Farzan M, Chen B, Choe H. Mechanisms of SARS-CoV-2 entry into cells. Nat Rev Mol Cell Bio (2022) 23:3–20. doi: 10.1038/s41580-021-00418-x
6. Harrison AG, Lin T, Wang P. Mechanisms of SARS-CoV-2 transmission and pathogenesis. Trends Immunol (2020) 41:1100–15. doi: 10.1016/j.it.2020.10.004
7. Shah P, Canziani GA, Carter EP, Chaiken I. The case for S2: The potential benefits of the S2 subunit of the SARS-CoV-2 spike protein as an immunogen in fighting the COVID-19 pandemic. Front Immunol (2021) 12:637651. doi: 10.3389/fimmu.2021.637651
8. Amanat F, Thapa M, Lei T, Ahmed SMS, Adelsberg DC, Carreno JM, et al. SARS-CoV-2 mRNA vaccination induces functionally diverse antibodies to NTD, RBD, and S2. Cell (2021) 184:3936–48.e10. doi: 10.1016/j.cell.2021.06.005
9. Abdalla M, El-Arabey AA, Jiang X. Progress in research on the s protein as the target of COVID-19 vaccines. Expert Rev Vaccines (2021) 20:769–72. doi: 10.1080/14760584.2021.1918003
10. Hou YJ, Chiba S, Halfmann P, Ehre C, Kuroda M, Dinnon KH 3rd, et al. SARS-CoV-2 D614G variant exhibits efficient replication ex vivo and transmission in vivo. Science (2020) 370:1464–8. doi: 10.1126/science.abe8499
11. Janik E, Niemcewicz M, Podogrocki M, Majsterek I, Bijak M. The emerging concern and interest SARS-CoV-2 variants. Pathogens (2021) 10:633. doi: 10.3390/pathogens10060633
12. He X, Hong W, Pan X, Lu G, Wei X. SARS-CoV-2 omicron variant: Characteristics and prevention. MedComm (2020) (2021) 2(4):838–845. doi: 10.1002/mco2.110
13. VanBlargan LA, Adams LJ, Liu Z, Chen RE, Gilchuk P, Raju S, et al. A potently neutralizing SARS-CoV-2 antibody inhibits variants of concern by utilizing unique binding residues in a highly conserved epitope. Immunity (2021) 54:2399–2416.e6. doi: 10.1016/j.immuni.2021.08.016
14. Chi X, Yan R, Zhang J, Zhang G, Zhang Y, Hao M, et al. A neutralizing human antibody binds to the n-terminal domain of the spike protein of SARS-CoV-2. Science (2020) 369:650–5. doi: 10.1126/science.abc6952
15. Weinreich DM, Sivapalasingam S, Norton T, Ali S, Gao H, Bhore R, et al. REGN-COV2, a neutralizing antibody cocktail, in outpatients with covid-19. N Engl J Med (2021) 384:238–51. doi: 10.1056/NEJMoa2035002
16. Baum A, Fulton BO, Wloga E, Copin R, Pascal KE, Russo V, et al. Antibody cocktail to SARS-CoV-2 spike protein prevents rapid mutational escape seen with individual antibodies. Science (2020) 369:1014–8. doi: 10.1126/science.abd0831
17. Koenig PA, Das H, Liu H, Kummerer BM, Gohr FN, Jenster LM, et al. Structure-guided multivalent nanobodies block SARS-CoV-2 infection and suppress mutational escape. Science (2021) 371(6530):eabe6230. doi: 10.1126/science.abe6230
18. Xiang Y, Nambulli S, Xiao Z, Liu H, Sang Z, Duprex WP, et al. Versatile and multivalent nanobodies efficiently neutralize SARS-CoV-2. Science (2020) 370:1479–84. doi: 10.1126/science.abe4747
19. Huo J, Bas AL, Ruza RR, Duyvesteyn HME, Mikolajek H, Malinauskas T, et al. Author correction: Neutralizing nanobodies bind SARS-CoV-2 spike RBD and block interaction with ACE2. Nat Struct Mol Biol (2021) 28:326. doi: 10.1038/s41594-021-00566-w
20. Xu J, Xu K, Jung S, Conte A, Lieberman J, Muecksch F, et al. Nanobodies from camelid mice and llamas neutralize SARS-CoV-2 variants. Nature (2021) 595:278–82. doi: 10.1038/s41586-021-03676-z
21. Du Y, Shi R, Zhang Y, Duan X, Li L, Zhang J, et al. A broadly neutralizing humanized ACE2-targeting antibody against SARS-CoV-2 variants. Nat Commun (2021) 12:5000. doi: 10.1038/s41467-021-25331-x
22. Fedry J, Hurdiss DL, Wang C, Li W, Obal G, Drulyte I, et al. Structural insights into the cross-neutralization of SARS-CoV and SARS-CoV-2 by the human monoclonal antibody 47D11. Sci Adv (2021) 7(23):eabf5632. doi: 10.1126/sciadv.abf5632
23. Baig AM. Targeting neuroinvasion by SARS-CoV-2: Emerging trends in drug and antibody delivery to combat COVID-19. ACS Chem Neurosci (2021) 12:2555–7. doi: 10.1021/acschemneuro.1c00372
24. Smail SW, Saeed M, Twana A, Khudhur ZO, Younus DA, Rajab MF, et al. Inflammation, immunity and potential target therapy of SARS-COV-2: A total scale analysis review. Food Chem Toxicol (2021) 150:112087. doi: 10.1016/j.fct.2021.112087
25. Harvala H, Mehew J, Robb ML, Ijaz S, Dicks S, Patel M, et al. Convalescent plasma treatment for SARS-CoV-2 infection: analysis of the first 436 donors in England, 22 April to 12 may 2020. Euro Surveill (2020) 25(28):2001260. doi: 10.2807/1560-7917.ES.2020.25.28.2001260
26. Salazar E, Kuchipudi SV, Christensen PA, Eagar T, Yi X, Zhao P, et al. Convalescent plasma anti-SARS-CoV-2 spike protein ectodomain and receptor-binding domain IgG correlate with virus neutralization. J Clin Invest (2020) 130:6728–38. doi: 10.1172/JCI141206
27. Boonyaratanakornkit J, Morishima C, Selke S, Zamora D, McGuffin S, Shapiro AE, et al. Clinical, laboratory, and temporal predictors of neutralizing antibodies against SARS-CoV-2 among COVID-19 convalescent plasma donor candidates. J Clin Invest (2021) 131(3):e144930. doi: 10.1172/JCI144930
28. Casadevall A, Joyner MJ, Pirofski LA. SARS-CoV-2 viral load and antibody responses: the case for convalescent plasma therapy. J Clin Invest (2020) 130:5112–4. doi: 10.1172/JCI139760
29. Nguyen D, Simmonds P, Steenhuis M, Wouters E, Desmecht D, Garigliany M, et al. SARS-CoV-2 neutralising antibody testing in Europe: towards harmonisation of neutralising antibody titres for better use of convalescent plasma and comparability of trial data. Euro Surveill (2021) 26(27):2100568. doi: 10.2807/1560-7917.ES.2021.26.27.2100568
30. Hayes C, Rubenstein W, Gibb D, Klapper E, Tanaka J, Pepkowitz S. Blood group O convalescent plasma donations have significantly lower levels of SARS-CoV-2 IgG antibodies compared to blood group a donations. Transfusion (2021) 61:2245–9. doi: 10.1111/trf.16524
31. Zeng F, Dai C, Cai P, Wang J, Xu L, Li J, et al. A comparison study of SARS-CoV-2 IgG antibody between male and female COVID-19 patients: A possible reason underlying different outcome between sex. J Med Virol (2020) 92:2050–4. doi: 10.1002/jmv.25989
32. Li C, Yu D, Wu X, Liang H, Zhou Z, Xie Y, et al. Twelve-month specific IgG response to SARS-CoV-2 receptor-binding domain among COVID-19 convalescent plasma donors in wuhan. Nat Commun (2021) 12:4144. doi: 10.1038/s41467-021-24230-5
33. Zimmerli A, Monti M, Fenwick C, Eckerle I, Beigelman-Aubry C, Pellaton C, et al. Case report: Stepwise anti-inflammatory and anti-SARS-CoV-2 effects following convalescent plasma therapy with full clinical recovery. Front Immunol (2021) 12:613502. doi: 10.3389/fimmu.2021.613502
34. Cui Z, Li D, Yi S, Guo Y, Dong G, Niu J, et al. Equine immunoglobulin F(ab')2 fragments protect cats against feline calicivirus infection. Int Immunopharmacol (2019) 75:105714. doi: 10.1016/j.intimp.2019.105714
35. Cunha LER, Stolet AA, Strauch MA, Pereira VAR, Dumard CH, Gomes AMO, et al. Polyclonal F(ab')2 fragments of equine antibodies raised against the spike protein neutralize SARS-CoV-2 variants with high potency. iScience (2021), 24(11):103315. doi: 10.1016/j.isci.2021.103315
36. Wilson BZ, Larsen J, Smelski G, Dudley S, Shirazi FM. Use of crotalidae equine immune F(ab')2 antivenom for treatment of an agkistrodon envenomation. Clin Toxicol (Phila) (2021) 59:1023–6. doi: 10.1080/15563650.2021.1892718
37. Racine T, Denizot M, Pannetier D, Nguyen L, Pasquier A, Raoul H, et al. In vitro characterization and In vivo effectiveness of Ebola virus specific equine polyclonal F(ab')2. J Infect Dis (2019) 220:41–5. doi: 10.1093/infdis/jiz068
38. Ainsworth S MS, Pleass RJ. Animal derived antibodies should be considered alongside convalescent human plasma to deliver treatments for COVID-19. Wellcome Open Res (2020) 5:115–6. doi: 10.12688/wellcomeopenres.15990.1
39. Li E, Chi H, Huang P, Yan F, Zhang Y, Liu C, et al. A novel bacterium-like particle vaccine displaying the MERS-CoV receptor-binding domain induces specific mucosal and systemic immune responses in mice. Viruses (2019) 11(9):799. doi: 10.3390/v11090799
40. Sia SF, Yan LM, Chin AWH, Fung K, Choy KT, Wong AYL, et al. Pathogenesis and transmission of SARS-CoV-2 in golden hamsters. Nature (2020) 583:834–8. doi: 10.1038/s41586-020-2342-5
41. Chen RE, Winkler ES, Case JB, Aziati ID, Bricker TL, Joshi A, et al. In vivo monoclonal antibody efficacy against SARS-CoV-2 variant strains. Nature (2021) 596:103–8. doi: 10.1038/s41586-021-03720-y
42. Planas D, Veyer D, Baidaliuk A, Staropoli I, Guivel-Benhassine F, Rajah MM, et al. Reduced sensitivity of SARS-CoV-2 variant delta to antibody neutralization. Nature (2021) 596:276–80. doi: 10.1038/s41586-021-03777-9
43. H. The Lancet. The resurgence of convalescent plasma therapy. Lancet Haematol (2020) 7:e353. doi: 10.1016/S2352-3026(20)30117-4
44. Zylberman V, Sanguineti S, Pontoriero AV, Higa SV, Cerutti ML, Morrone Seijo SM, et al. Development of a hyperimmune equine serum therapy for COVID-19 in Argentina. Medicina (B Aires) (2020) 80 Suppl 3:1–6.
45. Chippaux JP, Massougbodji A, Stock RP, Alagon A, Investigators of African Antivipmyn in Benin. Clinical trial of an F(ab')2 polyvalent equine antivenom for African snake bites in Benin. Am J Trop Med Hyg (2007) 77:538–46. doi: 10.4269/ajtmh.2007.77.538
46. Quiambao BP, Dytioco HZ, Dizon RM, Crisostomo ME, Laot TM, Teuwen DE. Rabies post-exposure prophylaxis in the Philippines: health status of patients having received purified equine F(ab')(2) fragment rabies immunoglobulin (Favirab). PloS Negl Trop Dis (2008) 2:e243. doi: 10.1371/journal.pntd.0000243
47. Luczkowiak J, Radreau P, Nguyen L, Labiod N, Lasala F, Veas F, et al. Potent neutralizing activity of polyclonal equine antibodies against severe acute respiratory syndrome coronavirus 2 variants of concern. J Infect Dis (2022) 227:35–9. doi: 10.1093/infdis/jiac331
48. Moreira-Soto A, Arguedas M, Brenes H, Bujan W, Corrales-Aguilar E, Diaz C, et al. High efficacy of therapeutic equine hyperimmune antibodies against SARS-CoV-2 variants of concern. Front Med (Lausanne) (2021) 8:735853. doi: 10.3389/fmed.2021.735853
49. Nie J, Li Q, Wu J, Zhao C, Hao H, Liu H, et al. Establishment and validation of a pseudovirus neutralization assay for SARS-CoV-2. Emerg Microbes Infect (2020) 9:680–6. doi: 10.1080/22221751.2020.1743767
50. Xiong HL, Wu YT, Cao JL, Yang R, Liu YX, Ma J, et al. Robust neutralization assay based on SARS-CoV-2 s-protein-bearing vesicular stomatitis virus (VSV) pseudovirus and ACE2-overexpressing BHK21 cells. Emerg Microbes Infect (2020) 9:2105–13. doi: 10.1080/22221751.2020.1815589
51. Hoffmann M, Arora P, Gross R, Seidel A, Hornich BF, Hahn AS, et al. SARS-CoV-2 variants B.1.351 and P.1 escape from neutralizing antibodies. Cell (2021) 184:2384–2393.e12. doi: 10.1016/j.cell.2021.03.036
52. Planas D, Bruel T, Grzelak L, Guivel-Benhassine F, Staropoli I, Porrot F, et al. Sensitivity of infectious SARS-CoV-2 B.1.1.7 and B.1.351 variants to neutralizing antibodies. Nat Med (2021) 27:917–24. doi: 10.1038/s41591-021-01318-5
53. Zhou D, Dejnirattisai W, Supasa P, Liu C, Mentzer AJ, Ginn HM, et al. Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera. Cell (2021) 184:2348–61.e6. doi: 10.1016/j.cell.2021.02.037
54. Yu X, Wei D, Xu W, Li Y, Li X, Zhang X, et al. Reduced sensitivity of SARS-CoV-2 omicron variant to antibody neutralization elicited by booster vaccination. Cell Discovery (2022) 8:4. doi: 10.1038/s41421-022-00375-5
55. Lazarevic I, Pravica V, Miljanovic D, Cupic M. Immune evasion of SARS-CoV-2 emerging variants: What have we learnt so far? Viruses (2021) 13(7):1192. doi: 10.3390/v13071192
56. Garcia-Beltran WF, Lam EC, St Denis K, Nitido AD, Garcia ZH, Hauser BM, et al. Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity. Cell (2021) 184:2523. doi: 10.1016/j.cell.2021.04.006
57. Prevost J, Finzi A. The great escape? SARS-CoV-2 variants evading neutralizing responses. Cell Host Microbe (2021) 29:322–4. doi: 10.1016/j.chom.2021.02.010
58. Chen RE, Zhang X, Case JB, Winkler ES, Liu Y, VanBlargan LA, et al. Resistance of SARS-CoV-2 variants to neutralization by monoclonal and serum-derived polyclonal antibodies. Nat Med (2021) 27:717–26. doi: 10.1038/s41591-021-01294-w
59. Dups J, Middleton D, Long F, Arkinstall R, Marsh GA, Wang LF. Subclinical infection without encephalitis in mice following intranasal exposure to nipah virus-Malaysia and nipah virus-Bangladesh. Virol J (2014) 11:102. doi: 10.1186/1743-422X-11-102
60. Lyon GM, Mehta AK, Varkey JB, Brantly K, Plyler L, McElroy AK, et al. Clinical care of two patients with Ebola virus disease in the united states. N Engl J Med (2014) 371:2402–9. doi: 10.1056/NEJMoa1409838
61. Bournazos S, Gupta A, Ravetch JV. The role of IgG fc receptors in antibody-dependent enhancement. Nat Rev Immunol (2020) 20:633–43. doi: 10.1038/s41577-020-00410-0
62. Gao R, Sheng Z, Sreenivasan CC, Wang D, Li F. Influenza a virus antibodies with antibody-dependent cellular cytotoxicity function. Viruses (2020) 12(3):276. doi: 10.3390/v12030276
63. Winkler ES, Gilchuk P, Yu J, Bailey AL, Chen RE, Chong Z, et al. Human neutralizing antibodies against SARS-CoV-2 require intact fc effector functions for optimal therapeutic protection. Cell (2021) 184:1804–20.e16. doi: 10.1016/j.cell.2021.02.026
64. Ilinykh PA, Huang K, Santos RI, Gilchuk P, Gunn BM, Karim MM, et al. Non-neutralizing antibodies from a marburg infection survivor mediate protection by fc-effector functions and by enhancing efficacy of other antibodies. Cell Host Microbe (2020) 27:976–991.e11. doi: 10.1016/j.chom.2020.03.025
65. Howell KA, Brannan JM, Bryan C, McNeal A, Davidson E, Turner HL, et al. Cooperativity enables non-neutralizing antibodies to neutralize ebolavirus. Cell Rep (2017) 19:413–24. doi: 10.1016/j.celrep.2017.03.049
66. Abreu-Mota T, Hagen KR, Cooper K, Jahrling PB, Tan G, Wirblich C, et al. Non-neutralizing antibodies elicited by recombinant lassa-rabies vaccine are critical for protection against lassa fever. Nat Commun (2018) 9:4223. doi: 10.1038/s41467-018-06741-w
67. Casadevall A, Pirofski LA, Joyner MJ. The principles of antibody therapy for infectious diseases with relevance for COVID-19. mBio (2021) 12(2):e03372–20. doi: 10.1128/mBio.03372-20
68. Yamin R, Jones AT, Hoffmann HH, Schafer A, Kao KS, Francis RL, et al. Fc-engineered antibody therapeutics with improved anti-SARS-CoV-2 efficacy. Nature (2021) 599:465–70. doi: 10.1038/s41586-021-04017-w
Keywords: SARS-CoV-2, equine antibody, variants of concern, variants of interest, broad-spectrum neutralizing activity
Citation: Li E, Han Q, Bi J, Wei S, Wang S, Zhang Y, Liu J, Feng N, Wang T, Wu J, Yang S, Zhao Y, Liu B, Yan F and Xia X (2023) Therapeutic equine hyperimmune antibodies with high and broad-spectrum neutralizing activity protect rodents against SARS-CoV-2 infection. Front. Immunol. 14:1066730. doi: 10.3389/fimmu.2023.1066730
Received: 11 October 2022; Accepted: 06 February 2023;
Published: 17 February 2023.
Edited by:
Pragya Dhruv Yadav, ICMR-National Institute of Virology, IndiaReviewed by:
James Tristan Gordy, Bloomberg School of Public Health, Johns Hopkins University, United StatesGiorgia Zanetti, Columbia University Irving Medical Center, United States
Copyright © 2023 Li, Han, Bi, Wei, Wang, Zhang, Liu, Feng, Wang, Wu, Yang, Zhao, Liu, Yan and Xia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yongkun Zhao, zhaoyongkun1976@126.com; Bo Liu, llbb520@aliyun.com; Feihu Yan, yanfh1990@163.com; Xianzhu Xia, xiaxzh@cae.cn
†These authors have contributed equally to this work