Impact Factor 4.076

The 3rd most cited journal in Microbiology

Original Research ARTICLE

Front. Microbiol., 11 July 2017 | https://doi.org/10.3389/fmicb.2017.01308

Amino Acid Substitution K470R in the Nucleoprotein Increases the Virulence of H5N1 Influenza A Virus in Mammals

Lin Chen1,2, Chengmin Wang1, Jing Luo1, Meng Li1,2, Huimin Liu1,2, Na Zhao1,2, Jingjing Huang1, Xili Zhu3, Guoyao Ma1,2, Guohui Yuan1 and Hongxuan He1*
  • 1National Research Center for Wildlife-Borne Diseases, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
  • 2University of the Chinese Academy of Sciences, Beijing, China
  • 3State Key Laboratory of Stem Cell and Reproductive Biology, Institute of Zoology, Chinese Academy of Science, Beijing, China

H5N1 is a highly pathogenic influenza A virus (IAV) and poses a major threat to the public health. The nucleoprotein (NP) has a multiple functions during the viral life cycle, however, the precise role of NP mutants in viral replication and pathogenicity is not completely understood. Here, we attempted to identify five residues in NP that may contribute to viral replication or pathogenicity. Of these, K227R, K229R, and K470R viruses were successfully rescued by reverse genetic, but the K91R and K198R viruses were not viable. A mini-genome assay demonstrated that the NP mutations K91R and K198R significantly decreased the polymerase activity. Moreover, these two mutations resulted in disrupted cellular localization in mammalian cells. Importantly, mutation at position 470 of NP significantly increased its virulence in vitro and in vivo. These findings demonstrated that the NP protein plays a major role in influenza virulence and pathogenicity, which adds to the knowledge of IAV virulence determinants and may benefit IAV surveillance.

Introduction

H5N1, a highly pathogenic avian influenza A (HPAI) virus, is a serious public health threat (El-Shesheny et al., 2017; Ly et al., 2017). Up to 5 April 2017, more than 856 humans were infected the H5N1 virus, with at least 452 deaths1. Influenza A viruses (IAVs) are enveloped viruses consisting of a negative-sense RNA genome composed of eight segments (Murphy and Webster, 1996), and encode up to 17 proteins: polymerase basic protein 1(PB1), PB1-F2 (Chen et al., 2001), PB1-N40, polymerase basic protein 2 (PB2), polymerase acidic protein (PA), PA-X (Jagger et al., 2012), PA-N155, PA-N182 (Muramoto et al., 2013), hemagglutinin (HA), nucleoprotein (NP), neuraminidase (NA), nonstructural protein1 (NS1), NS2, NS3 (Selman et al., 2012), matrix protein 1 (M1), M2 and M42 (Wise et al., 2012). The IAV particle contains three major components: the envelope, viral ribonucleoprotein (vRNP) core, and M1. The vRNP core consists of NP and the viral RNA (vRNA) polymerase that includes PB1, PB2, and PA, which together constitute the functional templates for replication and transcription of the viral genome (Honda et al., 1988).

Multiple amino acids residues in proteins of influenza A virus have been identified as the critical determinants for virulence in mammalian hosts. Amino acids 627K and 701D in the PB2 protein increased the polymerase activity and the virulence of the influenza virus in mammalian hosts (Tarendeau et al., 2008; Bussey et al., 2010). PB1 residues 473V and 598P increased the polymerase and virulence in mammalian cells (Xu et al., 2012). PA mutations P103H and S659L reduced viral replication in mammalian cells and attenuated pathogenicity in mice (Zhao et al., 2016). In addition to the polymerase complex, the HA is another major determinant of the host range of influenza A virus (Bottcher-Friebertshauser et al., 2014), and several basic amino acids at the HA cleavage sites contribute to high pathogenicity and induce systematic spread of the viruses in mammals (Chen et al., 1998). Moreover, T at position 215(T215) of the NS1 protein of the H3N2 virus is important for viral replication (Hale et al., 2009). All recombinant viruses containing M1 mutants (Y132A/F/D) were unable to be rescued (Wang et al., 2013). These studies suggested that viral pathogenicity was affected by amino acid mutations in viral proteins.

The NP is highly conserved among influenza A viruses. This protein is a major structural component of the viral particle and encapsulates the viral RNA (Portela and Digard, 2002). It contains a RNA-binding region at its N terminus and two domains, which are responsible for NP-NP self-interaction, at residues 189 to 358 and 371 to 465 (Ng et al., 2009). NP is a shuttle protein between the nucleus and cytoplasm (Neumann et al., 1997). In the early phases of infection, NP accumulates in the nucleus. During viral assembly and maturation, NP is exclusively distributed in cytoplasm. Several studies have reported that NP protein is associated with the virulence in mammalian cells. To assess the conserved amino acids, (Li et al., 2009) used a comprehensive analysis to identify 15 amino acid residues that are crucial for viral-genome replication, and 17 amino acid substitutions had no or little effect on viral growth in vitro. Several post-translational modifications affected the NP functions. K184 is an ubiquitination site on NP and this residue was crucial for virus RNA replication (Liao et al., 2010). K4 and K7 of the NP protein are the SUMO acceptor sites, and the mutant WSN-NP K4,7R virus was highly attenuated and had aberrant nuclear/cytoplasmic trafficking dynamics (Han et al., 2014). Phosphorylation of the NP protein was first observed many years ago (Kistner et al., 1989). It was reported that S9, Y10, S165, and Y296 were the phosphorylation sites. S9 and Y10 controlled nuclear import of NP, and mutation of these two sites dramatically decreased viral titers in cell culture and modulation of the viral polymerase activity in a mouse model (Zheng et al., 2015). Although many efforts have already been made to understand the role of NP residues, the precise role of some NP mutants in viral replication and pathogenicity remains unclear.

In our previous studies, five residues (lysine 91, 198, 227, 229, and 470) of NP were shown to be modified by acetylation using mass spectrometry (unpublished data). In the current study, we focused on the five residues of NP in viral replication and pathogenicity. Using reverse genetics, we successfully rescued K227R, K229R, and K470R mutant viruses and analyze the virulence of the recombinant viruses in vitro and in vivo. For the K91R and K198R mutations that did not support viral growth, we studied the cellular distribution in mammalian cells.

Materials and Methods

Cells and Viruses

MDCK cells, 293T cells and A549 cells were cultured in Dulbecco’s minimal essential medium (DMEM) supplemented with 10% fetal bovine serum (FBS), 100 IU/ml penicillin and 100 μg/ml streptomycin. The A/Qing Hai/environment/2005 (H5N1) strain was propagated in 10-day-old embryonated chicken eggs.

Biosafety and Ethics Statement

All experiments involving H5N1 viruses were performed in a Biosafety Level 3(BSL-3) containment laboratory in the Research Center for Wildlife Diseases, which was approved by the Chinese Academy of Science. All animal experiments were performed in compliance with the Guide for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of the People’s Republic of China. The protocol was approved by the Committee on the Ethics of Animal Experiments of the Institute of Zoology, Chinese Academy of Sciences (approval number: IOZ-15042). Mice were euthanized with CO2 to minimize animal suffering if the weight loss of more than 20% of the body weight was observed.

Plasmids Construction

The full-length NP (A/Qing Hai/environment/2005(H5N1)) was cloned into the pcDNA3.1-Myc and pHW2000 vectors. Lysine (K) to arginine (R) mutations in the NP genes were introduced using the QuikChange site-directed mutagenesis kit (Stratagene) with primers designed following the manufacturer’s protocol. All primers used in this study are shown in Table 1. All plasmid constructs were sequenced to confirm the sequence integrity.

TABLE 1
www.frontiersin.org

TABLE 1. Primers used for site-directed mutagenesis.

Generation of Recombinant Influenza Viruses by Reverse Genetics

Recombinant viruses were rescued by eight-plasmid reverse genetics followed by two-round plaque purification and propagation on MDCK cells as previously described (Hoffmann et al., 2000). Briefly, the NP gene from A/environment/Qinghai/1/2008 (H5N1) (as wild type) or the NP gene with the mutation K91R, K198R, K227R, K229R, or K470R was co-transfected with the plasmids encoding the other seven genes from the A/environment/Qinghai/1/2008 (H5N1) virus into 293T cells and MDCK cells. The medium was replaced 24 h later with DMEM plus 1 μg/ml tosylsulfonyl phenylalanyl chloromethyl ketone (TPCK)-treated trypsin. At 72 h after transfection, the supernatant containing the recombinant viruses was harvested and then centrifuged at 2,000 × g for 10 min to remove the cell debris.

Plaque Assay

Monolayer cultures of MDCK cells in 6-well plates were washed twice with phosphate-buffered saline (PBS). The cells were infected with serial 10-fold dilutions of virus for 1 h 37°C and shaken at every 15 min. The virus inoculums were then overlaid with agar overlay medium (DMEM supplemented with 0.6% low-melting-point agarose and 1 μg/ml TPCK-treated trypsin) and incubated at 37°C. After 3 days, the cells were fixed with 10% formaldehyde and stained with 0.5% crystal violet solution.

Viral Replication Kinetics in Mammalian Cells

A549 cells were inoculated with the indicated virus at a multiplicity of infection (MOI) of 1/0.1 plaque-forming units per cell. After 1 h of incubation at 37°C, cells were washed twice with 1× PBS and were maintained in fresh medium for virus replication. The supernatants were collected at 3, 6, 9, 12, 24, 36, and 48 h after virus infection. Virus titers in the supernatants were detected by TCID50.

Mini-genome Assay

The 293T cells in 6-well plates were transfected with plasmid (PB2/pHW2000, PB1/ pHW2000, PA/ pHW2000, and NP/pHW2000 or its mutants) DNA (1 μg each) and a luciferase RNA expression (vNS-luc/pHH21) vector and pRL-TK(50 ng), which expresses Renilla luciferase as an internal control. The cells were lysed at 36 h post-transfection, and luciferase activity was measured using the Dual Luciferase Assay System (Promega) according to the manufacturer’s protocol. Polymerase activity is calculated by normalizing the firefly luciferase activity to the Renilla luciferase activity. Polymerase activity of wild-type was set to 100%.

Confocal Microscopy

First, 293T cells were seeded on coverslips and grown to 80–90% confluence. Cells were then cotransfected with 1 μg of NP WT or mutants plasmids. At 9, 24, and 30 h post-transfection, the coverslips carrying cells were washed with PBS and then fixed with 4% paraformaldehyde and permeabilized with 0.1% Triton X-100. Cells were then blocked with 1% bovine serum albumin (BSA) and stained with NP antibodies for 2 h at RT, washed 3 times with PBS, and incubated with an Alexa Fluor 488-conjugated goat secondary antibody against rabbit IgG (1:500 in PBS) for 1 h. The cells were then incubated with PBS containing DAPI (4, 6-diamidino-2-phenylindole) for 5 min and then washed three times with PBS. The excitation and emission wavelengths were set at 488 nm and 519 nm, respectively. Prepared samples were observed under a confocal laser scanning microscope (Leica 780).

RNA Isolation, Reverse Transcription, and Viral mRNA, cRNA and vRNA Quantification

Total RNA from infected cells was extracted using TRIzol (Invitrogen) according to the manufacturer’s instructions. The cDNAs of the corresponding virus mRNA, cRNA and vRNA were first reverse-transcribed by using the GoScript Reverse Transcription System (Promega) with the tagged primers shown in Table 2. The analysis of relative vRNA, cRNA, and mRNA expression was performed using an ABI 7500 machine. The cycling conditions comprised an initial denaturation step of 3 min at 95°C, followed by 40 two-step cycles (95°C for 5 s and 60°C for 60 s). Dissociation curve analysis was performed after each assay to ensure specific target detection. The cellular β-actin mRNA in the infected cells was used as an internal control.

TABLE 2
www.frontiersin.org

TABLE 2. Primers used for reverse transcription and Real-time PCR.

Western Blot

Cells were washed once with PBS and lysed in lysis buffer (50 mM Tris-HCl pH 7.6, 150 mM NaCl, 0.1% SDS, 1% Triton-100) further supplemented with protease inhibitor (Roche). The whole cell lysates were resolved by SDS-PAGE. Proteins were transferred to nitrocellulose membrane and blocked with 5% nonfat milk in TBST. The membrane was probed with primary antibody followed by secondary antibody at room temperature for 2–4 h. The membrane was washed at least three times with TBST (0.1% Tween 20) after incubation with each antibody and developed with the Dura chemiluminescent kit (Millipore).

Animals

Female 4- to 6-week-old BALB/c mice (Stanford, China) were lightly anesthetized with CO2 and inoculated intranasally with a 50 μl volume of the indicated dose. The mice in the control group were inoculated with PBS. Mice were euthanized on day 3 post inoculation, and lungs were collected and titrated for viral titer in MDCK cells and hematoxylin and eosin (HE) staining. The remaining mice (n = 4) were monitored daily for 14 days and assessed for body weight and mortality.

Statistical Analysis

Data are expressed as the mean ± SDs. Student’s t-test was used for statistical comparisons; p < 0.05 was considered significant.

Results

Generation of NP K227R, K229R, and K470R Mutant Viruses

To determine the function of five lysines (K91, K198, K227, K229, and K470) in NP, we performed site-directed mutagenesis to replace the lysine (K) residues with arginine (R) and generated the K91R, K198R, K227R, K229R, and K470R recombinant viruses using plasmid-driven reverse genetics. MDCK and 293T cells were co-transfected with NP mutants and the other seven plasmids, and at 72 h post transfection, viral proteins within transfected cells were detected by western blot. The expression of viral protein (PB1, NA and PA) was almost similar in all of the wild-type (WT) and mutant NP-transfected cells, although the NP expression level in K91R, K198R and K227R was lower than that of wild-type (Figure 1A). The NP rWT, rK227R, rK229R, and rK470R viruses were successfully rescued. The plaque sizes of the rK227R, rK229R and rK470R mutant viruses were similar to that of the rWT virus (Figure 1B). The number of plaque in rK470R was higher than those of rWT. In contrast, the K91R and K198R mutant viruses were not viable in three independent experiments, suggesting the two mutations in NP were critical for the viral life cycle.

FIGURE 1
www.frontiersin.org

FIGURE 1. Rescue mutant viruses by reverse genetics. (A) MDCK and 293T cells were cotransfected with eight plasmids of the H5N1. At 72 h p.t., the transfected cells were lysed for western blot analysis. PA, NP, PB1, and NA were detected with respective antibodies. (B) The plaque morphology of NP-rWT, NP-rK227R, NP-rK229R, and NP-rK470R viruses.

Effect of NP Mutants on the Polymerase Activity

The polymerase complex of influenza A virus, composed of PA, PB1, PB2, and NP, is essential for viral transcription and replication (Honda et al., 1988). Minireplicon assays were performed to evaluate the impact of the NP substitutions K91R, K198R, K227R, K229R, and K470R on the polymerase activity in human embryonic kidney (293T) cells at 37°C. As shown in Figure 2, the K91R and K198R mutants displayed substantially lower luciferase activity than that of WT, while the K227R, K229R, and K470R mutants increased the polymerase activity by 1.4-, 1.5-, and 1.9-fold, respectively. Our results indicated that K91R and K198R mutants likely lost the ability to support viral polymerase activity, resulting in the growth defects of the mutant viruses.

FIGURE 2
www.frontiersin.org

FIGURE 2. Nucleoprotein (NP) mutants affected the influenza A virus genome RNA replication using the mini-replicon reporter system. HEK293T cells were transfected with pPolI–Luc and plasmids for the expression of the viral PB2, PB1, PA and NP (wild-type or mutants). Renilla luciferase was used as an internal control. Luciferase assays were performed at 24 h after transfection. The data were normalized relative to the values detected in the cells transfected with wild-type NP. Values are the mean ± SEM of three separate experiments. P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001.

K91 and K198 Disrupted the Cellular Distribution of NP

Nucleoprotein is a shuttle protein between the nucleus and the cytoplasm. Mutations of several residues, such as the R8A and Y296F mutations, cause cytoplasmic accumulation of NP (Zheng et al., 2015). We hypothesized that K91 and K198 would regulate the cellular distribution of NP. To test our hypothesis, imumunofluorescence assays (IFAs) were carried out to examine the cellular localization of NPs. Plasmids encoding NP WT, K91R or K198R were transfected into 293T cells, and the transfected cells were fixed at 9, 24, and 30 h p.t. Imaging revealed that NP WT exhibited nuclear localization at 9 and 24 h p.t. (Figures 3A,B). However, K91R and K198R mutants displayed cytoplasmic localization at 24 and 30 h p.t., respectively (Figures 3A,B). In addition, the results of cytoplasmic and nuclear fractionation experiments also demonstrated that K91R and K198R mutations increased the cytoplasmic distribution of NP (Figure 3C). We speculated that disrupting NP cellular distribution was associated with the lower RNP activity. Together, our results showed that NP-K91R/K198R disrupted the cellular location of NP.

FIGURE 3
www.frontiersin.org

FIGURE 3. K91R and K198R mutants disrupted the cell distribution of NP. (A) The localization of WT NP and its mutants (K91R/K198R) was determined using IFAs. The 293T cells were transfected with plasmids expressing mutant NPs or WT NP. At 9, 24, and 30 h p.t., the transfected cells were fixed and stained with anti-NP antibody (green). The nucleus was stained with DAPI (blue). Scale bar: 10 μm. (B) At least 200 cells in each group for three independent assays were scored as predominantly nuclear (N), nuclear and cytoplasmic (N+C), or predominantly cytoplasmic (C). P-values were determined using Student’s t-test by comparing to the rWT virus group. P < 0.05; ∗∗P < 0.01. (C) 293T cells were transfected with a plasmid expressing WT NP, K91R or K198R mutant and harvested at 30 h p.t. The cells were separated into the nuclear fraction (Nuc) and the cytoplasmic fraction (Cyto). Each fraction was examined by Western blotting.

Effects of K227R, K229R, and K470R on Virus Growth In Vitro

To evaluate the impact of the substitutions on virus growth, we inoculated recombinant viruses rK227R, rK229R, and rK470R into A549 cells at an MOI of 1 and harvested supernatants at 3, 6, 9, and 12 h post-infection (hpi). As shown in Figure 4A, viral titers of rK470R were 100-fold and 50-fold higher than those of rWT at 9 and 12 hpi, respectively. Viral titers of rK229R were lower than that of rWT at 12 hpi, while those of rK227R were similar with the rWT. Therefore, the K470R mutants enhanced virus growth of H5N1 in mammalian cells in a single cycle growth curve. For multiple cycle growth curves, A549 cells were inoculated at an MOI = 0.1, and the supernatants were harvested at 12, 24, 36, and 48 hpi. The rK227R and rWT viruses had the similar growth kinetics in A549 cells. In contrast, the titer of rK470R was significantly increased compared to that of the rWT virus at 36 and 48 hpi (20-fold and 16-fold, respectively), and the rK229R virus exhibited significantly reduced replication compared to the rWT virus (Figure 4B). Together, our results indicated that the K470R mutant increased the virus growth in mammalian cells, while the K229R suppressed the virus growth and the K227R had no effect.

FIGURE 4
www.frontiersin.org

FIGURE 4. The characteristics of the K227R, K229R and K470R mutant recombinant viruses in vitro. (A) A549 cells were infected with rK227R, rK229R, rK470R or rWT virus at an MOI of 1. The supernatants were collected at 3, 6, 9 and 12 h after infection and then titrated in MDCK cells by TCID50. (B) A549 cells were infected with rK227R, rK229R, rK470R or rWT virus at an MOI of 0.1. The supernatants were collected at 12, 24, 36 and 48 h after infection, and then titrated in MDCK cells by TCID50. Values are the mean ± SEM of three separate experiments. P-values were determined using Student’s t-test by comparing to the rWT virus group. P < 0.05; ∗∗P < 0.01.

Effects of K227R, K229R, and K470R on the Synthesis of Viral mRNA, cRNA, and vRNA

In influenza virus-infected cells, the vRNPs direct two types of RNA synthesis: mRNA synthesis (transcription) and vRNA amplification (replication). The vRNA template is copied to form full-length positive-stranded RNA (cRNA), which forms the template for synthesis of vRNA segments for amplification of mRNA synthesis and packaging into progeny virions (York and Fodor, 2013). We investigated whether these substitutions would affect viral vRNA, cRNA, and mRNA synthesis. A549 cells were infected with the rWT, rK227R, rK229R or rK470R virus at an MOI of 1 for 6 h, and the synthesis of viral NP mRNA, cRNA, and vRNA was quantified. The results showed that the rK227R, rK229R, or rK470R viruses produced higher levels of mRNA (Figure 5A), cRNA (Figure 5B), and vRNA (Figure 5C) than those of rWT. Among these substitutions, K470R produced the highest levels of mRNA, cRNA, and vRNA, which were 10.8-fold, 2.6-fold, and 5.1-fold greater than those of rWT, respectively (Figure 5). Taken together, these results indicated that K470R enhanced viral transcription and replication by increasing the synthesis of viral mRNA, cRNA, and vRNA in human cells.

FIGURE 5
www.frontiersin.org

FIGURE 5. Quantification of viral mRNA, cRNA and vRNA. A549 cells were infected with rK227R, rK229R, rK470R, or rWT virus at an MOI of 1. At 6 h after infection, total cellular RNA was analyzed by two-step real-time RT-qPCR to discriminate and quantify viral NP mRNA (A), cRNA (B) and vRNA (C). The reactions were carried out in triplicate and normalized to the levels of cellular β-actin mRNA. The value for rWT was set to 1. Data represent the mean ± SD values. P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001, compared with that of rWT by Student’s t-test.

Effects of K227R, K229R, and K470R on Viral Pathogenicity in Mice

To assess whether the phenotypes of mutant viruses in vitro manifested in vivo, we evaluated the growth and pathogenicity of these viruses in a mouse model. Seven BALB/c mice each were intranasally inoculated with 1 × 104 TCID50 of rWT or mutant rK227R, rK229R, or rK470R virus. Control mice were inoculated with PBS. Three mice from each group were euthanized on day 3 post-infection, and the remaining four mice were monitored daily for 14 days for changes in body weight and for death. Mice infected with the rWT virus rapidly lost weight, and two of four infected mice died by day 8. The weight change pattern and the survival rate of rK227R were similar to that of rWT-infected mice. However, only one of four mice infected with rK229R died by day 10. In contrast, four rK470R virus-infected mice died by day 8 (Figures 6A,B), consistent with the increased the virus growth of the K470R in mammalian cells. We next determined lung index and viral titers on day 3 post-infection. The lung index and the virus titer of rK470R-infected mice were higher than other groups (Figures 6C,D).

FIGURE 6
www.frontiersin.org

FIGURE 6. The characteristics of the K227R, K229R, and K470R mutant recombinant viruses in vivo. (A,B) BALB/c mice (eight/group) were inoculated intranasally with rWT, rK227R, rK229R, and rK470R viruses at a dose of 1 × 104 TCID50. Infected animals were monitored daily for body weight changes (A) and survival (B). (C,D) Three mice from each group were euthanized on day 3 postinfection to determine the lung index (C) and lung viral titer (D). The lung index was calculated as the lung wet weight/body weight × 100. Values are the mean ± SEM. of three separate experiments. P-values were determined using Student’s t-test by comparing to the rWT virus group. P < 0.05; ∗∗P < 0.01.

To further compare the pathogenicity of these viruses in vivo, we performed histological analysis of lungs from mice infected with the rWT, rK227R, rK229R, or rK470R virus on day 3 post-infection (Figure 7). Control group mice displayed thin alveolar walls and large air spaces, which indicated normal lung function (Figure 7A). Mice infected with rWT showed inflammation with immune infiltrates in the alveolar interstitial spaces and mild hemorrhage (Figure 7B). Similarly, infection with rK227R caused mild hemorrhage (Figure 7C). In contrast to the rWT, infection with rK229R exhibited mild inflammation (Figure 7D). In addition, infection with rK470R led to a significant increase of cellular infiltration and severe hemorrhage (Figure 7E). Collectively, these findings demonstrated the enhanced effect of rK470R on the pathogenicity of subtype H5N1 avian influenza in mice.

FIGURE 7
www.frontiersin.org

FIGURE 7. Histopathologic changes in the lungs of infected mice. The images shown are hematoxylin- and eosin-stained sections of lungs from BALB/c mice infected with PBS (A), rWT (B), rK227R (C), rK229R (D), or rK470R (E) virus. Original magnification, ×20.

Discussion

Potential mutations in influenza A virus that increase pathogenicity have not been thoroughly examined. In this study, we found that several lysine-to-arginine substitutions of the NP protein markedly affected viral virulence in vitro and in vivo. Compared with the rWT virus, the rK470R mutant virus showed enhanced virulence. The crystal structure of NP of H5N1 has been determined (Ng et al., 2008). The amino acid at position 470 was exposed on the surface of the trimeric NP, indicating the K470 residue had a critical role. The replacement of lysine with arginine could affect the interaction of NP with other viral proteins. Therefore, K470R enhanced the polymerase activity in minireplicon assays, and the vRNA, cRNA, and mRNA levels were higher than those of WT. By sequence alignment and analyses, we showed that the amino acid at position 470 of NP is polymorphic among influenza A viruses, but the 91, 198, 227, and 229 residues were conserved (Table 3). Our inspection of 13773 influenza A virus NP sequences (obtained from the Influenza Research Database2) revealed that most NP proteins encoded NP-470K. Interestingly, 92.6% of human H1N1 isolated before March 2009 presented NP-470R, whereas 99.7% of avian-origin H1N1 viruses harbored NP-470K (Table 4). In March 2009, a novel influenza A virus [influenza A (H1N1) 2009] emerged in Mexico and rapidly spread worldwide. The virus is more pathogenic than seasonal influenza viruses and caused more than 18,000 deaths (Scalera and Mossad, 2009; Stech et al., 2010). During the pandemic, K was substituted with R. We speculated that the NP-470 R to K substitution in influenza A (H1N1) 2009 may be one of the contributing factors that account for the higher pathogenicity and mammalian adaptation. The roles of NP-K470R in influenza H1N1 2009 in viral life cycle, pathogenicity and mammal transmissibility will be examined in our further work. In this study, we revealed that K470R substitution in H5N1 facilitated viral replication and pathogenicity in vitro and in vivo. Our findings may be helpful for the surveillance of IAVs.

TABLE 3
www.frontiersin.org

TABLE 3. Conservation of NP residues 91, 198, 227, 229, and 470 in influenza A viruses.

TABLE 4
www.frontiersin.org

TABLE 4. NP residues 470 in avian and human influenza A viruses.

The influenza virus RNA polymerase is an important virulence determinant (Watanabe et al., 2013). In the current study, the RNP activity of K227R and K229R was higher compared with WT, but the mutants did not show increased the virus replication or virulence in mammals. These results suggested that in addition to RNP activity, other factors were also associated with the pathogenicity of influenza A, such as viral entry, trafficking to the cell surface, viral assembly and host immune response (Neumann and Kawaoka, 2006; Keleta et al., 2008; Shapira et al., 2009; Bottcher-Friebertshauser et al., 2014). Recent studies demonstrated that G158 mutation of the HA protein enhanced viral productivity and induced stronger host immune, and thereby contributed to the high pathogenicity of H5N1 virus in mice (Zhao et al., 2017). Mutations in the M1 protein also contributed to the virulence of H5N1 influenza viruses in mice (Fan et al., 2009). Though there are many important data in this area, the mechanism of the virulent determinants of influenza A viruses needs to be further explored.

Nucleoprotein can shutter between the nucleus and the cytoplasm because of its nuclear localization signals (NLSs) and nuclear export signals (NESs). Three NESs have been found in NP, NES1 (E24IRASVGKMIDGIGRFYIQMCTELKL49) and NES2 (V183KGVGTMVMELIRMI197), which do not rely on CRM1 for export, and NES3 (P248GNAEFEDLIFL ARSALILRGSVAHKS274), which does (Yu et al., 2012). The three NESs are closely linked to nuclear export of vRNP, and mutation of the NP NESs inhibits the replication of the virus. In addition, there are three NLSs in NP. An unconventional NLS (M1ASQGTKRSYEQM13) is at the N-terminus of the NP. The second NLS (K198RGINDRNFWRGFNGRRTR216) is located in the central region of NP. The third NLS is located between amino acids 320 and 400 (Neumann et al., 1997). The K198R mutation, located in the second NLS, affected the NP nuclear location; polymerase activity was lower, and the recombinant virus K198R was not viable.

K91 is located in the body domain of NP. NP was shown to interact with PB2 through residues 1-161 and 255-465. Viruses with single point mutations at K91R could not be rescued, indicating that the residue is essential for the survival of the virus, probably because of its important roles in PB2 binding.

NP, a major component of vRNP, plays a role in RNA synthesis in the nucleus. Although PB1, PB2 and PA proteins are sufficient to synthesize short RNAs, NP is necessary for synthesis of longer RNA (Honda et al., 1988). NP is known to interact with several viral proteins, such as PB2, PB1, NS1, M1 and NP itself. NP is also believed to be a key adaptor for virus and host-cell interaction. Many efforts have been made to identify host proteins associated with NP (Generous et al., 2014; Sun et al., 2015). The host protein cluster (CLU) was identified as a novel interacting partner with influenza A viral NP during virus infection. Overexpression of CLU attenuated IAV replication in mammalian cells (Tripathi et al., 2013). NP interacted with F actin and caused cytoplasmic retention of NP late in infection (Avalos et al., 1997). The TREX complex adaptor protein Aly/REF interacted with NP, and siRNA knockdown of Aly/REF had little effect on export of NP and no significant reduction in viral titer (Balasubramaniam et al., 2013). Using Gal4-based yeast two-hybrid(Y2H) assay, phosphoribo-sylaminoimidazole succinocarboxamide synthetase (PAICS) was shown to interact with NP. PAICS was found to be significantly up-regulated during an IAV infection (Generous et al., 2014). The study of the interaction between host and NP can provide new ideas for the development of clinical drug target and the prevention of the disease (Kao et al., 2010).

Because influenza viruses have a segmented RNA genome, exchange of segments between different viruses and mutation can occur when different influenza viruses infect the same cell (Wu et al., 2013; Gerber et al., 2014). In our previous studies, we showed that diverse wild bird H9N9 influenza virus genes elevated the virulence of a H7N9 virus (Su et al., 2015). Recently, cases of human infections with newly reasserted influenza A (H7N9) virus have been continuously reported in China (Sha et al., 2016; Shen and Lu, 2017). The emergence and transmission of novel avian IAVs to humans are complicated, which can be due to a series of events, including the reassortment of the virus genome and accumulation of adaptive mutations.

In this study, we analyzed the effect of NP mutants on the virulence of the H5N1 in mammals. Of these, K91R and K198R disrupted the NP cellular location and decreased the polymerase activity. Compared with the rWT, the rK229R mutant exhibited an attenuated virulence, while the rK227R showed the same virulence as that of the rWT. Importantly, the K470R mutant enhanced viral replication and pathogenicity in mammals. The identification of these substitutions may contribute to our understanding of the biology of IAVs. In addition, it may also facilitate the prediction and controlling of future IAV epidemics and pandemics.

Author Contributions

HH and LC designed this experiment, and wrote this manuscript; LC, JL, ML, JH, XZ, GY, and GM performed the experiments; LC, CW, HL and NZ participated in the data analysis; CW and ML provided some useful suggestion on the project and manuscript revising.

Funding

This study was funded by grants from the National Natural Science Foundation of China (31602049, 31301710), and the State Forestry Administration of the People’s Republic of China, Chinese Academy of Sciences (CZBZX-1).

Conflict of Interest Statement

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.

Acknowledgment

The authors are grateful to Wen Su, Qingxun Zhang, and Lv Lu for helpful discussion.

Footnotes

  1. ^http://www.who.int
  2. ^http://platform.gisaid.org

References

Avalos, R. T., Yu, Z., and Nayak, D. P. (1997). Association of influenza virus NP and M1 proteins with cellular cytoskeletal elements in influenza virus-infected cells. J. Virol. 71, 2947–2958.

PubMed Abstract | Google Scholar

Balasubramaniam, V. R., Hong Wai, T., Ario Tejo, B., Omar, A. R., and Syed Hassan, S. (2013). Highly pathogenic avian influenza virus nucleoprotein interacts with TREX complex adaptor protein Aly/REF. PLoS ONE 8:e72429. doi: 10.1371/journal.pone.0072429

PubMed Abstract | CrossRef Full Text | Google Scholar

Bottcher-Friebertshauser, E., Garten, W., Matrosovich, M., and Klenk, H. D. (2014). The hemagglutinin: a determinant of pathogenicity. Curr. Top. Microbiol. Immunol. 385, 3–34. doi: 10.1007/82_2014_384

CrossRef Full Text | Google Scholar

Bussey, K. A., Bousse, T. L., Desmet, E. A., Kim, B., and Takimoto, T. (2010). PB2 residue 271 plays a key role in enhanced polymerase activity of influenza A viruses in mammalian host cells. J. Virol. 84, 4395–4406. doi: 10.1128/JVI.02642-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Lee, K. H., Steinhauer, D. A., Stevens, D. J., Skehel, J. J., and Wiley, D. C. (1998). Structure of the hemagglutinin precursor cleavage site, a determinant of influenza pathogenicity and the origin of the labile conformation. Cell 95, 409–417. doi: 10.1016/S0092-8674(00)81771-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W. S., Calvo, P. A., Malide, D., Gibbs, J., Schubert, U., Bacik, I., et al. (2001). A novel influenza A virus mitochondrial protein that induces cell death. Nat. Med. 7, 1306–1312. doi: 10.1038/nm1201-1306

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Shesheny, R., Mostafa, A., Kandeil, A., Mahmoud, S. H., Bagato, O., Naguib, A., et al. (2017). Biological characterization of highly pathogenic avian influenza H5N1 viruses that infected humans in Egypt in 2014-2015. Arch. Virol. 162, 687–700. doi: 10.1007/s00705-016-3137-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, S., Deng, G., Song, J., Tian, G., Suo, Y., Jiang, Y., et al. (2009). Two amino acid residues in the matrix protein M1 contribute to the virulence difference of H5N1 avian influenza viruses in mice. Virology 384, 28–32. doi: 10.1016/j.virol.2008.11.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Generous, A., Thorson, M., Barcus, J., Jacher, J., Busch, M., and Sleister, H. (2014). Identification of putative interactions between swine and human influenza A virus nucleoprotein and human host proteins. Virol. J. 11:228. doi: 10.1186/s12985-014-0228-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerber, M., Isel, C., Moules, V., and Marquet, R. (2014). Selective packaging of the influenza A genome and consequences for genetic reassortment. Trends Microbiol. 22, 446–455. doi: 10.1016/j.tim.2014.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Hale, B. G., Knebel, A., Botting, C. H., Galloway, C. S., Precious, B. L., Jackson, D., et al. (2009). CDK/ERK-mediated phosphorylation of the human influenza A virus NS1 protein at threonine-215. Virology 383, 6–11. doi: 10.1016/j.virol.2008.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Q., Chang, C., Li, L., Klenk, C., Cheng, J., Chen, Y., et al. (2014). Sumoylation of influenza A virus nucleoprotein is essential for intracellular trafficking and virus growth. J. Virol. 88, 9379–9390. doi: 10.1128/JVI.00509-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffmann, E., Neumann, G., Kawaoka, Y., Hobom, G., and Webster, R. G. (2000). A DNA transfection system for generation of influenza A virus from eight plasmids. Proc. Natl. Acad. Sci. U.S.A. 97, 6108–6113. doi: 10.1073/pnas.100133697

PubMed Abstract | CrossRef Full Text | Google Scholar

Honda, A., Ueda, K., Nagata, K., and Ishihama, A. (1988). RNA polymerase of influenza virus: role of NP in RNA chain elongation. J. Biochem. 104, 1021–1026. doi: 10.1093/oxfordjournals.jbchem.a122569

PubMed Abstract | CrossRef Full Text | Google Scholar

Jagger, B. W., Wise, H. M., Kash, J. C., Walters, K. A., Wills, N. M., Xiao, Y. L., et al. (2012). An overlapping protein-coding region in influenza A virus segment 3 modulates the host response. Science 337, 199–204. doi: 10.1126/science.1222213

PubMed Abstract | CrossRef Full Text | Google Scholar

Kao, R. Y., Yang, D., Lau, L. S., Tsui, W. H., Hu, L., Dai, J., et al. (2010). Identification of influenza A nucleoprotein as an antiviral target. Nat. Biotechnol. 28, 600–605. doi: 10.1038/nbt.1638

PubMed Abstract | CrossRef Full Text | Google Scholar

Keleta, L., Ibricevic, A., Bovin, N. V., Brody, S. L., and Brown, E. G. (2008). Experimental evolution of human influenza virus H3 hemagglutinin in the mouse lung identifies adaptive regions in HA1 and HA2. J. Virol. 82, 11599–11608. doi: 10.1128/JVI.01393-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Kistner, O., Muller, K., and Scholtissek, C. (1989). Differential phosphorylation of the nucleoprotein of influenza A viruses. J. Gen. Virol. 70, 2421–2431. doi: 10.1099/0022-1317-70-9-2421

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Watanabe, T., Hatta, M., Watanabe, S., Nanbo, A., Ozawa, M., et al. (2009). Mutational analysis of conserved amino acids in the influenza A virus nucleoprotein. J. Virol. 83, 4153–4162. doi: 10.1128/JVI.02642-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, T. L., Wu, C. Y., Su, W. C., Jeng, K. S., and Lai, M. M. (2010). Ubiquitination and deubiquitination of NP protein regulates influenza A virus RNA replication. EMBO J. 29, 3879–3890. doi: 10.1038/emboj.2010.250

PubMed Abstract | CrossRef Full Text | Google Scholar

Ly, S., Horwood, P., Chan, M., Rith, S., Sorn, S., Oeung, K., et al. (2017). Seroprevalence and transmission of human influenza A(H5N1) virus before and after virus reassortment, Cambodia, 2006-2014. Emerg. Infect. Dis. 23, 300–303. doi: 10.3201/eid2302.161232

PubMed Abstract | CrossRef Full Text | Google Scholar

Muramoto, Y., Noda, T., Kawakami, E., Akkina, R., and Kawaoka, Y. (2013). Identification of novel influenza A virus proteins translated from PA mRNA. J. Virol. 87, 2455–2462. doi: 10.1128/JVI.02656-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, B. R., and Webster, R. G. (1996). Orthomyxoviruses. Fields Virol. 3, 1397–1445.

Google Scholar

Neumann, G., Castrucci, M. R., and Kawaoka, Y. (1997). Nuclear import and export of influenza virus nucleoprotein. J. Virol. 71, 9690–9700. doi: 10.3201/eid1206.051336

CrossRef Full Text | Google Scholar

Neumann, G., and Kawaoka, Y. (2006). Host range restriction and pathogenicity in the context of influenza pandemic. Emerg. Infect. Dis. 12, 881–886.

Google Scholar

Ng, A. K., Wang, J. H., and Shaw, P. C. (2009). Structure and sequence analysis of influenza A virus nucleoprotein. Sci. China C Life Sci. 52, 439–449. doi: 10.1007/s11427-009-0064-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ng, A. K., Zhang, H., Tan, K., Li, Z., Liu, J. H., Chan, P. K., et al. (2008). Structure of the influenza virus A H5N1 nucleoprotein: implications for RNA binding, oligomerization, and vaccine design. FASEB J. 22, 3638–3647. doi: 10.1096/fj.08-112110

PubMed Abstract | CrossRef Full Text | Google Scholar

Portela, A., and Digard, P. (2002). The influenza virus nucleoprotein: a multifunctional RNA-binding protein pivotal to virus replication. J. Gen. Virol. 83, 723–734. doi: 10.1099/0022-1317-83-4-723

PubMed Abstract | CrossRef Full Text | Google Scholar

Scalera, N. M., and Mossad, S. B. (2009). The first pandemic of the 21st century: a review of the 2009 pandemic variant influenza A (H1N1) virus. Postgrad. Med. 121, 43–47. doi: 10.3810/pgm.2009.09.2051

PubMed Abstract | CrossRef Full Text | Google Scholar

Selman, M., Dankar, S. K., Forbes, N. E., Jia, J. J., and Brown, E. G. (2012). Adaptive mutation in influenza A virus non-structural gene is linked to host switching and induces a novel protein by alternative splicing. Emerg. Microbes Infect. 1, e42. doi: 10.1038/emi.2012.38

PubMed Abstract | CrossRef Full Text | Google Scholar

Sha, J. P., Chen, X. W., Ren, Y. J., Chen, H. J., Wu, Z. Q., Ying, D., et al. (2016). Differences in the epidemiology and virology of mild, severe and fatal human infections with avian influenza A (H7N9) virus. Arch. Virol. 161, 1239–1259. doi: 10.1007/s00705-016-2781-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Shapira, S. D., Gat-Viks, I., Shum, B. O., Dricot, A., de Grace, M. M., Wu, L., et al. (2009). A physical and regulatory map of host-influenza interactions reveals pathways in H1N1 infection. Cell 139, 1255–1267. doi: 10.1016/j.cell.2009.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y. Z., and Lu, H. Z. (2017). Global concern regarding the fifth epidemic of human infection with avian influenza A (H7N9) virus in China. Biosci. Trends 11, 120–121. doi: 10.5582/bst.2017.01040

PubMed Abstract | CrossRef Full Text | Google Scholar

Stech, J., Beer, M., Vahlenkamp, T., and Harder, T. (2010). The pandemic influenza virus H1N1/2009: a review of the molecular biology, phylogeny, history of reassortments, and parameters of host switching. Bundesgesundheitsblatt Gesundheitsforschung Gesundheitsschutz 53, 1231–1237. doi: 10.1007/s00103-010-1166-0

PubMed Abstract | CrossRef Full Text

Su, W., Wang, C., Luo, J., Zhao, Y., Wu, Y., Chen, L., et al. (2015). Testing the effect of internal genes derived from a wild-bird-origin H9N2 influenza A virus on the pathogenicity of an A/H7N9 virus. Cell Rep. 12, 1831–1841. doi: 10.1016/j.celrep.2015.08.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, N., Sun, W., Li, S., Yang, J., Yang, L., Quan, G., et al. (2015). Proteomics analysis of cellular proteins co-immunoprecipitated with nucleoprotein of influenza A virus (H7N9). Int. J. Mol. Sci. 16, 25982–25998. doi: 10.3390/ijms161125934

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarendeau, F., Crepin, T., Guilligay, D., Ruigrok, R. W. H., Cusack, S., and Hart, D. J. (2008). Host determinant residue lysine 627 lies on the surface of a discrete, folded domain of influenza virus polymerase PB2 subunit. PLoS Pathog. 4:e1000136. doi: 10.1371/journal.ppat.1000136

PubMed Abstract | CrossRef Full Text | Google Scholar

Tripathi, S., Batra, J., Cao, W., Sharma, K., Patel, J. R., Ranjan, P., et al. (2013). Influenza A virus nucleoprotein induces apoptosis in human airway epithelial cells: implications of a novel interaction between nucleoprotein and host protein Clusterin. Cell Death Dis. 4:e562. doi: 10.1038/cddis.2013.89

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Zhao, Z., Bi, Y., Sun, L., Liu, X., and Liu, W. (2013). Tyrosine 132 phosphorylation of influenza A virus M1 protein is crucial for virus replication by controlling the nuclear import of M1. J. Virol. 87, 6182–6191. doi: 10.1128/JVI.03024-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, T., Tisoncik-Go, J., Tchitchek, N., Watanabe, S., Benecke, A. G., Katze, M. G., et al. (2013). 1918 Influenza virus hemagglutinin (HA) and the viral RNA polymerase complex enhance viral pathogenicity, but only HA induces aberrant host responses in mice. J. Virol. 87, 5239–5254. doi: 10.1128/JVI.02753-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Wise, H. M., Hutchinson, E. C., Jagger, B. W., Stuart, A. D., Kang, Z. H., Robb, N., et al. (2012). Identification of a novel splice variant form of the influenza A virus M2 ion channel with an antigenically distinct ectodomain. PLoS Pathog. 8:e1002998. doi: 10.1371/journal.ppat.1002998

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, A., Su, C., Wang, D., Peng, Y., Liu, M., Hua, S., et al. (2013). Sequential reassortments underlie diverse influenza H7N9 genotypes in China. Cell Host Microbe 14, 446–452. doi: 10.1016/j.chom.2013.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, C., Hu, W. B., Xu, K., He, Y. X., Wang, T. Y., Chen, Z., et al. (2012). Amino acids 473V and 598P of PB1 from an avian-origin influenza A virus contribute to polymerase activity, especially in mammalian cells. J. Gen. Virol. 93, 531–540. doi: 10.1099/vir.0.036434-0

PubMed Abstract | CrossRef Full Text | Google Scholar

York, A., and Fodor, E. (2013). Biogenesis, assembly, and export of viral messenger ribonucleoproteins in the influenza A virus infected cell. RNA Biol. 10, 1274–1282. doi: 10.4161/rna.25356

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, M., Liu, X., Cao, S., Zhao, Z., Zhang, K., Xie, Q., et al. (2012). Identification and characterization of three novel nuclear export signals in the influenza A virus nucleoprotein. J. Virol. 86, 4970–4980. doi: 10.1128/JVI.06159-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, D., Liang, L., Wang, S., Nakao, T., Li, Y., Liu, L., et al. (2017). Glycosylation of the hemagglutinin protein of H5N1 influenza virus increases its virulence in mice by exacerbating the host immune response. J. Virol. 91:e2215-16. doi: 10.1128/JVI.02215-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H. J., Chu, H., Zhao, X. Y., Shuai, H. P., Wong, B. H. Y., Wen, L., et al. (2016). Novel residues in the PA protein of avian influenza H7N7 virus affect virulence in mammalian hosts. Virology 498, 1–8. doi: 10.1016/j.virol.2016.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, W., Li, J., Wang, S., Cao, S., Jiang, J., Chen, C., et al. (2015). Phosphorylation controls the nuclear-cytoplasmic shuttling of influenza A virus nucleoprotein. J. Virol. 89, 5822–5834. doi: 10.1128/JVI.00015-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: highly pathogenic avian influenza, H5N1, nucleoprotein, mutation, virulence, mammalian cells

Citation: Chen L, Wang C, Luo J, Li M, Liu H, Zhao N, Huang J, Zhu X, Ma G, Yuan G and He H (2017) Amino Acid Substitution K470R in the Nucleoprotein Increases the Virulence of H5N1 Influenza A Virus in Mammals. Front. Microbiol. 8:1308. doi: 10.3389/fmicb.2017.01308

Received: 07 April 2017; Accepted: 28 June 2017;
Published: 11 July 2017.

Edited by:

Sunil Kumar Lal, Monash University, Australia

Reviewed by:

Amelia Nieto, Consejo Superior de Investigaciones Científicas (CSIC), Spain
Haider Abdul-Lateef Mousa, University of Basrah, Iraq
Shashank Tripathi, Icahn School of Medicine at Mount Sinai, United States

Copyright © 2017 Chen, Wang, Luo, Li, Liu, Zhao, Huang, Zhu, Ma, Yuan and He. 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) or licensor 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: Hongxuan He, hehx@ioz.ac.cn